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Abstract 



The vanishing moment method was introduced by the authors in |37j as a reh- 
able methodology for computing viscosity solutions of fully nonlinear second order 
partial differential equations (PDEs), in particular, using Galerkin-type numeri- 
cal methods such as finite element methods, spectral methods, and discontinuous 
Galerkin methods, a task which has not been practicable in the past. The crux of 
the vanishing moment method is the simple idea of approximating a fully nonlinear 
second order PDE by a family (parametrized by a small parameter e) of quasilinear 
higher order (in particular, fourth order) PDEs. The primary objectives of this 
book are to present a detailed convergent analysis for the method in the radial 
symmetric case and to carry out a comprehensive finite element numerical analy- 
sis for the vanishing moment equations (i.e., the regularized fourth order PDEs). 
Abstract methodological and convergence analysis frameworks of conforming finite 
element methods and mixed finite element methods are first developed for fully 
nonlinear second order PDEs in general settings. The abstract frameworks are 
then applied to three prototypical nonlinear equations, namely, the Monge- Ampere 
equation, the equation of prescribed Gauss curvature, and the infinity-Laplacian 
equation. Numerical experiments are also presented for each problem to validate 
the theoretical error estimate results and to gauge the efficiency of the proposed 
numerical methods and the vanishing moment methodology. 
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CHAPTER 1 



Prelude 



1.1. Introduction 

Fully nonlinear partial differential equations (PDEs) are those equations which 
are nonlinear in the highest order derivative(s) of the unknown function(s). In the 
case of the second order equations, the general form of fully nonlinear PDEs is given 

by 

(1.1) F{D^u,\/u,u,x) ^0, 

where D^u{x) and Vu{x) denote respectively the Hessian and the gradient of u 
at a: g 17 C R". Here, F is assumed to be a nonlinear function in at least one 
of its entries of D^u. Fully nonlinear PDEs arise from many scientific and engi- 
neering fields including differential geometry, optimal control, mass transportation, 
geostrophic fluid, meteorology, and general relativity (cf. |18L I19L 1421 14H [58j and 
the references therein). 

Examples of such equations include (cf . |42j ) 

• The Monge-Ampere equation 

(1.2) det{D^u) = /. 

• The equation of prescribed Gauss curvature 

(1.3) dct{D^u)=IC(l + \\/u\^)'^. 

• The Bellman equation 

(1.4) inf(Leu-/e) = 0. 

Here, det{D'^u{x)) denotes the determinant of the Hessian D^u at x, and {Lg} 
denotes a family of second order linear differential operators. 



Because of the full nonlinearity in (1.1 1, the standard weak solution theory 



based on the integration by parts approach does not work and other notions of weak 
solutions must be sought. Progress has been made in the latter half of the twentieth 
century concerning this issue after the introduction of viscosity solutions. In 1983, 
Crandall and Lions [24j introduced the notion of viscosity solutions and used the 
vanishing viscosity method to show existence of a solution for the Hamilton- Jacobi 
equation: 

(1.5) ut + H{S/u,u,x) ^0 inR"x(0,oo). 

The vanishing viscosity method approximates the Hamilton- Jacobi equation by 
the following regularized, second order quasilinear PDE: 

(1.6) < - eAu- -I- i^(V^t^ ^t^ x) = in R" x (0, cx)). 
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It was shown that 24j there exists a unique solution to the regularized 
Cauchy problem that converges locally and uniformly to a continuous function u 
which is defined to be a viscosity solution of the Hamilton- Jacobi equation (1.5 1. 
However, to establish uniqueness, the following intrinsic definition of viscosity so- 
lutions was also proposed 



in VL, 
on dVl. 



Definition 1.1. Let H : R" x Rx O — > R and g : dVl — > R be continuous functions, 
and consider the following problem: 

(1.7) H{Vu,u,x) = Q 

(1.8) u = g 



(i) u e C°(r2) is called a viscosity suhsolution of (1.7)-(1.8) if uj^^ = (7, and 



for every function i^ix) such that u — ip has a local maximum at Xq G 
there holds 



(ii) u e C*'(ri) is called a viscosity supersolution of (1.7|-(1.8| if u\ 



dn 



9, 



and for every function (p{x) such that u — (p has a local minimum at 
xq E n, there holds 

H{V(p{xq),u{xo),xo) > 0. 



(iii) u £ C"(r2) is called a OTscosiiy so/ittion of (1.7)-(1.8| if it is both a viscosity 
subsolution and supersolution. 

Clearly, the above definition is not variational, as it is based on a "differentiation 
by parts" approach (a terminology introduced in \24\ I25j). In addition, the word 
"viscosity" loses its original meaning in the definition. However, it was shown 
|24L 125] that every viscosity solution constructed by the vanishing viscosity method 
is an intrinsic viscosity solution (i.e., a solution that satisfies Definition |l.l[ ). Besides 
addressing the uniqueness issue, another reason to favor the intrinsic differentiation 
by parts definition is that the definition and the notion of viscosity solutions can 
be readily extended to fully nonlinear second order PDEs as follows (cf. 18 ): 

R and g : dfl — >■ R be continuous 



in J7, 
on dil. 



Definition 1.2. Let F : R"^" x R" x R x 
functions, and consider the following problem: 

(1.9) F{D'^u,Wu,u,x) ^0 

(1.10) u^g 

(i) u G C^{n) is called a viscosity subsolution of ( 1.9 1-( 1.10) if u 

and for every function (p{x) such that u ~ ip has a local maximum at 
xq G^l, there holds 

u{xq),X(^) < 



idn 



9, 



(ii) u G C'^{Q,) is called a viscosity supersolution of ( 1.9 )-( 1.10) if 



and for every function ip{x) such that u 
xo £ ^, there holds 

F{D^^{xo),^v{xo), u{xq),Xq) > 



I an 



= 9, 



Lp has a local minimum at 



(iii) u E C'^{il) is called a viscosity solution of ( 1.9 1-( 1.10) if it is both a 
viscosity subsolution and supersolution. 
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Figure 1. A Geometric interpretation of viscosity solutions 



Remark 1.3. Without loss of generality, we may assume that u{xo) = (p(xo) 
whenever u — tp achieves a local maximum or local minimum at Xq € in Definition 
|1.2[ Therefore, in an informal setting, u is a viscosity solution if for every smooth 
function ip that "touches" the graph of u from above at xq (see Figure [T]) there 
holds 

F{D^(p{xo),W(p{xo),(p{xQ),xo) < 0, 
and if (p "touches" the graph of u from below at xq, then 
F{D'^Lp{xQ),W<fi{xo),ip{xo),XQ) > 0. 

In case of the fully nonlinear first order PDEs, tremendous progress has been 
made in the past three decades in terms of PDF analysis and numerical methods. 
A profound viscosity solution theory has been established (cf. |24l 1251 1261 I41| l 
and a wealth of efficient and robust numerical methods and algorithms have been 
developed and implemented (cf. [m [T4l[23l[27l[56l[66l[67l[68l[72l [79l [80 ] ). 
However, in the case of fully nonlinear second order PDFs, the situation is strikingly 
different. On the one hand, there have been enormous advances in PDF analysis 
in the past two decades after the introduction of the notion of viscosity solutions 
by M. Crandall and P. L. Lions in 1983 (cf. [HI [26l |44] ) . On the other hand, in 
contrast to the success of the PDF analysis, numerical solutions for general fully 
nonlinear second order PDFs is a relatively untouched area. 

There are several reasons for this lack of progress in numerical methods. First, 
the most obvious difficulty is the full nonlinearity in the equation. Second, solu- 
tions to fully nonlinear second order equations are often only unique in a certain 
class of functions, and this conditional uniqueness is very difficult to handle nu- 
merically. Lastly and most importantly, it is extremely difficult (if all possible) to 
mimic the differentiation by parts approach at the discrete level. As a consequence, 
there is little hope to develop a discrete viscosity solution theory. Furthermore, it 
is impossible to directly compute viscosity solutions using Galerkin-type numerical 
methods including finite element methods, spectral Galerkin methods, and discon- 
tinuous Galerkin methods, since they are all based on variational formulations of 
PDFs. In fact, this is clear from the definition of viscosity solutions, which is not 
based on the traditional integration by parts approach, but rather is defined by the 
differentiation by parts approach. 



4 



To explain the above points, consider the Dirichlet problem for the Monge- 
Ampere equation as an example: 

(1.11) det{D'^u) = / inn, 

(1.12) u^g ondn, 

which corresponds to F{D'^u,Vu,u,x) = f(x) — det(D^u). It is known that for a 
non-strictly convex domain f2, the above problem does not have classical solutions 
in general even if f,g, and dfl are smooth [42 . Results of A. D. Aleksandrov 
state that the Dirichlet problem with / > has a unique generalized solution 
(which is also the viscosity solution) in the class of convex functions [2l I44j . The 
reason to restrict the admissible set to be the set of convex functions is that the 
Monge- Ampere equation is only elliptic in that set |42L I44j . It should be noted 



that in general, the Dirichlet problem ( 1.11 )-( 1.12 ) may have other nonconvex 
solutions even when / > 0. It is easy to see that if one discretizes ( |1.11[ ) directly 
using a standard finite difference method, not only would the resulting algebraic 
system be difficult to solve, one immediately loses control on which solution the 
numerical scheme approximates - and this is assuming that the nonlinear discrete 
problem has solutions! Furthermore, the situation is even worse if one tries to 
formulate a Galerkin-type numerical method because there is not a variational or 
weak formulation in which to start. 

Nevertheless, a few recent numerical attempts and results have been known in 
the literature. In j 65| Oliker and Prussner proposed a finite difference scheme for 
computing Aleksandrov measure induced by D^u (and obtained the solution u of 



(1.2) as a by-product) in two dimensions. The scheme is extremely geometric and 
difficult to generalize to other fully nonlinear second order PDEs. In [6] Barles 
and Souganidis showed that any monotone, stable, and consistent finite difference 
scheme converges to the viscosity solution provided that there exists a comparison 
principle for the limiting equation. Their result provides a guideline for construct- 
ing convergent finite difference methods, but they did not address how to construct 
such a scheme. In [4], Baginski and Whitaker proposed a finite difference scheme 



for the equation of prescribed Gauss curvature (1.3) in two dimensions by mimick- 
ing the unique continuation method (used to prove existence of the PDE) at the 
discrete level. The method becomes very unstable when the homotopy is domi- 
nated by the fully nonlinear equation. Oberman |64j constructed a wide stencil 
finite difference scheme for fully nonlinear elliptic PDEs which can be written as 
functions of eigenvalues of the Hessian matrix and proved that the scheme satisfies 
the convergence criterion established by Barles and Souganidis in [6]. In a series 
of papers [291 1301 131] Dean and Glowinski proposed an augmented Lagrange mul- 



tiplier method and a least squares method for problem (1.2) and Pucci's equation 
(cf. (18], I42j ) in two dimensions by treating the nonlinear PDEs as a constraint 
and using a variational criterion to select a particular solution. However, as the 
admissible set is contained in H^{n), it could become empty if all solutions of the 
underlying fully nonlinear PDE are not differentiable. Finally, Bohmer [15j re- 
cently introduced a projection method using finite elements for approximating 
classical solutions of a certain class of fully nonlinear second order elliptic PDEs. 
However, the issue of how to reliably compute a selected solution (the resulting 
discrete problem often has multiple solutions) was not addressed and still remains 
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an open question. Numerical experiments were reported in \65\ IH I64|, 129^ I30|,l31j. 

however, convergence analysis was not addressed except in |64| . 

In addition, we like to remark that there is a considerable amount of literature 
available on using finite difference methods to approximate viscosity solutions of 
fully nonlinear second order Bellman-type PDEs arising from stochastic optimal 
control (cf. [6|, [71 150|, I53j ). However, due to the special structure of Bellman- type 
PDEs, the approach used and the methods proposed in those papers could not be 
extended to other types of fully nonlinear second order PDEs since the construction 
of those methods critically relies on the linearity of the operators Lg . 

The first goal of this book is to present a general framework for the vanishing 
moment method and the notion of moment solutions for fully nonlinear second 
order PDEs. The vanishing moment method is very much in the same spirit of the 
vanishing viscosity method introduced in 24 , and the notion of moment solutions 
for fully nonlinear second order PDEs is a natural extension of the (original) notion 
of viscosity solutions for fully nonlinear first order PDEs. This methodology was 
first introduced by the authors in [37] as a reliable way for computing viscosity 
solutions of fully nonlinear second order PDEs, in particular, using Galerkin-type 
numerical methods. The crux of this new method is to approximate a fully nonlinear 
second order PDE by a family of quasilinear fourth order PDEs. The limit of the 
solutions of the fourth order PDEs (if it exists) is defined as a moment solution 
of the original fully nonlinear second order PDE. As moment solutions are defined 
constructively, they can be readily computed by existing numerical methods. In 
the case of Monge- Ampere-type equations, extensive numerical experiments in [371 
1381 1391 162] suggest that the moment solution coincides with the viscosity solution 
as long as the latter exists. In this book, we shall present a detailed convergence 
theory for the vanishing moment method in the radial symmetric case. This then 
provides a theoretical foundation for the method and for the numerical results of 
(an [38, 39, 62 . 

The second goal of this book, which is the bulk of the book's content, is to carry 
out a comprehensive finite element numerical analysis for the vanishing moment 
method. Two abstract frameworks are developed for this purpose in a general 
setting. The first framework concerns conforming finite element approximations 
of the vanishing moment equations (i.e., the regularized fourth order equations). 
The second framework develops (Herman-Miyoshi) mixed finite element methods 
for the vanishing moment equations. Each of these two frameworks consists of the 
formulation of the respective numerical methods, proving existence and uniqueness 
of numerical solutions, and deriving error estimates for the numerical solutions. Due 
to the strong nonlinearity of the PDEs, the standard numerical analysis techniques 
for finite element methods do not work here. To overcome the difficulty, we combine 
a fixed point argument with a linearization technique. After having completed both 
abstract frameworks, we apply them to three prototypical nonlinear equations, 
namely, the Monge-Ampere equation, the equation of prescribed Gauss curvature, 
and the infinity-Laplacian equation. The three equations are chosen because they 
present three different and interesting scenarios, that is, their linearizations are 
respectively coercive, indefinite, and degenerate. It is shown that our abstract 
frameworks are rich enough to cover all three scenarios. 
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The remainder of the book is organized as follows. Chapter [2] represents the 
formulation of the vanishing moment method and its informal insights. The mate- 
rial of this chapter has a large overlap with that of 37". Chapter |3] is devoted to the 
convergence analysis of the vanishing moment method for the Monge- Ampere equa- 
tion in the radial symmetric case. The main tasks of the chapter are to analyze the 
vanishing moment equations and to derive uniform (in e) estimates for its solutions. 
The chapter also contains a convergence rate estimate result for the regularized so- 
lutions in the case that the viscosity solution of the Monge- Ampere equation belongs 
to W'^'°°{fl) n H^{n). Chapter [4] and [5] develop, respectively, the abstract frame- 
works for the two types of finite element (i.e., conforming and mixed finite element) 
approximations of the vanishing moment equations under some structure assump- 
tions on the nonlinear differential operator F. Chapter[6]presents applications of the 
abstract frameworks of Chapter [4] and [5] to three prototypical nonlinear equations: 
the Monge- Ampere equation, the equation of prescribed Gauss curvature, and the 
infinity-Laplacian equation. For each equation, we formulate its vanishing moment 
approximations, subsequent finite element and mixed finite element methods, and 
obtain their error estimates by fitting the equation into the abstract frameworks. 
For the Monge- Ampere equation, besides some slight improvements, we essentially 
recover the early results reported in [38, 39 . On the other hand, the results for 
the equation of prescribed Gauss curvature and the infinity-Laplacian equation are 
new. In fact, to the best of our knowledge, no comparable results are known in the 
literature. Numerical experiments are also presented for each problem to validate 
the theoretical (error estimate) results, and to gauge the efficiency of the proposed 
numerical methods and the vanishing moment methodology. Finally, we end the 
book with a few concluding remarks in Chapter [7j 



1.2. Preliminaries 

Standard space notation is adopted in this book, we refer the reader to [13|, 142^ 
I22j for their exact definitions. In addition, J7 denotes a bounded convex domain 
in R". (•, •) and {■,-)dn denote the L^-inner products on n and on dil, respec- 
tively. The unlabeled constant C is used to denote generic e- and /i— independent 
positive constants that may take on different values at different occurrences, where 
as labeled constants denote e-dependent (but /i-independent) constants. Further- 
more all constants, labeled and unlabeled, are chapter-independent unless otherwise 
specified. 

Throughout this book we assume that 

F : R"^" xR"xRxf2 — ^R 

is a differentiable function in all its arguments. For a given (small) constant e > 0, 
we define 



G,{r,p,z,x) -.^ eAtT{r) + F{r,p,z,x) Vr e R"^", p e R", z e R, xefl. 



7 



For a given scalar function v and an n x n matrix- valued function /z = [fJ-ij] we 
set 

Fp[r,p,z,x]{v) 
F^[r,p, z,x]{v) 

F'[r,p,z,x]{fi,v) := Fr[r,p, z,x]{fi) + Fp[r,p, z,x]{v) + F^[r,p, z,x]{v), 
G'^[r,p, z, x]{^, v) :— eAtr(^) + F'[r,p, z, x]{fi, v). 

We also define, with a slight abuse of notation, for a scalar function w and an 
n X n tensor function k — [nij], the following short-hand notation, which will be 
extensively used when developing mixed finite element methods in Chapter [5j 

(1.13) F{k,w) -.^ F{k,\/w,w,x), 

Fr[K,w]{^J.) -.^ Fr[K,\'w,W,x]{p), 

Fp[K,w]{v) := Fp[K,Vw,w,x]{v), 
Fz[k, w]{v) :— Fz[k, Vw, w, x]{v), 

:= Fr[K,w]{p) + Fp[K,w]{v) + F^[k,w]{v), 
Ge{K, w) eAtr(K) + F{k, w), 
Gg[K, w] {p, v) :~ eAtr(/i) + F'[k, w] (/i, v). 

For notation used in Chapter [4j we overload the operators F,G^,F', and 
once again and define the additional short-hand notation: 

(1.14) F{w) := F{D'^w,w), 
Fr[w]{v) := Fr[D^w,w]iD^v), 
Fp[w]iv) ■.^Fp[D^w,w]iv), 
F,[w]{v) -.^ F,[D^w,w]{v), 

F'[w]{v) := Fr[w]{v) + Fp[w]{v) + F,[w]iv), 
G^{w) G^{D^w, w) = eA'^w + F{w), 
G',[w]{v) ■.^eA^v + F'[w]{v). 

We conclude this section and chapter by citing a divergence-free row property 
of the cofactor matrix of the gradient of a vector- valued smooth function (a special 
case of Piola's identity) . This property will be used many times in the later chapters 
of the book. A proof of this property can be found in |321 page 440] . 

Lemma 1.4. Given a vector-valued function v = {vi,V2, ■ ■ ■ ,fn) '■ ^ ^ R"- As- 
sume V G [G^(r2)]". Then the cofactor matrix cof(Z?v) of the gradient matrix Dw 



^In an effort to clarify notation, we mostly use Greek letters to represent matrix-valued 
functions, and Roman letters to represent scalar functions throughout the book 



dF 



Fr-.fJ,^ ^ -----{r,p,z,x)fi,j{x), 



^ 9F dv 

Fp-Vv = y --—{r,p,z,x} — {x), 
^ dpi dxi 

dF 

Fz-v = —{r,p,z,x)v{x), 



of V satisfies the following row divergence-free property: 

div(cof(L'v))i = —{coi{Dv))ij =0 for i = 1, 2, • • • , n, 

where (cof(Z)v)), and {co{{Dv))ij denote respectively the ith row and the {i,j)-entry 
o/cof(£)v). 



CHAPTER 2 



Formulation of the vanishing moment method 



In this chapter we shall present the formulation of the vanishing moment 



method for fully nonlinear second order PDE ( 1.1 ). We also explain how the method 
was conceived and give some informal insights about the method. We note that 
the material of this chapter has a large overlap with that of [37 . 

For the reasons and difficulties explained in Chapter [T] as far as we can see, it is 
unlikely (at least very difficult if at all possible) that one can directly approximate 
viscosity solutions of general fully nonlinear second order PDEs using available 
numerical methodologies such as finite difference methods, finite element methods, 
spectral and discontinuous Galerkin methods, meshless methods, etc. In particular, 
the robust and popular Galerkin-type methods (such as finite element methods, 
spectral, and discontinuous Galerkin methods) for solving linear and quasilinear 
PDEs become powerless when facing fully nonlinear second order PDEs. From a 
computational point of view, the notion of viscosity solutions is, in some sense, an 
"inconvenient" notion for fully nonlinear second order PDEs because it is neither 
constructive nor variational. In searching for a "better" notion of weak solutions 
for fully nonlinear second order PDEs, we are inspired by the following simple but 
crucial observation: the crux of the vanishing viscosity method for the Hamilton- 
Jacohi equation and the original notion of viscosity solutions is to approximate a 
lower order fully nonlinear PDE by a family of quasilinear higher order PDEs. 

It is exactly this observation which motivates us to apply the above quoted idea 
to fully nonlinear second order PDE in |37| . To this end, we take one step 



further and approximate fully nonlinear second order PDE (1.1 1 by the following 
fourth order quasilinear PDEs [SfQ 

(2.1) eA^u'' + F(£|2u^, Vw^w^a;) in fl, e > 0. 

Here and for the continuation of the paper, we only consider the Dirichlet problem 
for ([TT]), so we suppose that 

(2.2) u = g on dfl. 
It is then obvious that we need to impose 

(2.3) ^ g on dn. 



However, the Dirichlet boundary condition (2.3) is not sufficient for well-posedness, 
and therefore an additional boundary condition must be used. Several boundary 
conditions could be used for this purpose, but physically, any additional boundary 
condition will introduce a so-called "boundary layer" . A better choice would be 



^Other higher order linear operators may be used in the place of A^n^, we refer the reader to 
|37| for more discussions on the choices of the regularizing operators. Here, we implicitly assume 
that —F is elliptic in the sense of 1421 Chapter 17], otherwise, | |2.1| needs to be replaced by 
-eA'^u^ + F{D'^u^,Wu^,u^,x) = 0. 
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one which minimizes the boundary layer. Based on some heuristic arguments and 
evidence of numerical experiments, we propose to use one of the following additional 
boundary conditions: 



(2.4) 



on dCt, 



(2.5) 



dAu" 
dv 



on dVl, 



(2.6) 



on 9f2, 



where v denotes the outward unit normal to d^l. 

We note that another valid boundary condition is the following Neumann 
boundary condition: 

— — = £ on oil. 

Of 

However, since this is an essential boundary condition, it produces a larger bound- 
ary layer than the other three boundary conditions, and therefore, we do not rec- 
ommend the use of this boundary condition. 



The rationale for picking boundary condition (2.4) is that we implicitly impose 
an extra boundary condition 



^Au' 
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on dil, 



which is a higher order perturbation of the original Dirichlet boundary condition 
(2.2). Intuitively, we expect that the extra boundary condition converges to the 



original Dirichlet boundary condition as e tends to zero for sufficiently large positive 
integer m. 



Remark 2.1. (a) We note that boundary conditions (2.4) and (2.5), which are 



natural boundary conditions for equation (2.1), have an advantage in PDE conver- 



gence analysis. Also, both boundary conditions are better suited for conforming 
and nonconforming finite element methods \39\ I62j , where as boundary condition 
in ( |2.6[ ) fits naturally with the mixed finite element formulation [38j . 

(b) From the PDE analysis viewpoint, the reason why high order boundary 
conditions such as (2.4)-(2.6) work better may be explained as follows. Since vis- 



cosity solutions generally do not have second or higher order (weak) derivatives, 
we do not expect to converge to u in H^{n) for s > 2 in general. Therefore, it 
is possible that errors in higher order derivatives, which could be big, would have 
small effects on the convergence of in the lower order norms if u'^ is constructed 
appropriately. Also, as we shall see later, the reason we do not impose homoge- 



neous boundary conditions in (2.4)-(2.6) is that the regularized solution inherits 



favorable properties such as strict convexity. 

To summarize, the vanishing moment method consists of approximating the 
(given) nonlinear second order problem 



(2.7) 
(2.8) 



in fi, 
on dil, 



11 



by the following quasilinear fourth order boundary value problem: 



(2.9) 
(2.10) 



,x) 



in ri, 
on dfl. 



(2.11) 



Au" = e, 



dv 



= e, or D u^v ■ v = e 



on 



on. 



Since equation (2.9) is quasilinear, we can then define the notion of a weak 



solution using the usual integration by parts approach. 



Definition 2.2. We define u'^ e H^{fl) with w|^^ = g to be a solution of (2.9| 
( pUj li if for aU V e H^{n) n H^ifl) 

2 ^\ 



(2.12) 



e(A^t^ Av) + {F{D^u', Vu^ ^^^ x), w) 



We now are ready to define the notion of moment solutions for (2.7)-(2.8) 



a weak {resp. strong) moment solution to p 
holds in i?^-weak {resp. H^-weak) topology. 



2.9 


-( 


2.11 


)i. 


em ( 


2.7 


|-( 


2.8 



imj_j.o+ '^'^ is called 
if the convergence 



Remark 2.4. (a) The terminologies "moment solutions" and "vanishing moment 
method" were chosen due to the following consideration. In two-dimensional me- 
chanical applications, often stands for the vertical displacement of a plate, and 
D^u^ is the moment tensor. In the weak formulation, the biharmonic term becomes 
e{D'^u'^ , D^v) which should vanish as e ^ 0+. This is the reason we call lim£_^o+ '"'^ 
(if it exists) a moment solution and call the limiting process the vanishing moment 
method. 

(b) Since weak moment solutions do not have second order weak derivatives in 
general, they are difficult (if at all possible) to identify. On the other hand, since 
strong moment solutions do have second order weak derivatives, they are naturally 



expected to satisfy equation (2.7) almost everywhere and to fulfill the boundary 



condition (2.10). In the remainder of this book, moment solutions will always 



mean weak moment solutions. 



As problem ( 2.9 1-( 2.11) is a quasilinear fourth order problem, one can com- 
pute its solutions using literally any well-known numerical methods, in particular, 
Galerkin-type methods such as finite element methods, spectral and discontinu- 
ous Galerkin methods. We note that ( |2.12[ ) provides a variational formulation for 
(2.9)-(2.11)i. Indeed, developing finite element numerical methods is one of two 



main goals of this book. In Chapter |4] and [5] we shall present comprehensive finite 



element and mixed finite element analysis for problem (2.7)-(2.8 

However, a natural and larger question is whether the vanishing moment method- 
ology will work. There are two ways to address this question. First, one can do 
many numerical experiments to see if the methodology works in practice. We indeed 
have done so (and beyond) in a series of papers |37L l38 , 39, 40, 62' (also see |61| ) 
for the Monge- Ampere equation. All numerical experiments of these papers show 
that the vanishing moment methodology works effectively. Second, one can give a 
definitive answer to the question by laying down its theoretical foundation, namely, 
proving the convergence (and rates of convergence if it is possible) (cf. |36j ) of the 
vanishing moment method. Partially accomplishing this goal is in fact the second 
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main objective of this book. In the next chapter, we shah give a detailed con- 
vergence theory for the vanishing moment method apphed to the Monge-Ampere 
equation in the radial symmetric case. We refer the interested reader to |36| for 
the convergence analysis in more general cases. 

We conclude this chapter by mentioning another intriguing property of the 
vanishing moment method, which was reported in [37 and discovered numerically 
by accident. When constructing the vanishing moment approximation (2.9), we 
restrict the parameter e to be positive (and drive it to zero from the positive side). 
An interesting question is what happens if we allow e to be negative (and drive 
it to zero from the negative side). In order words, we want to know the limiting 
behaviour as e \ 0+ of the following problem: 

(2.13) -eA^u'' + F{D^7f,\7u^,u^,x) ^0 in fl, 

(2.14) = 5 on dfl, 

dAu^ 

(2.15) Au" = ~e, or — — = ~£, or D'^u''v-v = -e on dVL. 

ov 

The numerical experiments of |37| (also see |61j ) indicate that in the case of 
the two-dimensional Monge-Ampere equation (cf. Chapter 6), that is, 

F{D'^v,\Iv,v,x) = f -Aet{D^v), / > 0, 



converges to the concave solution of the Dirichlet problem (2.7)-(2.8|! In the 
next chapter, we shall also give a proof for this numerical discovery in the radial 
symmetric case. 



CHAPTER 3 



Convergence of the vanishing moment method 

The primary goal of this chapter is to present a detailed convergence analysis 
for the vanishing moment method applied to the Monge-Ampere equation in the 
n-dimensional radial symmetric case. Such a result then puts down the vanishing 
moment method on a solid footing and provides a (partial) theoretical foundation 
for the numerical work to be given in the remaining chapters. 



3.1. Preliminaries 

Unless stated otherwise, throughout this chapter V, = Bfi{0) C R" {n > 2) 
stands for the ball centered at the origin with radius R. We do not assume is the 
unit ball because many of our results will depend on the size of the radius R. 

Suppose that / = f{r),f ^ and g ^ g{r) in ( |l.ll[ )- pri2l ), that is, / and g 



are radial. Then the solution u of (1.11 |-( |l.l2[ ) is expected to be radial, namely. 



u{x) is a function of r :— \x\ = ^Y^^j^ixj. We set u{r) :— u{\x\) — u{x), and for 

the reader's convenience, we now compute Au, A^u and det(£'^u) in terms of u (cf. 
[MllIO]). Trivially, 



dr Xj d /I 



dxj r ' dxj \r 
By the chain rule we have 

du(x) ^ , s dr ^ , ^Xj 

-7. = Mr)^ = Ur{r)^, 

oxj oxj r 

d'^u{x) Xjd^^^,^^^dfXj\^ ^ ( \\ , '^r{i") 



-(r) +Ur{r)——{ — ] ^ - [ -Ur[r) XiXj H d.y. 

oxi \ r / r \ r / r r 



dxjdxi r dxi dxi \ r ' r\r I r r 

Here, the subscripts stand for the derivatives with respect to the subscript variables. 

On noting that D'^u{x) is a diagonal perturbation of a scaled rank-one matrix 
xx"^ , and since the eigenvalues of xx^ are (with multiplicity n — 1) and \x\'^ — 
r^ (with multiplicity 1 and corresponding eigenvector x), then the eigenvalues of 
D^u{x) are 



4-Ur{r)\ = Urr{r) (with multiplicity 1), 
\r / r 



r 

A2 := — (with multiplicity n — 1). 
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Thus, 



Au{x) = Ai + (rt — 1)A2 = Urr{r) + 



n — 1 



Ur{r) 



det(L>^u(a;)) = Ai(A2) 



?i-i 



tij.r(r) 



Abusing the notion to denote u{r) by u{r), then problem ( |l.ll[ )"( [l.l2 ) becomes 
seeking a function u = u{r) such that 

(3.1) ^{i^rYX^f in(0,i?), 

(3.2) u{R)^giR), 

(3.3) M^(0) = 0. 



We remark that boundary condition (3.3) is due to the symmetry of u = u{r). 

Lemma 3.1. Suppose that r"-^f e L^{{0,R)) and f > a.e. on {0,R). Then 
there exists exactly one real solution if n is odd and there are exactly two real 
solutions if n is even, to the boundary value problem (3.1|-(3.3). Moreover, the 
solutions are given by the formula 



(3.4) 



u{r) 



5(i?)±/, (nL/(s))^ ds 
g{R)-jf {nLf{s)f^ ds 



if n is even, 
if n is odd 



for r e (0, R). Where 
(3.5) 



Lfis) 



ji-i 



fit)dt. 



Since the proof is elementary (cf. |59L I70j ). we omit it. Clearly, when n is even, 
the first solution (with "+" sign) is concave and the second solution (with "— " sign) 
is convex because Ur and Urr simultaneously positive and negative respectively in 
the two cases. When n is odd, the real solution is convex. 

Remark 3.2. The above theorem shows that u is at a point tq e (0, R) as long 
as / is at tq and Lf{rQ) ^ 0. Also, u is smooth in (0, R) if / is smooth in (0, R). 
We refer the reader to |59l 170) for the precise conditions on / at r = to ensure 
the regularity of m at r = 0, extensions to the complex Monge-Ampere equation, 
and generalized Monge-Ampere equations in which / ~ f(S/u,u,x). 

Similarly, it is expected that = u^{r) is also radial, and the vanishing moment 
approximation (2.9)-(2.11)i then becomes (cf. Chapter |6]) 



(3.6) 
(3.7) 

(3.8) 
(3.9) 



u'iR)^g{R), 



<(0) = 0, |<,(0)|<(X3, |<^^(r)| =o(— ^) 



as r 



ul^{R) 



n-1 
R 



u%R) 
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Later in this chapter, we shah analyze problem (3.6)-(3.9) which includes prov- 



ing its existence and uniqueness as well as regularities. After this is done, we then 



(3.1H3.3) 



show that the solution of (|3.6|)-(|3.9|) converges to the unique convex solution of 
Integrating over (0, r) after multiplying (3.6 ) by r"^^, using boundary condition 

n — 1 



(3.8) and 



lim r" 

r->-0+ 



n-1 







we get 
(3.10) 



er 



"-i(-^(r"-V),)^ + ^Kr=% m(0,i?). 

r"^^u^(r). A direct calculation shows 



Introduce the new function w^{r) 
that satisfies 



(3.11) 



er 



1 



in {0,R). 



Converting the boundary conditions (|3.8|)-(|3.9| to we have 
(3.12) 
(3.13) 

In addition, since 



w%0) = <(0) = 0, 



< = r"-^<, + (n - l)r"^"<, 

<^ = r"-i<^^ + 2(71 - l)r"-2„e^ + (n - l)(7i - 2)r"-3ye^ 



we have 
(3.14) 



(0) = o(^:;^) for < J < min{2, n ~ 1}. 



So we have derived from (3.6) a reduced equation (3.11), which is only of second 



order, hence, it is easier to handle. After problem (3.11 )-(3.13) is fully understood. 



we then come back to analyze problem (3.6)-(3.9). 



3.2. Existence, uniqueness, and regularity of vanishing moment 

appr oximat ions 



We now prove that problem (3.11 1-(3.13) possesses a unique nonnegative clas- 



sical solution. First, we state and prove the following uniqueness result. 



(3.11)-(3.13). Let 



Theorem 3.3. Problem (3.11)-(3.13) has at most one nonnegative classical solu- 
tion. 

Proof. Suppose that and are two nonnegative classical solutions to 

a+f3=n-l 
a,P>Q 

(wf )" - (u;^)" 

otherwise. 



and 
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Subtracting the corresponding equations satisfied by and W2 yields 

(3-15) -er-^{^c^l)^ + ^J^,^'^' = in (0,i?), 

(3.16) 0^(0) = 0^(0) = 0, 

(3.17) </)^(i?) = 0. 

On noting that > in [0, i?], by the weak maximum principle [32, Theorem 
2, page 329] we conclude 

max|0^(r)| = max{|0^(O)|, W(R)\\ - max{0, W{R)\\. 

If 4>'^{R) — 0, then 0'^ = 0. If (t>^{R) 7^ 0, then 0^ takes its maximum or minimum 
value at r = i?. However, the strong maximum principle |691 Theorem 4, page 
7] implies that <j)l{R) ^ 0, which contradicts with boundary condition (j)^{R) = 0. 
Hence, (/)^ = or = wf. The proof is complete. □ 

Remark 3.4. A more direct way to prove 0^ = is given as follows. Multiplying 
(3.15) by 0^, integrating by parts, and using the boundary conditions (3.16)-(3.17) 
yield 

(3.18) e / |0^(r)pdr + e / 0^(r)0"(r)dr 

Jo Jo '"' 



R ^ 







nr' 



On noting that 



" — l,^.^,^,^, n— 1 



— w'{rW{r)\'dr = Q. 



r=R n-1 



4>l{rW{r)dr=—{<t>^{r)y _ + / _(^^(r))^dr 



r=o ./q 2r 



n - 1 . ..„.^2 f" n-l 



(0^(i?)) + / (^^(r))^dr. 



2R ' " Jo 2r2 

so each term on the left-hand side of ( |3.18[ ) is nonnegative, hence, they all must 
be zero. The first term then gives = 0. Then 0^ = const. Hence, = by 



Next, we prove that the existence of nonnegative solutions to problem (3.11 1- 



( |3.13 ). 

Theorem 3.5. Suppose r^^^f e L^{{0,R)) and f >0 a.e. in {0,R), then there is 
a nonnegative classical solution to problem (3.11 1 -(3.13). 

Proof. We divide the proof into three steps. 

Step 1: Let V'" G C2([0, R]) be nonnegative and satisfy V'°(0) = V'°(0) = and 
i)^{R) = £i?"~^. One such an example is ^'^(r) = ^"r". We then define a sequence 
of functions {'0''}fc>o recursively by solving for A: = 0, 1, 2, • • • 



= Lf{r) / s^-'f{s)ds in (0,i?), 

(3.20) ^/;'=+i(0) =1/^,^+1(0) -0, 

(3.21) i/;^+i(i?) =£7?"-!. 
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We first show by induction that for any such sequence satisfying (3.19H(3.21 ), 
there holds ■0'^ > in [0, R] for all fc > 0. Note that ■0° > by construction. 
Suppose that V'' > in [0, R]. Since / > and / ^ 0, then 



Hence, ^p'''^^ is a supersolution to the linear differential operator on the left-hand 
side of (3.191. By the weak maximum principle |321 Theorem 2, page 329] we have 

min 7/'''+^ (r) > min{0, i;''+\0) , (R)} = niin{0, tp''+\R)}. 

[0,R] 

If ijj''~^^{R) < 0, and since = ei?"^^ > 0, the strong maximum prin- 

ciple [69'j Theorem 4, page 7] implies that 1/;'^+^ = const in [0, i?], which leads to 
a contradiction as tp'^^^{0) — 0. Thus, we must have ip''^^{R) > 0, and therefore, 
^fc+i > Q [0,i?]. By the induction argument, we conclude that ip'' > in [Q,R] 
for all k>0. 

It then follows from the standard theory for linear elliptic equations (cf. |32|, 
I42j ) that (3.19|-(3.21 1 has a unique classical solution tp^^^. Hence the (fc + l)th 
iterate ip''^^ is well defined, and therefore, so is the sequence {ip''}k>o- 

Step 2: Next, we shall derive some uniform (in fc) estimates for the sequence 
{'4'''}k>o- To this end, we first prove that ip''~^^{R) can be bounded from above 
uniformly in fc. Multiplying (3.19) by and integrating by parts yield 

(3.22) - eij^+^ (r)i^'=+i (r) '""^-t-e ^ (r) p dr 

r=0 . n 



n-l 



0^i(r)^'=+i(r)dr+ / , (0^^))"' Wl' 



nr 



^''+'{r)Lf{r)dr. 



It follows from boundary conditions (|3.20|) and (|3.21|) that 
(3.23) - eij^+\r)ij;''+\r) 



r=R. 



r=0 



e^i?"-V+^(i?). 



Integrating by parts gives 



(3.24) 



-^;f+i(r)V/'+'(r)dr = ^(^/.'=+i(r))' ' ^ 
r 2r r=o 



2r2 



{^p''+\r)) dr 



1 

2R 



2r2 



{iP''+\r)f dr. 



By Schwarz, Poincarc, and Young's inequalities, we get 



(3.25) 



i;''+\r)Lf{r)dr < 



\2p';+\r)\Ur + 



2 Jo 



CI 

2e 



2 pR 



[Lj{r)fdr 



for some positive constant Ci = Ci{R). 
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Combining ( 3.22 )-( 3.25) we obtain 

(3.26) -2£2i?"-iv/+i(i?) + £ r' \^';:+\r)\^ dr + '-^^^-^{i,''+\R)y 

Jo ^ 



£ Jo 



Let 



ri — 1 

Then from (3.26) we have 



e 



z'^ - bz - c < 0, 



which in turn imphes that 

b - Vb'^ + 4c b + Vb^ + Ac 
zi<z<Z2, where zi = ^ , Z2 = ^ . 

Since — Z2 < zi, the above inequahty then infers that \z\ < z^. Thus, there exists a 
positive constant C2 ~ C2{R,Lf) such that 



(3.27) 



Substituting (3.27) into the first term on the left-hand side of (3.26) we also get 



(3.28) e r \i:';+\r)\^ dr + 



— {^'^{r))'"'\^'^+'ir)\'dr 



r 



nr 



<C,:=^{Lf{R)y + 2e'R^^-^C2. 



Now using the pointwise estimate for linear elliptic equations [42] Theorem 3.7] 
we have 



(3.29) 



max|V^'=+i(r)| < f ^ + 
[OM] \ e e 



Next, we show that V'r^^ is also uniformly bounded (in k) in [0, i?]. To this 
end, integrating (3.191 over (0,r) after multiplying it by 7-"("~^\ and integrating 
by parts twice in the first term yield 

(3.30) ^^^(r) = -^^^^^^^gfey^ f V'^+^C.) ds 



er' 



i(n-l) 



s"("-i)%(s)ds Vre(0,i?). 
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Using L'Hopital's rule it is easy to check that the hmit as r — > 0+ of each term on 
the right-hand side of (3.301 is zero, hence, each term is bounded in a neighborhood 
of r = 0. Moreover, on noting that tp'^ > 0, by Schwarz inequahty, we have 



(3.31) 



< 



„n{n- 



i(n-l) 



Now in view of (3.28 1-( 3.31 1 we conclude that there exists a positive constant 

C4 



C4 = Ci{R,Lf) such that 
(3.32) 



inax|V'r^(r)|< 

10, -K] £ 2 



By ( |3.19[ ) we get 
(3.33) ^';+\r) 



--Lf{r) Vre(0,i?). 



Again, using L'Hopital's rule and (3.141 it is easy to check that the limit as r — ?> 0+ 
of each term on the right-hand side of (3.33) exists, and therefore, each term is 
bounded in a neighborhood of r = 0. Hence, it follows from (3.291 and (3.32) that 
there exists a positive constant C5 = C5(i?, Lf) such that 



(3.34) 



max \ ib: 

[0,R] 



rr 



r) < 



-n+l ■ 



To summarize, we have proved that ||?A'^^"'^||cj([o,fl]) !i C'(e, i?, n, Lj) for j = 
0, 1, 2 and the bounds are independent of k. Clearly, by a simple induction argument 
we conclude that these estimates hold for all fc > 0. 

Step 3: Since ||V''^llc2([o,i?]) is uniformly bounded in k, then both {4>''}k>o and 
{'4'r}k>o are uniformly equicontinuous. It follows from Arzela-Ascoli compactness 



theorem (cf. [321 page 635]) that there is a subsequence of {4'''}k>o (still denoted 
by the same notation) and ip G C^([0, R]) such that 



uniformly in every compact set E C (0, R) as k ^ 00, 
uniformly in every compact set E C (0, i?) as fc — >■ 00. 



Testing equation (3.19) with an arbitrary function x G Cq{{0,R)) yields 



71-1 



^^{4'Hr)r V'+'(r)x(Orfr-- / Lfir)xir)dr 



nr 

Setting fc — > CXI and using the Lebesgue Dominated Convergence Theorem, we get 



(3.35) 



'4'rir)Xr{r) dr + e 



1 



A{r)x{r) dr 



+ nr^(n~i) (^(^)) X{r)dr= j Lf{r)x{r)dr. 
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Since V' S C'^{[0,R]), we are able to integrate by parts in the first term on the 
left-hand side of (3.35), yielding 
rR , 



e{n — 1) 



^'r{r) + --^Mr)T -Lf{r) x{r) dr = 



r nr'^^ 
for all X e Cq ((0, i?)). This then implies that 

1 



, , , e(n — 1) , , , 



^(V(r))"-L/(r) = Vre(0,i?), 



i(n-l) 



that is, 



-er 



Thus, ip satisfies (3.111 pointwise in (0, i?). 



Finally, it is clear that -0 > in [0,i?], and it follows easily from (3.201 and 

V'(O) = Vr(0) = 



(3.211 that 



and Vr(i?) = ei?""^ 
So we have demonstrated that ip S C^([0,i?]) is a nonnegative classical solution to 
problem ( 3.11| )-(3.13). The proof is complete. □ 

Remark 3.6. (a) The proof at the beginning of Step 2 gives an estimate for the 
Neumann to Dirichlet map: 'ip'^~^^{R) ijj^^'^iR). 

(b) We note that the a priori estimates derived in the proof are not sharp in 
£. Better estimates will be obtained (and needed) in the next section after the 
positivity of Am^ is established. 

The above proof together with the uniqueness theorem. Theorem |3.3[ and the 
strong maximum principle immediately give the following corollary. 

Corollary 3.7. Suppose r"-'^f G L^{{0,R)) and f > a.e. in {0 , R), th en there 
exists a unique nonnegative classical solution to problem (3.11)-( [3.13 ). More- 
over, >0 m {0,R), e C3((0,i?)) if f e C°((0,i?)), and is C°° provided 
that f isC°°. 



Recall that — r" where and w'^ are solutions of (3.6)-(3.9) and 



(3.11 )-(3.13). Let w'^ be the unique solution to (3.11 )-(3.13l, as stated in Corollary 
|3.7[ define 



(3.36) 



u'{r) :^g{R) 



-—^w%s)ds Vre(0,i?). 



We now show that is a unique monotone increasing classical solution of 
problem (Is^el-p^gl). 



Theorem 3.8. Suppose f e C°((0,i?)) and f > in (0,i?), then problem ([3l3|- 
(3.9) has a unique monotone increasing classical solution. Moreover, is smooth 
provided that f is smooth. 



Proof. By direct calculations one can easily show that defined by (3.36) 
satisfies (3.6)-(3.9|. Since > in (0,i?), then is a monotone increasing 
function. Hence, the existence is done. 

To show uniqueness, we notice that is a monotone increasing classical so- 
lution of problem (3.6)-(3.9) if and only if is a nonnegative classical solution 
of problem (3.11 )-(|3.13 ). Hence, the uniqueness of (3.6)-(3.9) follows from the 
uniqueness of (3.11|-(3.13). The proof is complete. □ 
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3.3. Convexity of vanishing moment approximations 



The goal of this section is to analyze the convexity of the solution whose 
existence is proved in Theorem |3.8[ We shall prove that is strictly convex either 
in (0,i?,) or in (0, i? — cqe) for some e-independent positive constant cq. From 
calculations of Section 3.1 we know that D^u'^ only has two distinct eigenvalues 
Ai = u^j. (with multiplicity 1) and A2 — -^u^ (with multiplicity n — 1), and we 
have proved that A2 > in (0,i?), so it is necessary to show Ai > in {0,R) or 
in (0, R — cqe). In addition, in this section we derive some sharp uniform (in e) a 
priori estimates for the vanishing moment approximations u'^ , which will play an 
important role not only for establishing the convexity property for u'^ but also for 
proving the convergence of in the next section. 

First, we have the following positivity result for Au^. 

Theorem 3.9. Let be the unique monotone increasing classical solution of prob- 



lem (^3.6|-(3.9) and define 



Then 



(i) w'^ > Q in (0, R), consequently, Au"^ > in (0, R), for all e > 0. 

(ii) For any tq G (0,i?), there exists an eg > such that > er" 
Au^ > e in (ro, R) for e G (0, eo)- 



and 



Proof. We split the proof into two steps. 

Step 1: Since is monotone increasing and differentiable, then > in 



[0, R\ . From the derivation of Sect 
nonnegative classical solution of (3 
the Laplacian A we have 



ionJO 

3.iiD-((: 



we know that 



-(|3.13). Let ip^ 



" is the unique 
w^. By the definition of 



(3.37) ip' = < = r"-i<^ + {n- l)r" 

So > in (0, R) infers Au^ > in {0,R). 



'An". 



To show ip'^ > 0, we differentiate (3.11 1 with respect to r to get 

-e(n-l)(n-2) (w^)"-i 



+ e(n-l)."-(-^<)^ + 



i(n-l) 



1 



r.(n-l)2+n 



(u.-^)" = r"-V(r) in(0,i?). 



From (3.11), we have 



er 



1 



j^j.(n— l)^+n 

Combining the above two equations yields 



r 



(3.38) 



-ewi 



e[n - 



\f + 



1 



Substituting 
(3.39) -eip'„ 



into the above equation we get 
l)(n-2) (w^)"-!- 



-\f + 



1 



-L 



^n(n— 1) 

(n - 



,.(n-l)2+. 



-(^^^)">0 in(0,i?), 
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since f,Lf,w^ > in (0,R). This means that ip'^ is a supersolution to a hnear 
uniformly elhptic differential operator. By the weak maximum principle we get (cf. 
[321 page 329]) 

min ip%r) > min{0, ip^iO), ip%R)} = min{0, 0, i?"-^£} = 0. 
[o,R] 

Here we have used the fact that (p^{R) — R"^^Au^{R) = i?"^^£. Hence, 1^9^ > in 
[0,R], so Am^ > in [0,R]. 

It follows from the strong maximum principle (cf. j69|, Theorem 4, page 7]) 
that ip"^ can not attain its nonpositive minimum value at any point in (0,R). 
Therefore, (p'^ > in (0, i?), which implies that Am^ > in {Q,R). So assertion (i) 
holds. 

Step 2: To show (ii), let := w^ — er'^~^ = r'^~^{Au'^ — e). Using the identities 



we rewrite (3.38) as 



(3.40) -eVrr + 



e(n-l)(n-2) {w^) 



r 



n-l 



Hence, ■0'^ satisfies a linear uniformly elliptic equation. 

Now, on noting that > by (i), for any € (0, R) (i.e., rg is away from 0), 
it is easy to see that there exists an £1 > such that the right-hand side of (3.40) 
is nonnegative in {rQ,R) for all s G (0,ei). Hence, "0^ is a supersolution in (ro,i?) 
to the uniformly elliptic operator on the right-hand side of (3.40). By the weak 
maximum principle we have (cf. |32[ page 329]) 

min '0^(r) > min{0, (ro) , (R)} = min{0, Au^(ro) — e, 0}. 

Again, here we have used the fact that Au^{R) — e. 

Since Au'^(ro) > 0, choose = minjei, iAu'^(ro)}, then ip'^{ro) = Au^(ro) — 
£ > ^Au^{ro) > for £ € (0,£o). Thus, minf^^.j?] '0^('') > for £ e (0,£o). 
Therefore, > £r"^^, consequently, Au^ > £ in [ro,i?] for e G (0,£o). 

Finally, an application of the strong maximum principle (cf. [691 Theorem 4, 
page 7]) yields that > £r"~^, hence Au^ > £, in (ro,i?) for e G (0,£o). The 
proof is complete. □ 

Remark 3.10. The proof also shows that £0 decreases (resp. increases) as ro 
decreases (resp. increases), and := Au^ takes its minimum value e in [ro,i?] at 
the right end of the interval r = R. 

With help of the positivity of Au^ , we can derive some better uniform estimates 
(in e) for and u^. 

Theorem 3.11. Suppose f € C°{{0, R)) and f > in {0, R) . Let u" be the unique 
monotone increasing classical solution to problem (3.6) -(3.9). Define ~ r"^^uf, 
and — Am^ = u^^ -I- ^^u^. Then there holds the following estimates (at least 
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for sufficiently small e > 0): 



(i) h1lco([o,/?]) + / Krdr<Co, 

(ii) Wu^WcHlQM])^ 

(iii) |K||co([o,fl]) < 



(ii) Wu^WcmoM]) + ll^«''llc°([o,_R]) < C*!, 
C2 



(iv) ||«1|co(N,fl]) < VO<ro<i?, 



C4 



(v) |k^||c»([ro.fl]) < („„i)2 VO<ro<i?, 



er, 



(vi) / |w^(r)|2dr+ / r2("-i)|t;-'(r)|2dr < — , 



(vii)£ / r"-^^"|w^(r)|^dr 



— «(r))"i;^(r)dr < — Va < n - 1, 



(viii) e / r''~'\vt{r)Y dr + {ui{r)y'-'\v' {r)Y dr < — for n > 3, 

Jo Jo e 

pR pR (-i 

(ix) e r^-"\v'^{r)\^ dr + r^-"u'^{r)\v^ {r)f dr < — for n ^ 2, a < 1, 

Jo Jo ^ 

where Cj = Cj{R, f,n) > for j — 0, 1, 2, • • • ,8 are e-independent positive con- 
stants. 



Proof. We divide the proof into five steps 
Step 1: Since is monotone increasing, 



(3.4f) 



maxu^(r) < u^(R) — q(R). 
[o,R] ^ ' ~ ^ ^ ^ ' 



We note that the above estimate also follows from Aw^ > and the maximum 
principle. 

On noting that satisfies equation (3.11), integrating (3.111 over (0,i?) and 
using integration by parts on the first term on the left-hand side yield 



-ewl{R)+e{n-l) I -w''{r)dr + 



1 



1 



w^{r) 



dr 



Lf{r) dr. 



Because w^{R) — ei?" ^ and > 0, the above equation and the relation 

= r^' 



,,6 _ iniply that 



(3.42) 



w^(r) 



R 



dr= / \u%r)\''dr<nR[Lf{R) + e^R''-^]. 



It then follows from ( |3.41[ ), ( |3.42p and ( |3.36[ ) that 

(3.43) g{R) ~ nR[Lf{R)+e^R"-^]" < u^{r) < g{R) Vr G [0,R]. 

Hence, is uniformly bounded (in e) in [0,i?], and (i) holds. 
Step 2: Let 

n — ^ 1 



24 



By (|3.6[) we have 



(3.44) 



-{r"-'v'r)r + -(«)"),, = r"-\f in {0,R). 



It was proved in the previous theorem that > e in i?) for sufficiently small 
e > and it takes its minimum value e at r — R. Hence we have v^{R) < 0[^ 
Integrating (3.441 over (0, i?) yields 



r=R 1 r=R 
r=0 n r=0 



hence, 

therefore, 
(3.45) 



«(i?))" = nLf{R) + £?ii?"-iw^(i?) < nLf{R), 
<,(i?) = |<(i?)| < (nLfiR))". 



Here we have used boundary condition ( 3.8 ) and the fact that (R) < and > 0. 
By the definition of w^{r) := r"^^u^(r) we have 



(3.46) 



w'iR) = \w'{R)\ < R"-^uUR)\ < R^-\nLf{R)y 



Using the identity 



we get 

Hence, 
(3.47) 



n — 1 

w^(r) = AM''(r) = <,,(r) H w^(r), 



7? — 1 77 — 1 

<,(i?) = A^^(i?) - — <(i?) = e - —KiR). 



77 — ' 1 , , 77 — 1 1 

\uUR)\ < e + ^^\ul{R)\ <£+^^{nLf{R))'\ 



R 



R 



Step 3: From Theorem 3.9 we have that w^(r) > in (0, i?), and hence, is 
monotone increasing. Consequently, 

(3.48) maxw^(r) = max|w=(r)| <w''{R) < R'"''^ (nL f(R))" . 

[OM] [0,B] I ' 'I \ '/ 



Evidently, (3.48) and the relation u'^(r) — r" ^u^{r) as well as limj._>.o+ u'^{r) 
imply that there exists rp > such that 



(3.49) 



1 i 
maxw!;(r) = max uj(r) < - + ( — ) (nLflR)) " . 
lo,R] [o,Ry - 2 \roJ ^ " 



Hence, (ii) holds. 



In addition, since satisfies the linear elliptic equation (3.38 ) , by the pointwise 



estimate for linear elliptic equations 42, Theorem 3.7] we have 

(r)| < + + („ _ l)||r-i%|U 

(n- 1)2 



(3.50) maxlw: 



{<r\[ 



^This is the only place in the proof where we may need to require e to be sufHciently small. 
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Since = r" ^Au^ =: r" ^v^, it follows from (3.50) that for any ro > there 
holds 

(3.51) max Iv^r)] = max |Au^(r)| 

[ro,K] lro,R] 



< e 



^— 1 ,S \7l 



Thus, (iii) and (iv) are true. 

Integrating (3.44) over (0,r) yields 



(3.52) 



er"-^< + -«)" = i/ in(0,i?). 



By ( |3.52 ) and (3.49) we conclude that for any tq > there holds 



(3.53) max w^(r) = max 



[ro.R] 



(A«^M)^| <-(!+(-) )^ 



So (v) holds. 



Step 4'- Testing (3.11) with and integrating by parts twice on the first term 

e(n- 1) 



on the left-hand side, we get 

rB. 



2 Jo 

e(n — 1) , \,9 , 



2R 



— \yj%r)r+Ur= / Lfir)w%r)dr. 



Combing the above equation and ( |3.48 1 we obtain 



(3.54) 



Consequently, 



u;^(r)pdr + 



2r2 



|w"(r)|"dr + 



1 



\w%r)\"+' dr 



(3.55) 



<i?[e2i?"-2+i;(i?)] (ni/(i?))^ 



\r"-^Au%r)\^dr+^-^^^^—^ ^ r^^''-^^\uf.{r)\^ dr 



r u, 



;(r)|"+idr < R[e^R''-^ +Lf{R)] {nLf{R)) 



Hence, (vi) holds. 

Step 5: For any real number a < n — I, testing (3.52) with r~°'v'^ and using 
w^(i?) = £ we get 



(3.56) 



, ejn-l-a) 
+ 2 



„n-2-ai .£/ \|2 



w''(r)|^dr 



+ / — «(r))'V(r)dr= / _L/(r)«-(r) dr. 



On noting that > 0, > 0, and 
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it follows from (3.56) that 
(3.57) 



e{n-l-a) „-2-qi er m2 J 



1 1 



< — i?" 



^„+3-a|jj||2 



ein?{n — 1 — a) (n + 3 — a) 



-{ul{r)rv%T)dT 

dr Va < n — 1. 



This gives (vii) 
Recall that 



F AF F If 

V :— Au — Uj.^ H M^, 



and therefore, we can rewrite (3.44 1 as follows 



n — 1 



-e(r"-iw^)^ + = r"- V + — — «)" in (0, R). 

Testing the above equation with r^w^ for /3 > 1 — ?i and using v^{R) — e, we get 



e/3(n + /3-2) 



r"+^-3|„^(r)|2dr+ / r'^«(r))"->"(r)p dr 



n — 1 



Hence, 
(3.58) 



£/3(n + /3-2) 



< 



-""^+''|^;^(r)|2dr + / r"~i+'^|/(r)|2 dr 



2e 



Jo ' 

To continue, we consider the cases n — 2 and n > 2 separately. First, for n > 2, 
since v^{R) < 0, it follows from ( 3.57| ) with a — 1 and (3.58) with = that 



(3.59) 



2 7o 



r"-i|<(r)|2dr+ / «(r))"-^|t-"(r)|" dr 



< 



1 



r"-i|/(r)|2dr + 



2 ^ e(n2 - 4) 



When n = 2, we note that a = 1 is not allowed in (3.57). Let a < 1 be fixed in 
(3.57), set /3 = 1 — a in ( |3.58 ) we get 



(3.60) 



r^'"|i;^(r)|^dr + / r^-"<(r)|v'^(r)|^ dr 
Jo 

rR r- d4|| fl\2 



< 



1 

2eJo 



r^-°'\f{r)\^dr + 2R' 



£3 + ^ 11-^1'^°° dr 

e(l — q;)(5 — a) 



Hence, (viii) and (ix) hold. The proof is complete. 



□ 
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We now state and prove the following convexity result for the vanishing moment 
approximation . 

Theorem 3.12. Suppose f G C'^((0,i?)) and there exists a positive constant /q 
such that / > /o on [O.R]. Let denote the unique monotone increasing classical 
solution to problem (3.6|-(3.9). 

(i) If n — 2, 3, then either it^ is strictly convex in (0, R) or there exists an 
e-independent positive constant Cq such that is strictly convex in (0, R~ 
cqe). 

(ii) // n > 3, then there exists a monotone decreasing sequence {sj}j">o C 
(0, i?) and two corresponding sequences {£j}j">o C (0,1), which is also 
monotone deceasing, and {''jlj^o C (0,i?) satisfying sj \, 0^ as j —> oo 
and u^^{sj) > and R — r*= 0{e) such that for each j > 0, is strictly 
convex in {sj,r*) for all e € (0,ej). 

Proof. We divide the proof into three steps. 

Step 1: Let := r'^~^u'^ and w"^ :— Au^ = + ^^u^ = ''^r be same as 
before, and define 77^ := r"~^u^^. On noting that 



(3.44) can be rewritten as 

(3.61) -ev'„+[^+^-^^\v' ^r'^-'f + s{3~n)r--\' in(0,i?). 

So 77^ satisfies a linear uniformly elliptic equation. 

Clearly, 77^ (0) — 0. We claim that there exists (at least one) ri e (0, i?] such 
that rj'^{ri) > 0. If not, then ry*^ < in {0,R], so is m^^. This implies that is 
monotone decreasing in (0, i?]. Since ^^(0) = 0, hence, < in (0, i?]. But this 
contradicts with the fact that > in (0, R]. Therefore, the claim must be true. 

Due to the factor (3 — n) in the second term on the right-hand side of (3.61 1, 
the situations for the cases n < 3 and tt. > 3 are different, and need to be handled 
slightly different. 

Step 2: The case n = 2, 3. Since w'^ > 0, hence, 

(3.62) -erf„^\-^ + ^-^^\rf>Q m{0,R). 

Therefore, rj'^ is a supersolution to a linear uniformly elliptic differential operator. 
By the weak maximum principle (cf. 32, page 329]) we have 

min ?7^(r) > minjO, 77^(0), ry'^(ri)} ~ niinjO, ri^{ri)} = 0. 

Let — max{ri g (0,i?]; rj'^{ri) > 0}. By the above argument and the 
definition of we have 77"^ > in [0, r.^,], 77"^ (r^) = if 7^ i?, and rj^ < in (r*, R\. 
If = R, then 77^ > in [0, i?]. An application of the strong maximum principle 
to conclude that Ty*^ > in {0,R). Hence, uf.^ > in (0, i?). Thus, is strictly 
convex in (0,i?). So the first part of the theorem's assertion is proved. 

On the other hand, if < R, we only know that u"^ is strictly convex in (0, r*). 
We now prove that R ^ r^, = 0{e), which then justifies the remaining part of the 
theorem's assertion. 
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By (3.61) and the above setup we have 



in {r^,,R). 



Integrating the above inequahty over {r^ , r) for r < R and noting that rjf. (r* ) < 
we get 

-ev'r>-{r''-K) in(n,i?). 
n 

Integrating again over (r* , R) and using the fact that rj^ {r.^,) = and the following 
algebraic inequality 

1 - rr+i 



n + 1 R — n + 1 



we arrive at 



-eR'^-^UrriR) = -e7]%R) > 



m^{R-r,) 
n{n + 1) 



It follows from (3.47) that 

en{n + X)\urr{R)\ 



R-r^< 



< 



Rfo 

en{n + 1) 



(n- l)(nL/(i?))" 



CqS. 



Thus, 

(3.63) R-n^O{e), 
and is strictly convex in (0, R — cqe). 

Step 3: The case n > 3; First, By the argument used in Step 1, it is easy 
to show that rj'^ can not be strictly negative in the whole of any neighborhood of 
r = 0. Thus, there exists a monotone decreasing sequence {sj}j>a C (0, i?) such 
that Sj \ 0+ as j — T' oo and ri^{sj) > 0. 

Second, we note that 



e(3 - n)r"-S" = £(3 - n) 



n-l 



= e{i-n){^ + [n - l)r"-^<. 



Using this identity in (3.61), we have 



{n-l)e , K) 



if = r"-*[rV + £(3-n)(n- IX] in (0,i?). 



By (ii) of Theorem [3.11 we know that is uniformly bounded in [0, R]. Then for 
each Sj there exists an £j > (without loss of the generality, choose £j < £j-i) 
such that for e g (0, Sj) 

P/ + e(3-n)(n- IX] > in (sj,i?). 

Hence, rj"^ is a supersolution to a linear uniformly elliptic operator on {sj,R) for 
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Third, for each fixed j > 1, let r* = niax{r £ {sj,R]; rj'^{r) > 0}. Trivially, by 
the construction, r* > Sj^i > Sj. By the weak maximum principle (cf. [321 page 
329]) we have 



min f]^{r) > minlO, ?7'^(sj), 77^(r*)} > 0. 



Finally, repeating the argument of Step 2:, we conclude that is either strictly 
convex in (sj, R) or in {sj, r*) with R — r* = 0(e) for e £ (0, Sj). The proof is now 
complete. □ 



3.4. Convergence of vanishing moment approximations 

The goal of this section is to show that the solution u'^ of problem (3.61 



(3.9) converges to the convex solution u of problem (3.1)-(3.3|. We present two 



different proofs for the convergence. The first proof is based on the variational 
formulations of both problems. The second proof, which can be extended to more 
general non-radially symmetric case [36], is done in the viscosity solution setting 
|26j . Both proofs mainly rely on two key ingredients. The first is the solution 
estimates obtained in Theorem |3.11[ the second is the uniqueness of solutions to 
problem (Islll-p^l). 



Theorem 3.13. Suppose f <E C°((0,i?)) and there exists a positive constant fo 
such that f > fa in [0, i?]. Let u denote the convex (classical) solution to problem 
(|3.1|)-(3.3) and u"^ be the monotone increasing classical solution to problem (3.6)- 



(3.9) 



Then 
(i) 



converges to uniformly in 

,0 



(ii) 
(iii) 



u = lim£_^g+ u'^ exists pointwise and 

every compact subset of {0,R) as e 0^ . Moreover, is strictly convex 
in every compact subset, hence, it is strictly convex in [0, R] 
converges to weakly * in L°°((0, R)) as e — > 0^. 





u 



u. 



Proof. It follows from (ii) of Theorem 3.11 that ||M'^||ci([o,i?]) is uniformly 



bounded in e, then {u^}e>Q is uniformly equicontinuous. By Arzela-Ascoli com- 
pactness theorem (cf. |32[ page 635]) we conclude that there exists a subsequence 
of {w^}e>o (still denoted by the same notation) and m° € C^([0, i?]) such that 

u"^ — > u° uniformly in every compact set E C (0, R) as e — > 0^, 
< — > M° weakly * in i°°((0, R)) as e 0+, 

and u^iR) = g{R) implies that u°(fl) = g{R). 



In addition, by Theorem 3.12 we have that for every compact subset E C (0, R) 
there exists eo > such that E C (0, R—c^e) and is strictly convex in (0, R — cqe) 
for £ < £o- It follows from a well-known property of convex functions (cf. |48| ) 
that vP must be strictly convex in E and S Cj^'^, ((0,i?)). 



Testing equation (3.44) with an arbitrary function x G C'o((0,i?)) yields 



(3.64) e r r''-^vl{r)xr{r)dr-- ^ {u',{r))''xr{r) dr = ^ r^-^ f{r)x{r) dr, 
Jo nJo ^ ' Jo 



where as before — Au^ — ui^ + ^^^ui. 
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By Schwartz inequahty and (vi) of Theorem 3.11 we have 



r" '^v^{r)xr{r) dr 



n— 1 e / 



71-1 



el r 
Jo 

<e( / r2("-i)|t;^(r)|2dr)' ( 
R 1 
^ 







r 
R 



Xr{r) 



dr 



\Ax{r)\^dr 







as e 



0^ 



Setting e — 0+ in (3.64 1 and using the Lebesgue Dominated Convergence 
Theorem yield 

r-i? rR 



(3.65) 



KW)"x.(r)dr= / r^-'f{r)x{r)dr Vx G ^^((0, i?)). 



1 

Jo Jo 
It also follows from a standard test function argument that 

u°(0) = 0. 

This means that € C^([0, i?])nCio'J^((0, R)) is a convex weak solution to problem 
(3.1)-(3.3). By the uniqueness of convex solutions of problem (3.1)-(3.3), there 



must hold u = u. 

Finally, since we have proved that every convergent subsequence of {u'^}e>o 
converges to the unique convex classical solution u of problem (3.1 1-( [33| , the 
whole sequence {u^}e>o must converge to u. The proof is complete. □ 



Next, we state and prove a different version of Theorem |3.13[ The difference is 
that we now only assume problem (3.1 )-(3.3) has a unique strictly convex viscosity 



solution and so the proof must be adapted to the viscosity solution framework. 

Theorem 3.14. Suppose f G C'^((0,i?)) and there exists a positive constant /o 
such that f > fo on [0,i?]. Let u denote the strictly convex viscosity solution to 



problem (3.1 1-(3.3) and it^ he the monotone increasing classical solution to problem 



(3.6 1 -(3.9 1. Then the statements of Theorem 3.13 still hold. 



Proof. Except the step of proving the variational formulation (3.65), all other 



parts of the proof of Theorem 3.13 arc still valid. So we only need to show that m° is 
a viscosity solution of problem (3.1 )-( |3.3[ ), which is verified below by the definition 
of viscosity solutions. 

Let (j) G C^([0, i?]) be strictly conve^j^ and suppose that u'^ — (p has a local 
maximum at a point tq € (0, R), that is, there exists a (small) number Sq > such 
that (ro ~ So,ro + Sq) CC (0, R) and 

u%r) - 0(r) < u"(ro) - 0(ro) Vr G (ro - So,ro + So). 

Since (which still denotes a subsequence) converges to uniformly in [ro — 
Sqt^q + Sq], then for sufficiently small £ > 0, there exists r^ e (0, i?) such that 
r^ — ^ ro as £ — > 0+ and u"^ — <j> has a local maximum at r^ (see [321 Chapter 10] for 
a proof of the claim). By elementary calculus, we have 



^In fact, <j> can be taken as a convex quadratic polynomial (cf. | 18L I44p . 
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Because both u'^ and are strictly convex, there exists a (small) constant po > 
such that for sufficiently small e > 

Urr{r) < (l)rr{r) Vr G (ro - p, tq + p) , p<po, 

which together with an application of Taylor's formula implies that 

<(r) = (l)r{r) + 0{\r - r^l) Vr G (ro - p,rQ+ p), p < po- 



Let X G C§{{ro ^ P,ro + p)) with x > and x(?'o) > 0. Testing ( |3.44[ ) with x 
yields 



(3.66) 



1 

27ip 



{mr)r)rX{r)dr 



1 



2P 



6,(r))"-Vrr(r)x(r)dr 



> 



2p 



[«(r))"-i + 0{\r - r,|"-i)]<,(r)x(r) dr 



/-i-o+P 



i-o-p 



> 



2np 

Y rro+P 

2p 



(«(r))"),. + 0(|r-r,r-i)<,(r) x{r) dr 



r''-'[nr)x{r)+evt.{r)xr[r)\ dr 



ra+p 



K{'^)Xr{r) dr 



ra-p 



for some positive p-independent constant Cg. Here we have used the fact that 
u'^rr 5^ 0, X > in [ro — p, ro + p] to get the last inequality. 
From (vi) of Theorem 3.11 we have 

(3.67) e / r"-i«^(r)xr(r)dr 



r""^w''(r) Xrr{r) H Xr(7-) 

r 

ro-p 



< e 



R 



|w^(r)pdr 



\^x{r)?dr 



/ |Ax(r)|2dr 



(3.68) 



Setting e — 7> 0+ in (3.66) and using (3.67) we get 

^ rro+P 



2p 



.{r)r-'4>rr{r)x{r)dr 



ro-p 



> 



I rro+P 



2p 



r''-^f(r)x{r)dr-C,p- 



ro-p 



ro+P 



Ur{r)Xr{r) dr. 



ro-p 



Where we have used the fact that u^. converges to u° weakly * in L°° {{0, R)) to 
pass to the limit in the last term on the right-hand side. 

Finally, letting p — >■ O"*" in (3.68) and using the Lebcsgue-Besicovitch Differen- 
tiation Theorem (cf. [32]) we have 

{Mro)r''^rr{ro)xiro) > r^'' f iro)x{ro) ■ 
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Hence, 



4>rr{rQ) > /(ro), 



so M° is a viscosity subsolution to equation (3.1 ) 
Similarly, we can show that if u'^ - 



for a strictly convex function </> e Cq{(0, R)), there holds 



'^0 



assumes a local minimum at ro G (0, i?) 
!)), there holdf 



so u° is also a viscosity supersolution to equation (3.1 1, 
solution. The proof is complete. 



Thus, it is a viscosity 
□ 



3.5. Rates of convergence 

In this section, we derive rates of convergence for u'^ in various norms. Here 
we consider two cases, namely, the n-dimensional radially symmetric case and the 
general n-dimensional case, under different assumptions. In both cases, the lin- 
earization of the Monge- Ampere operator is explicitly exploited, and it plays a key 
role in our proofs. 



Theorem 3.15. Let u denote the strictly convex classical solution to problem (3.1 ) 



(3.3) and be the monotone increasing classical solution to problem (3.6) -(3.9 1 



Then there holds the following estimates: 
(3.69) 



( / e%r)\ur{r) - u%{r)\^ dry < e'i 



r"-i|AM(r)- Au"(r)pdr < ei di 



(3.70) 



where Cj — Cj{\\r" ^Aur\\L2) for j = 10,11 are two positive e -independent con- 
stants, and 



(3.71) r(r) 

Proof. Let 

V :— Au — Ur- 



K)" ^ iKY 

u — 



n — 1 



J2iMr)yKir)r-'-^ >0 in(0,i?]. 

j=0 



n — 1 



On noting that (3.6) can be written into (3.44), multiplying (3.11 by 



and 



subtracting the resulting equation from (3.44 1 yield the following error equation: 
(3.72) 



+ - [K)" - = in (0, R). 

Testing ( 3.72[ ) with using boundary condition e^(0) — e'^{R) = we get 



(3.73) 



r"-i<(r)e^(r) dr + 



1 



^(r)|e^(r)pdr = 0. 



where 9"^ is defined by (3.71 ) 
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Integrating by parts on the first term of (3.73) and rearranging terms we get 

.fl 1 
(3.74) e r"-^\Ae'{r)\^dr+- 9" {r)\eUr)\^ dr 
Jo n Jq 

= ei?"-iAe^(i?)e^(i?)-£ / r''~\r{r)el{r) dr. 

Jo 

We now bound the two terms on the right-hand side as foUows. First, for the 
second term, a simple appHcation of the Schwarz and Young's inequahties gives 



(3.75) e / r"-ii;^(r)e^(r)dr < 



1 

An 



O^r) 



\vr{r)\'^ dr. 



Second, to bound the first term on the right-hand side of (3.74), we use the bound- 
ary condition v'^{R) = e to get 



|Ae^(i?)| = \v{R) - v^{R)\ = \v{R) - e\ < \v{R)\ + 1 =: M, 

|i?"-ie^(i?)|^ = / ((r"-ie^(r))2)^dr = 2 / r2("-i)e^(r) Ae'^(r) dr 
Jo ^ Jo 

r''-^\Ae%r)\^dry U r3("-i)|ef.(r)|2 drl 



and 



< 2 

Hence by Young's inequahty, we get 
(3.76) |£i?"-iAe'^(i?)e^(i?)| 



rR \ - / 

<V2eM( r"~i|Ae^(r)|2dr)'( / r^^""!) |e^(r)|2 dr 



< 



< 



r"-i|Ae"(r)pdr + 2eAft( / r3("-i)|e^,(r)p dr 



An 



r"--^|Ae"(r)|"dr-|- — / ^ (r)|e^(r)|Mr -f e^r^M^C 



for some e-independent constant C = C{f,R,n) > 0. 
Combining ([3J4])-([3J6]) yields 

cR 



(3.77) 



r"-^|Ae^(r)|^dr + 



1 



< 2e^n 



^r)\eUr)\'dr 

R j,2(n-l) 



^(r) 



\vr{r)Ydr + e'inM^C. 



Thus, (3.69) and (3.70) follow from the fact that \\r 



LOO < OO. 



□ 



Corollary 3.16. Inequality (3.69) implies that there exists an e-independent con- 
stant C > such that 

(3.78) ^j\''-^\ur{r) - ul{r)\^ dry < CCio- 

Since the proof is simple, we omit it. 
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Theorem 3.17. Under the assumptions of Theorem \3.15[ there also holds the 
following estimate: 



(3.79) 



r''-'\u{r)-u^r)\Ur) <eC 



^12 



for some positive e-independent constant C12 = Ci2{R, n, u, Cn). 



Proof. Let O'^ be defined by (3.71 1, and e*^, v and v'^ be same as in Theorem 



3.15 Consider the foUowing auxihary problem: 



(3.80) 
(3.81) 
(3.82) 



0,(0) = 0. 



in (0,i?), 



We note that the left-hand side of (3.80) is the linearization of (3.11 at 9'^ . 

Since 0^ > in (0,i?], then ( |3.80[ ) is a linear elliptic equation. Using the fact 
that ci > r"^^(0^)^^ > Co > in [0, i?] for some e-indcpcndcnt positive constants 
Co and ci, it is easy to check that problem (3.80)-(3.82| has a unique classical 
solution 6. Moreover, 



(3.83) 



r"-Vrr(r)pdr-+ /" r"-^\<j)r{r)\^dr<C [ r"-i|e'^(r)p dr 
Jo Jo 



for some e-independent constant C = C{f,R,n,ca,ci) > 0. 

Testing ( |3.80[ ) by e^ using the facts that (t>r{0) = (f>{R) = 0, e^R) = and 
v^{R) = £ as well as error equation (3.72) we get 



(3.84) 



r"-^|e"(r)|"dr = - 



'(r)(/)r(r)e^(r) dr 



e I r'^~^vf.{r)4)r{r) dr 
Jo 

eR"-'^v'{R)(j)r{R)-£ 

Jq 



r"--^v=(r)A0(r) dr 



R 



£^i?"-Vr-(i?) +e / r"-i[w(r) -'y"(r)]A<?!)(r)dr 
Jo 

-e r''-^v{r)A(j){r)dr, 
Jo 



where we have used the short-hand notation Aip — r" ^[(firr + {n — l)r '^(f>r]- 
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For each term on the right-hand side of (3.84 1 we have the following estimates: 



R 



{r-^(j}r{r))rdr 



r"-i|(/)^(r)|2dr 



[v{r) — ti^(r)]A0(r) dr 



<e(^ r"-Xr)-i;'^(r)|2dr 
<eJCii(y" r"-i|A0(r)pdr 

rR. 



r"-i|A<?!>(r)|2dr 



r'^-^v{r)A(j){r)dr -CeU r"-Xr)|2 dr) ' ( / r"-i|A(/)(r)|2 dr) 





Substituting the above estimates into (3.84 1 and using (3.831 we get 



(3.85) 



r"-^|e-'(r)|^dr 




r"-i|</)^(r)pdr 



+ £(e^Cn+C„)(^ r"-i|A0(r)| = 
< 4e(£i?t + eV^R-^ + e^Cn + Cu)c(^ 



dr 



\e'{r)\^dr 

for some e-independent constant C„ = C{u) > 0. 

Hence, by ( |3.85| ) we conclude that p.79| holds with C12 = 4(£i?? +ev^i?"i + 
e^Cii + Cu)C. The proof is complete. □ 

Remark 3.18. The argument used in the above proof is so-called duality argu- 
ment, which has been the main technique and used extensively in the finite clement 
error analysis to derive error bounds in norms lower than the energy norm of the 
underlying PDE problem (cf. [13L I22j and the references therein). However, as far 
as we know, the duality argument is rarely (maybe has never been) used to derive 
error estimates in a PDE setting as done in the above proof. 

Since the proofs of Theorem |3.15| and |3.17| only rely on the ellipticity of the 
linearization of the Monge-Ampere operator, hence, the results of both theorems 
can be easily extended to the general Monge-Ampere problem ( |l.ll[ )-( [l.l2 1 and its 
moment approximation (2.9)-( 2jT| )i p| 

Theorem 3.19. Let u denote the strictly convex viscosity solution to problem 
(1.11)-(1.12) and he a classical solution to problem (2.9)-(2Tl li. Assume 
u £ W n H^(Q) and is either convex or "almost conve^' in £7. Then 



This observation was pointed out to the first author by Professor Haijun Wu of Nanjing 
University, China, and the proof for | |3.86| and ( |3.87[ l is essentially due to him. 

^"Almost convex" means that u'^ is convex in Q minus an e-neighborhood of dQ, see Theorem 

3.12 for a precise description. 
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there holds the following estimates: 
(3.86) 

(3.87) 

(3.88) 



ib 



Vu- Vii^pdxj ' < £1 Ci3, 

Au- Au"]'^ dxj ' < £i Ci4, 
1 



where Cj = Cj{\\V Au\\l2) for j = 13, 14, 15 are positive e-independent constants. 

Proof. Since the proof follows the exact same lines as those for Theorem |3.15[ 
we just briefly highlight the main steps. 



First, the error equation (3.72) is replaced by 

(3.89) eAv" + det{D'^u)-det{D'^u^) =0 

where v'^ = Au"^ . 

Next, equation (3.71) becomes 

(3.90) 



in r2, 



i"Ve" • Ve^ dx 



Ae^—dS-e / Vw Vendee, 
on ov 7n 

where 

(3.91) 9" ■.= co{{tD'^u + {l-t)D'^u'') for some t e [0, 1], 

now stands for the cofactor matrix of tD^u + (1 — t)D^u^. Since u is assumed to 
be strictly convex and is "almost convex" , then there exists a positive constant 
^0 such that (see Chapter |4| 



It remains to derive a boundary estimate that is analogous to (3.76). To the 
end, by the boundary condition w^|ao = e and the trace inequality we have 



(3.92) 



dn 



Ae''—dS<e{e\d^\ + \\Au\\l, 



(dn)) 



de" 



dv 



L2(an) 



<£M||Ve^|l£.(^)| 



< 2l|Ae^||i.(n)+M^e||Ve^||l.(^) 



,e||2 



L2(f2) 



f l|Ve11i.(0) + 



The desired estimates (3.86) and (3.871 follow from combining (3.90) and (3.92). 

Finally, (3.88) can be derived by using the same duality argument as that used 
in the proof of Theorem |3.17| We leave the details to the interested reader. □ 

Remark 3.20. The convergence rates proved in Theorem |3. 15^3.191 have been 
observed in numerical experiments. We refer the reader to Chapter [6] for details. 
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3.6. Epilogue 



We like to comment that the analysis of Section [3T 3.5 



can be easily extended 
to the cases of the other two boundary conditions in (|2.11). We note that in the 



case (2.11)2 boundary condition (3.9 1 should be replaced by 

- 1 



R 



and (3.131 should be replaced by 

wl^{R) - 



R 



We also reiterate an interesting property of the vanishing moment method 
which was briefly touched on at the end of Chapter [2] That is, the ability of 
the vanishing moment method to approximate the concave solution of the Monge- 
Ampere problem ( |l.ll| -( |l.l2 1. This can be achieved simply by letting e 0~ 
in (2.9)-(2.11)i. This property can be easily proved as follows in the radially 



symmetric case. 

Before giving the proof, we note that for a given / > in fJ, equation ( 1.11 1 does 
not a have concave solution in odd dimensions (i.e., n is odd) because det{D"'u) — f 
does not hold for any concave function u as all n eigenvalues of Hessian D^u of a 
concave function u must be nonpositive. On the other hand, in even dimensions 
(i.e., n is even), it is trivial to check that if u is a convex solution of problem ( |1.11[ )- 
(1.12) with (7 = 0, then — u, which is a concave function, must also be a solution of 
problem ( [ril] )-( [rT2l ). 

Next, by the same token, it is easy to prove that if is a convex or "almost 
convex" solution to problem (2.9|-(2.11 li, then — w^, which is concave or "almost 
concav^', must also be a solution of (2.9)-(2.11 )i. 

Finally, let n be a positive even integer, it is easy to see that changing to 
-u" in ([2^ 

let S := After replacing e by 



(2.11 )i is equivalent to changing e to - 
S and u'^ by :— 



e in (2.9)-(2.11)i. For e < 0, 



satisfies the same set of equations (3.6)-(3.9) with S{> 0) in place of s. Hence, 



m l3^l3:9l), we see that 



by the analysis of Section 3.2 3.5|we know that there exists a monotone increasing 



solution u° to problem (3.6)-(3.9) with e being replaced by S, which satisfies all 



the properties proved in Section 



3.2 3.5 Translating all these to 



conclude that problem ( |3.6[ )-( |3.9 1 for e < has a monotone decreasing solution 
which is either concave or "almost concave" in (0, R) and converges to the unique 
concave solution of problem ( |l.lip-(|1.12 ) as e 0^. In addition, satisfies the 
error estimates stated in Theorem 13. 151 and 13.171 

The final comment we like to make is about the possible but well-behaved 
boundary layer generated by the vanishing moment solution u^. In the worst case 
scenario, the boundary layer, where it^ may cease to be convex, is confined in an 
0(£)-neighborhood of the boundary dil. This nice behavior of the boundary layer 
can be exploited in numerical computations. Indeed, in Chapter [7] we propose 
an iterative surgical procedure to take advantage of this property of the (possible) 
boundary layer. We refer the reader to Chapter [7] for the detailed description of the 
procedure and numerical experiments which show the effectiveness of the proposed 
iterative surgical procedure. 



function (fi'^ is said to be "almost concave" in Q if it is concave in Q minus an 0(e)- 
neighborhood of the boundary 9f2 of f7. 



CHAPTER 4 



Conforming finite element approximations 



The goal of this chapter is to construct and analyze finite element approx- 



imations for the general fully nonlinear second order Dirichlet problem (2.7)-(2.8) 
based upon the vanishing moment methodology introduced in Chapter 2 and fur- 
ther analyzed in Chapter [s] Letting u"^ be the solution to problem (2.9)-(2.11 d, we 
construct and analyze conforming finite element methods to approximate using 
a class of finite elements such as Argyris, Bell, Bogner-Fox-Schmit, and Hsieh- 
Clough- Tocher elements (cf. |22) ). As a result, we obtain convergent numerical 
methods for fully nonlinear second order PDEs. 

We note that finite element approximations of fourth order PDEs, in particular, 
the biharmonic equation, were carried out extensively in the seventies for the two- 
dimensional case |22) . and have attracted renewed interests lately for generalizing 
the well-known two-dimensional finite elements to the three-dimensional case (cf. 
[73[ \77\ 178) ). Although all of these methods can be readily adapted to discretize 
problem (2.9)-(2.11 li, the convergence analysis does not come easy due to the 
strong nonlinearity of the PDE ( 2.9 1 . For example, to use the standard perturbation 



technique for deriving error estimates (a technique successfully used for linear and 
mildly nonlinear problems), we would have to assume very stringent conditions 
on the nonlinear differential operator F, which would rule out many interesting 
application problems, and hence, should be avoided. Instead, we assume very mild 
conditions on the operator (see Section 4.1 for details), and use a combined fixed- 



point and linearization technique to simultaneously prove existence and uniqueness 
for the numerical solution, and also derive error estimates. 

The remainder of the chapter is organized as follows. First in Section |4.H we 
give additional notation, and then define the finite element method based upon 
the variational formulation (2.12). Next, we make certain structure assumptions 



about the nonlinear differential operator F which will play an important role in our 
analysis. In Section 4.2 we show existence of solutions of the linearized problem 
and prove stability and convergence results of its finite element approximations. 
The main results of the chapter are found in Section |4.3[ where we use a fixed 
point argument to simultaneously show existence, uniqueness, and convergence of 



the finite element approximation of ( 2.9 1-( 2.11) 



4.1. Formulation of conforming finite element methods 

First, we introduce the following function space notation: 



Va 



{v e V; v\ 



g}- 



Let 7^ be a quasiuniform triangular or rectangular partition of f2, and let 
C be a conforming finite element space consisting of piecewise polynomials 
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of degree fc > 4 such that for any v E V D H^{n), we have 



(4.1) inf \\v-v^,\\h, <Ch'-^\v\\Hi 



j = 0,1,2, ^ = niin{s, fc + 1}. 



Let 
(4.2) 

Based on ( |2.12[ ), we define the finite element formulation of ( |2.9[ )-( [2jT ) as to 
find u% G such that 



on 



9} 



(4.3) e(Au|, Avh) + {FiD^ul, Vii|, u%, x), v^) 



2 dvh \ 



Let be the solution to (2.121 and let u| be a solution to (4.3). The primary 
goal of this chapter is to derive error estimates of ~ ""1, which then means we 
need to first prove that there exists m| € Vg solving (4.3), and that u| is unique. 
Clearly, we must assume some structure conditions on the nonlinear differential 
operator F to achieve any of these goals. Indeed, the assumptions that we make 
will play an important role in our results and in the techniques to derive them. We 
refer to Section |1.2| for the notation used in this chapter. 



Assumption (A) 

[Al] There exists eg G (0, 1) such that for all e € (0,eo]j there exists a locally 
unique solution to ([tjl)-( |2.11[ )i with G H^i^l) (s > 3). 

[A2] For £ e (0,£o], the operator (G^[u'^])* (the adjoint of G'^[u^]) is an iso- 
morphism from Vq to Vq . That is for all e V^* (the dual space of V^*), 
there exists w S Vb such that 



(4.4) 



Here, (•, •) denotes the dual pairing between Vq and . Furthermore, 
there exists positive constants Cq — Co(£), Ci = Ci{e) such that the 
following Carding inequality holds: 



(4.5) 



{G'Au']{v),v)>C^\\v\\]j.-C^\\v\\l. 
and there exists C2 — C2(£) > such that 

\\F'[u-]\\yy,<C2. 



where 

ll^'Hllyv. 



sup 

veVo 



\F>']{v)\\h- 



sup sup 

\v\\h^ vevowevo ll^'ll//^ |1 w||//2 

Moreover, there exists p > 2 and Cr = Cr^e) > such that ii tp E L^(ri) 
and V eVo satisfies ( |4.4[ ), then v S HP{il.) and 

\\v\\hp < Cb.\\ip\\l^- 

[A3] There exists a Banach space Y with norm || • ||y that is well-defined and 
finite on V'^, and a constant C > 0, independent of e, such that 



{F'[u%v),w) 



sup ■ 



\F'[y] 



< C. 



\y\\Y 



41 



[A4] There exists a constant C > independent of e such that 

WI'^u'Wy <C\\u'\\y, 

where X^u"^ G Vg denotes the finite element interpolant of u^. 
[A5] There exists a constant S = 6{e) G (0, 1), such that for any Wh & Vg with 

\\I^u'^ — Wh\\ff2 < S, there holds 

- F'[wh]\\yy, < Lih)\\u' - whWh^, 
where L{h) ~ h) may depend on both h and e and satisfies L(h) — 

Remark 4.1. (a) Conditions [A1]-[A5] are fairly mild, and a very large class of 
fully nonlinear second order differential operators satisfy these requirements (cf. 
Chapter [6]) . Clearly, we must assume [Al] in order for the finite element method 



(4.3) to have any significance, and the regularity requirements of are needed to 
obtain any meaningful error estimates. 

(b) Condition [A2] is naturally satisfied if ~F is elliptic at (cf. |42[ Chapter 
17]), and the regularity requirements are expected to hold provided that and dfl 
are sufficiently regular. 

(c) By standard interpolation theory |22l I13j , there holds 

(4.6) \\u^ -I''u''\\hj <Ch^-mu''\\Hi j^0,l,2, ^ = min{s, fc + 1}. 

(d) Condition [A5], which states that F' is locally Lipschitz near u^, is the 
strongest requirement among the five listed, and it is the authors' experience that 
this is the most difficult property to verify. As one may expect, this assumption 
plays an important role in the fixed point argument, which is needed in our analysis 
to obtain existence, uniqueness, and error estimates of the finite element method 



(4.3) 



4.2. Linearization and its finite element approximations 



To construct the necessary tools to analyze the finite element method (4.3), 



we first study finite element approximation of the linearization of (2.9). We note 
that the materials of this section have an independent interest within themselves. 
To the best of our knowledge, finite element error estimates for non-coercive linear 
fourth order problems have not been explicitly given in the literature before. 

4.2.1. Linearization. For given ip e Vq and ip E H^i{drt), we consider the 
following linear problem: 

(4.7) G'^[u'']{v) = ^ inn, 

(4.8) v^O on dn, 

(4.9) Aw = ^ on dn. 



Multiplying the equation (4.7 1 by w e Vq, integrating over f2, and integrating 
by parts, we obtain 



{G',[u']{v), w) = e{Av, Aw) + {F' [u']{v) , w) - e l^Av, 
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Based on this calculation, we define the weak formulation of (4.7)-(4.9) as to find 
w e Vq such that 



(4.10) 
where 



a%v, w) ^ {(p,w) + e{ 



dn 



a^(w,w) := e{Av,Aw) + {F'[u^]{v),w). 

In view of assumptions [A1]-[A2], we immediately have the following theorem. 

Theorem 4.2. Suppose assumptions [Al]-[A2] hold. Then there exists a unique 
solution V G Vo to (4.10). Furthermore, there exists C3 = Csi^) > such that 



(4.11) 



Proof. From the Garding inequality (4.5) and the fact {G'^[u^]) is injective 
on Vq, it follows that G'^[u^] is an isomorphism from Vq to Vq using a Fredholm 
alternative argument [U Theorem 8.5]. 

We now claim that there exists C5 = Csis) such that \\v\\l2 < Cs{\\'p\\h-^ + 
'^"^"i/-5(an))' would exist sequences {(Pm}m=i ^ V^o*' {^m}m=i ^ 

H^^{d^), and {vm\m=i ^ Vq such that 

/ dxv \ 

{G'^[u%V„^),w) = {ip,n,w) + e (tpm, ) 



on 



but 



Without loss of generality, we may as well suppose ||fm||L2 = 1 (and therefore 
ll<^m||//-2 +£|l^'nll^-|(gf^) ^- as m ^ oo). In light of ([45]), {v^}^^^ is bounded 
in Vq, and hence by a compactness argument, there exists a subsequence {vmj}m=i 
and V £ Vq such that 



(4.12) 
(4.13) 

Therefore, 



V weakly in Vq, 



{G',[u^]{v),w) ^0 yweVQ. 



Since G^[u'^] is an isomorphism, v = 0. However (4.13) implies that ||w||i,2 = 1, a 
contradiction. 

Hence there exists such that 
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and therefore by (4.5 1 and a trace inequality, we have 

C,\\v\\l, < e{Av,Av) + {F'[u']{v),v) + Co\\v\\ 
^a'{v,v)+Co\\v\\l2 



L2 



/ dv\ 

= + 'A, +Co\\v\\l2 

< c[MH-2+eM^_.^^^^+Co\\vh.)\\v\\H' 

< Cil + CoCs){Mh-^ + em^^i^^^^^)\\v\ 



Dividing by Ci||t;||jj2, we obtain ( [iHj ) with C3 = CCf ^(1 + CoCs). □ 

4.2.2. Finite element approximation. Let Vq C Vq be one of the finite 
dimensional subspaces of degree fc > 4 defined in Section |4.1[ Based on the varia- 



tional formulation (4.101, we define the finite element method for (4.7)-(4.9) as to 
find Vh € such that 

dwh " 



(4.14) 



dv 



Using a modification of the well-known Schatz's argument (cf. [13, Theorem 
5.7.6]), we obtain the following result. 

Theorem 4.3. Let assumptions [A1]-[A2] hold and suppose thatv G H''(il) {s > 3) 
is the unique solution to (4.101. Then for h < ho{e), there exists a unique solution 
vji e to (|4.14|), where 



(4.15) = \i ^ ^ 



*/Co^0, 
*/Co = 0, 



min{p, fc -|- 1}. 



Furthermore, there holds the following inequalities: 

(4.16) \\v-Vh\\m<Cih'-^v\\H^, 

(4.17) \\v-Vy,\\L-<C,h'+--^\\v\\H^, 

where 

C4 = C4(e) = CCf ^C2, C5 = CsCe) = CC^^CICr, I = min{s, fc + 1}. 

Proof. To show existence, we begin by deriving estimates for a solution to 
(4.14) that may exist. We start with the error equation: 

a%v-Vh,Wh)=Q ywh<GV^. 
Then using (4.5) and [A2], we have for any Wh € Vq 
Ci\\v - Vh\\]j2 

= a^{v - Vh,v ~ Vh) + Co\\v - 



= a^{v -Vh,v - Wh) + Co\\v - Vh\\] 

< e\\A{v - Vh)\\LA\^{v - Wh)\\L-- 



\F'\u' 



\vv 



.\\v - fhllif^ll?; - WhWrn + Co\\v - Vh\\ 



L2 



< CC2\\v ~ Vh\\HA\'" - WhWH-^ + Ca\\v - Vh\\ 



L2- 
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Thus, by (|4l]) 

(4.18) Cillf - vh\\j,2 < CC^'Cih^'-^v\\%, + Co\\v - VhWh. 

Next, we let w ^ Vod _ffP(ri) (p > 2) be the solution to the following auxiliary 
problem: 

{{G',[u'])* {w),z) ^{v~ Vh, z) Vz e Vo. 
By assumption [A2] , there exists such a solution w with 

(4.19) \\w\\hp <CR\\v-Vhh-■ 
We then have for any Wh £ Vq 

\\v-v,\\l.^{iG'Au^]r{w),iv~VH)) 
= {G',[u-]iv~VH),w) 
— (v — Vh, w) 

= &"(« ~Vh,W- Wh) 

< CC2\\v - t)/i|j^f2||w - Wft|l_ff2. 



Consequently from (4.1) and (4.19) 



\\v - Vh\\l2 < CC2hJ'-^\\v ~ vh\\hAM\h^ 

< CC2CRh''^^\\v - vhWh^Wv - v/i|1l2, 



and thus, 
(4.20) 



\\V - Vhh- < CC2CRh^-^\\v ^ Vh\\H- 



Applying the inequality (4.20) into (4.18) gives us 
C4v-Vh\\'H2<CC^'Cih'''^\ 



Thus, for h < ho 



and therefore 



< CC^'Cih''-^\\v\\l, + CCoClClh'^-^Wv - vhfn.. 



VhU- <CC^'C2h 



e-2\ 



\\v~Vh\\L^ <CC^^ClCRh 



So far, we have been under the assumption that there exists a solution Vh- We 
now consider the question of existence and uniqueness. First, since the problem 
under consideration is linear and in a finite dimensional setting, existence and 
uniqueness are equivalent. Now suppose (p = 0, ip = 0. In light of (4.11), we have 
w = 0, and therefore, (4.16) implies = as well provided that h is sufficiently 
small. In particular, this means that (4.14) has a unique solution for h < ho. □ 



Remark 4.4. (a) Because (4.4) is a fourth order problem, we expect p > 3. There- 
fore, since the polynomial degree k is strictly greater than four, we expect r — p in 
Theorem 1431 
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(b) In many cases, it is possible to get a relatively good idea of how the constant 
Cr depends on e. To see this, suppose that there exists a constant C2 > such 



that if V solves then e H-^{n) and 

\\{F'[u^]yiv)\\^_,<d2\\v\\H^. 

Here, denotes the adjoint operator of 

Now if p = 3 in [A2] , then 

((GiK])*(^;),A^;) = (^,A«), 

where (•,•) now denotes the dual pairing of i?o(ri) and H^^{il). Therefore, after 
integrating by parts 

e\\VAv\\h^({F'[u^]y{v),Av)~{^,Av) 

< [C2\\v\\m + Mh- 
Hence, by Poincare's inequality 

(4.21) ||VA«|U2 < Ce-'{c2\\v\\Hi + Ml^^ 
By the proof of Theorem |4.2[ it is apparent that 

\\v\\h2<CC3Ml-, 

and therefore 

(4.22) ||VAw|U2 < Ce-\C2C3 + 1)\Ml- < Cs-^C^CMl^. 

Furthermore, if is coercive on i?Q(il), that is, there exists a constant 

Ci > such that 



|Ai;| 



(4.23) 



then by ( |4.21 ) 
(4.24) 



In view of (4.22) or (4.24), we can expect that in the general case 

\\v\\H^<Ce-'d2C3M\L-, 



and if (4.23) holds 



Hence, for p = 3 we have Cr = Ce^^C2C3 in the general case andCfl = Ce-^C^^C2 
if (|42^ holds. 



Now we consider the case p — A, and for simplicity, we assume F'[u^] is self- 
adjoint. We then have 



(G:.K](i;),A\)-(^,AM, 
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and therefore 



< c(\Ml2 + \\F'[u']\U\v\\„2)\\^'v\\l2 
<c(i + C3\\F'[u']\\]Ml4A\\\l2, 



where we define 



:— max 

°° l<ij<n 



dF{u^) 



max 

l<i<n 



dFiif) 
dpi 



We then expect that in this case that 

\\v\\H^<CC3e-'\\F'[u']\\^ML2 
Hence, for p = 4 we have Cr = CC^e^^ 



dFiu") 
dz 



4.3. Convergence analysis of finite element approximation 

In this section, we give the main results of this chapter, where we establish 
existence and uniqueness, and derive error estimates for the finite element method 



(4.3). First, we define an operator Th : V^} i->- Vj* such that for a given Vh G V'', 



Th{vh) is the solution to the following linear problem: 
(4.25) a'{vh-Th{vh).Wh) 

= e(Aw,„ Awh) + {F{vh),Wh) - 



J dwh \ 
' dv I 



In view of Theorem 4.3 Th is well-defined provided that assumptions [Al]-[A2] 
hold and h < Hq. We note that the right-hand side of (4.25) is the residual of the 



finite element method (4.3), and therefore, any fixed point of T;i (i.e. T{vh) = Vh) is 



a solution to (4.3 1 and vice- versa. Our goal is to show that indeed, has a unique 



fixed point in a small neighborhood of u^. To this end, we define the following ball: 

B„(p) :- {vH e Vg^; \\l''u' - Vh\\m < p}, 

where the center of the ball I'^u'^ is the finite element interpolant of u'^ . 

For the continuation of this chapter, we let £ = min{s, k + I}, where we recall 
that k is the polynomial degree of the finite element space V'^ and s is defined in 
[Al]. The following lemma shows that the distance between the center of B;, and 
its image under Th is small. 



Lemma 4.5. Suppose assumptions [Al]-[A4] hold. Then for h < /io(e), 

(4.26) Wl'^u' - Th{l"u')\\^, < Ceh'-^u'Wn^, 

where 

Ce = Ce{e) = CC;^ \\u'\\Ymax{C;\cj Cr}. 
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Proof. To ease notation, set r'^ — X'^u'^ — . Using the definition of T/i(-) 
and the mean value theorem, we have for any z/j e Vq 

(4.27) a^(lV-r;,(lV),zJ 



an 



= e(Ar^Az;,) + (F'[yJ(r^),z„), 

where y/i = Z'^u'^ — 7r'^ for some 7 G [0, 1]. 

Setting z/i — X^u"^ — Th{X^u^) and making use of [A2]-[A4], we have 



+ c||u^||y||r^||H2||zV - r;,(iV)||^, + Co||xV - r;,(zV)||^„ 

and so by the Cauchy-Schwarz inequahty, 

(4.28) Ci||l"u^-r„(X''u^)||^, 

< C^^e^\\r'\\]j. + CC^^\\u'\\l\\r'\\]j, + Co||l''w^ ~ T^(I"m^)||^, 

< CC^^e'-^Wu'fyWu'f^, + Co||x"u^ - Tu{X''u')\\\,. 

Next, we let w G Vq n HP{il) {p > 2) be the solution to the following auxiliary 
problem: 

{{G',[u'])* {w),z) = {X'^u' - n{X''u'),z) \fz e Vo, 

with 

(4.29) \\w\\hp < Cn\\x''u' - r„(Z''w=)||^,. 
Then for any Zh G Vq'' we get 

||i"u^-r„(zV)||', 

= a(l''u'^ -rft(l''u^),w) 

= a(l''ii^ - r,,(Z''u^), - Zft) + £(Ar^ Az,,) + (F'[y,,](r"), z^) 
<CC2||xV-r^(lV)||^,||7«-Zft||^2 +e||Ar^|U2||Az;,|U2 
+ C||u"||y||r"||H2||zh||ff2. 



Taking z/j — X^w, we have from (4.6) and (4.29) 



< CCl(clh^''~"^\\X^u'' - Th{X^u')\\\^ + h 



Substituting this last bound into the inequality (4.28) we have 



It then follows that for h < ho 



X''u' - n{X''u')\\^, < CC, ^ (Ci ^ + C-'CR)h'-^u'\\Y\\u'\\ 
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which is the inequahty (4.26). The proof is complete. □ 

Lemma 4.6. Suppose assumptions [Al]-[A5] hold. Then there exists an hi — 
hi{e) > such that for h < min{ft,o, hi}, the operator Th is a contracting mapping 
in the hall B/j(po) with a contraction factor ^, that is 

\\Th{vh) - Th{wh)\\jj2 < - whWh^ '■^Vh,Wh eMhiPo), 

where 

Pq = min|5, CCf i-i(/i)min{Cf ,Co"^C^^}| , 
and hi is chosen such that 

hi^C (Cj;'i(/ii)max{Cj"5,c|Cfl}) . 

Proof. By the definition of T^, we have for any Vh,Wh G Mh{po), Zh e Vq , 

a^{Th{vh) - Th{wh),Zh) = a^ivh.Zh) - a^[wh,Zh) + e{A{wh - Vh),Azh) 

+ {F{wh)~F{vh),Zh) 

= {F'[u']ivh -Wh),z^,) + {F{wh) ~ F{vh), Zh). 

Using the mean value theorem, we obtain 

a'{Thivh) - Tn{wn).Zh) - (i^'KlK - wn), z^) + {F{wh) - F{vh), z^) 

= {[F'[u']-F'[yh]){vh-Wh),zu), 

where yh — Wh + "f{vh — wt) for some 7 S [0, 1]. Here, we have abused the notation 
of J//,,, defining it differently in two different proofs in this section. 
Using [A2] , [A5] , and the triangle inequality yields 

Ci\\Th{vh) ~ T,,(w,i)||^2 

< - F'[yh]\\yy,\\vh - Wh\\H2\\Th{vh) - Th{wh)\\jj2 

+ Co||7)i(wh) - Th{wh)\\'^2 
<L{h)\\u'-yh\\ 

+ CQ\\Th{vh) - Th{wh)\\^j^2 

< CL{h){h^-^\\u'\\H' + po)\\vh - Wh\\H4Th{vh) - Th{wh)\\^, 

II 1 1 2 

+ Co\\Th{vh) - Th{wh)\\^2- 

Thus, 

(4.30) Ci\\Th{vh) - Thiwh)\\l2 < Co\\Thivh) - T;,K)||', 

+ CCi'L\h){h^'-4u^\le + pl)\\Vh - Wh\\j,2. 



Next, employing a duality argument similar to the one used in Lemma |4.5[ we 
let w eVof] HP{n) {p > 2) satisfy 

{{G'^[u'])* (w), z) = {Tn{vu) - n(wh), z) Vz e Vo, 

with 



(4.31) \\w\\hp < CR\\nivh) - n{wh)\ 



L2- 
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2 



Then using the same methods as in Lemma |4.5[ we conclude 

\\Th{vh)-Th{wh)\^^^ < c(^Lih){h'-^\\u'\\Hi + po)\\vh - whWh^ 
+ C'2h''~^\\Thivh) - Th{wh)\\jj2^ \\w\\hp 
< CCR(^L{h){h'-^u'\\Hi + Po)\\vh - whWh- 

+ C2h''^^\\Thivh) - Th{wh)\\jj2^ \\Th{vh) - Th{wh) 

and therefore 

\\Th{vh) - n{wh)\\l, < c(clL\h){h^'-''\\u'\\j,, + pl)\\vh - WhW 

+ c^cih^^-^niv,,) ~ mw,,)\\i,). 

Using this last inequality in (4.301 gives us 

II 1 1 2 

Ci\\Th{vh) - Th{wh)\\fj2 

< c(L\h){C^^ + CoCl) {h^'-^u'Wl, + pl)\\vh - WhW^H. 

+ CoClClh^^'^-^\\Th{vh)~TH{wh)\\\,. 
Therefore, for h < Hq 

\\Thivh) - Th{wh)\\j^2 

<cc;h{h){c;'^ + Co'CB^){h'-^\\u'\\He + po)\\vh-wh\\m. 

It then fohows from the definition of po and hi that for h < min{ft,Oi 
\\Th{vh) - Th{wh)\\jj2 < ^\\vh ~ WhWrn- 



□ 



With these two lemmas in hand, we can now derive the main results of this 
chapter. 



Theorem 4.7. Under the same hypotheses of Lemma 4-6, there exists /12 — ^'2(2) > 
such that for h < mm{ho, ft.2}> there exists a locally unique solution to (4.3), where 
/i2 is chosen such that 

h2 = C (Ce\\u'\\H^ina^{s-\C\h{h2)uiax{c;\cl^CR}}y'' . 

Furthermore, there holds the following error estimate: 

(4.32) \\u'-ul\\H2<Crh'-^u^\\H^, 
with 

C7 = Cj{e) = CC;'^\\u'\\Yma.x{C;\cjcR}. 
Moreover, there exists h^ — /13(e) > such that for h < min{/io, /i2, /13} 

(4.33) \\u'-uI\\l2 < Cs[c2h'+"-*\\u'\\H> + L{h)Crh^'-^\\u'\\l,y 
where 

h3 = C{C76-^\\u'\\H^.)^' , Cs = CC7Cr, r = mm{p,k + l}. 
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Proof. Let pi :— 2Ceh^ and note that for h < /12, there holds 

Pi ^ Po- Thus for h < min{/ioi^2} and noting /12 ^ hi, we use Lemmas 4.5 and 



4.6 to conclude that for any Vh € B/i(pi), 



2" 



- vhWh^ 



Thivh)\ 



^ Pi , Pi 
< h — = pi. 



2 2 

Hence, T/j maps Mh{pi) into B/i(pi). Since T/i is continuous and a contraction 
mapping in Mh{pi), by Banach's Fixed Point Theorem |42] T/j has a unique fixed 
point it| G B;j(pi), which is the unique solution to (4.3). To derive the error 
estimate (4.32), we use the triangle inequality to obtain 

\\u'-uI\\h^ < \\u'-I''u^\h2 + \\I''u'-uI\\h2 

< Ch'^-^\\u^\Hi + Pi < C7h^-^\\u^\H^. 

To obtain the error estimate (4.33), we start with the error equation: 

(Ae^ ^Vh) + {F{u') - F{ul),vh) = V«,, G V^^ , 

where e"^ := u'^ — u\. Using the mean value theorem, we obtain 

(4.34) (Ae%At-^) + (F'[y^](e^),«^) =0 ^vhEV^, 

where ?//,, — u"^ — 76^ for some 7 € [0, 1]. Again, we have abused the notation of y/j, 
defining it differently in different proofs. 

Next, let w £ H'p{VI) n Vq (p > 2) be the solution to the following auxiliary 
problem: 

((G^K])*H,z) = (e^z) Vzel/o, 

with 

(4.35) Ikllffp <Cfl|le"|U2. 
Using (4.34), we then have for any Wh € 



(4.36) 



{{G',[u^]nw),e^) 
a^(e^, w) 

a'^(e% w-Wh) + e(Ae^ Aw^) + (F'[u^](e^), wh) 
a=(e^ w~Wh) + { {F'[u'] - F'[yn]) {e'),Wh) 



< CC2V\\hA\w - WhWm + \\F'[u' 
Then by (|432| for < /ig 



F'[yh 



wv 



\m\\'Wh\\H2- 



Therefore, setting Wh ~ 1-hW in (4.36), we have for h < min{ft,o, ^2, ft.3}, 



<CCn(C2h--^\\e'\\H2+Lih)\\e'\\l,)\\e'\ 



L2. 
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Thus, 



||e1U2 < CCH(C2h--^\\e'\\H2 + mWe'Wl,) 

< CC^CR{c2h'+^'-^u'\\He + Lih)Cjh^'-^\\u'\\l,y 



□ 



Remark 4.8. (a) Noting 2i-A > £ for £ > 4 and /c > 4, Theorem 4.7 requires p > 4 
to obtain optimal order error estimates in the i^-norm. This regularity condition 
is expected provided that the domain CI is smooth and solution u' is sufficiently 
regular. 

(b) If {G'^{v)) is coercive on Vq, that is Co = in the inequality (4.5), then 
Cr = ^ 

Ci = 0{e) in such cases, and therefore (4.32) reads 

'2| 



Y in the error bound (4.32). Furthermore, it is expected that 



^h\\H^- ^ '-'fc" "■ II'"" \\y\\U \\h(- 

(c) We note that the constants C2,Cj,Cs appeared in the error bounds of 
Theorem |4.7| all depend on some negative powers of e, which is expected. The 
dependence of C2,Cj, Cs on we derived are the worst-case scenarios, they are far 
from being sharp (in particular, in the 3-D case) althrough the proved convergence 
rates in h are optimal. In Section |6] we shall present a detailed numerical study 
about the sharpness of the dependence of the error bounds on e~^. Our numerical 
experiments suggest that the error bounds only grow in in some small power 
orders, which are considerably better than the theoretical estimates indicate. 



CHAPTER 5 



Mixed finite element approximations 



The goal of this chapter is to construct and analyze a family of Hermann- 
Miyoshi mixed finite element methods for general fully nonlinear second order 
problem (2.7)-(2.8) based on the vanishing moment method (2.9)-(2.11 13. The 
mixed formulation is based on rewriting (2.9) as a system of two second order 
PDEs by introducing an additional variable. By decoupling (2.9) as a system, we 
are able to approximate (2.9|~( 2.11 )3 using only C° finite elements, opposed to 

finite elements used in Chapter [Jj which can be computational expensive and 
complicated. 

We note that the theory of mixed finite element methods, such as Hermann- 
Miyoshi mixed methods, has been extensively developed in the seventies and eighties 
for biharmonic problems in two dimensions (cf. [22, 13j). It is straightforward to 
formulate these methods for the fourth order quasilinear PDE (2.9 1 in two and three 
dimensions. Although it is a simple task to define mixed finite element methods 
for problem ( 2.9 )-( 2.11 13, proving existence of solutions and obtaining convergence 
rates are quite difficult. As is now well-known, proving existence and deriving error 
estimates for mixed methods relies heavily on the so-called inf-sup condition, and 
naturally, this is the starting point in our analysis. However, due to the strong 
nonlinearity in (2.9), the inf-sup condition is not sufficient for our purposes, and 
therefore, we must look for other techniques to obtain existence, uniqueness, and 
error estimates. To this end, we use a combined fixed-point and linearization tech- 
nique that is in the same spirit as in the previous chapter. 

The chapter is organized as follows. In Section [5T| we define the mixed formu- 
lation of ( 2.9 )-( 2.11 13, and then define the Hermann-Miyoshi mixed finite element 
method based upon this formulation. We then make certain structure assumptions 
on the nonlinear differential operator F, which will be used frequently in the anal- 
ysis of the mixed finite element method. The assumptions are generally mild and 
are very similar to those in Chapter [4] In Section 5.2 we prove convergence results 
of the mixed finite element method for the linearized problem (4.7)-(4.9). In Sec- 



tion [5^ we obtain our main results, where we obtain existence and uniqueness for 
the proposed Hermann-Miyoshi mixed finite element method and also derive error 
estimates. 



5.1. Formulation of mixed finite element methods 

There are several popular mixed formulations for fourth order problems. How- 
ever, since the Hessian matrix appears in (2.9) in a nonlinear fashion, we cannot 
use Am^ as an additional variable. This observation then rules out the family of 
Ciarlet-Raviart mixed finite element methods. On the other hand, this observation 
motivates us to try Hermann-Miyoshi mixed elements which use tr"^ := D^u'^ as 
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an additional unknown, and so, in this chapter, we will only focus on developing 
Hermann-Miyoshi type mixed methods for problem (2.9 1-(2.11 )3. 

In addition to the notation introduced in Section |L2| we also define the follow- 
ing space notation: 



Wo : -- 



\dn 



0- 



Recall that we use Greek letters to represent tensor functions and Roman letters 
to represent scalar functions throughout the paper. 

To define the mixed variational formulation for problem ( 2.9 )-( 2.11 (3, we rewrite 
the PDE into a system of two second order equations as follows: 

(5.1) a^-D^u'^O, 

(5.2) eAtr(cr^) +F(c^^^i^) = 0, 



where F{a^,u^) is defined in (1.131. 
Testing (5.1) with fi G Wo, we get 



(5.3) 

where {ti(x), , 
X, and 



(CT^M) + (div(/i),Vu") = 



E 



dg\ 



i=i ^ - ' / an 

., r„_i(x)} denotes the standard basis of the tangent space to dVl at 



(o•^/^) = / (j' ■■ ^J. dx ^ Y] / cr; 



dx. 



Next, multiplying ( |5.2[ ) with w £ Qo and integrating over gives us 
(5.4) -£(div(cr''), Vw) + (F(o•^ M^), w) = 0. 



Based on ( 5.3 H( |5.4[ ), we define the mixed formulation of ( 2.9 )-( 2.11 13 as fol- 
lows: find (cr^,^^) &We X Qg such that 

(5.5) {a',K) + b{K,u^) = Gin) yneWo, 

(5.6) 6(cr^^;) -e"lc(CT^u^^;) = Vw G Qo, 
where for e M^, v,w <E Q 

b{iJ.,v) : = (div(^), Vw), c{fi,w,v) 



(5.7) 



G(m) : = 



E 



dn/ 



an 



Next, let 7^1 be a quasiuniform triangular or quadrilateral partition of Q if 
n — 2, and tetrahedral or hexahedra mesh if rt = 3 parameterized by /i S (0, 1). 
Let Q'^ C Q he the Lagrange finite element space consisting of globally continuous, 
piecewise polynomials of degree k (> 2) associated with the mesh Th- 

We then define the following finite element spaces: 



Qo 



O'^ngo, 
[Q''] n Wo, 



^^-1 n xn 
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and define the norms ||(-, ■)\\^ , |||(-, ■)\\\^ : W x Q i-^ R+ such that for any {fi,v) e 
WxQ, 

\Kt^,v)\l: = hM^,+\\{^,,v)l, 

and Ki is defined by [B2] below. 

Based on (5.5)-(5.6), we define the Herman-Miyoshi-type mixed finite element 
method as follows: find {al,u%) G x Qg such that 

(5.8) ialK^,) + b{Kh,ul)^G{Kh) yKheW^\ 

(5.9) biaf^,Zh) - cial ul,Zh) = Vz„ G 

The main goal of this chapter is to prove existence and uniqueness for problem 
(5.8)-(5.9) and to also derive error estimates for — a\ and u*^ — u\. As a first 
step, we state the following inf-sup condition for the finite element pair (VFJjQq). 
The proof can be found in 



Lemma 5.1. For every Wh G Qq, there exists C > independent of h, such that 
(5.10) sup - — >C\\wh\\Hi- 



Remark 5.2. By [351 Proposition 1], Lemma 5.1 implies that there exists a linear 
operator ll'* : i->- W'^ such that 

(5.11) 6(M-nV,w'/0 =0 Vu-^gQJ}, 
and for ;u G n [iJ'*(r2)]"''" , s > 1, there holds 

(5.12) ll/z- n''/i||^^. < C/i^-^ |Ai|i/^ J =0,1, 1 < ^ < min{s,fc + 1}. 

Next, we assume the following structure conditions on the nonlinear differential 
operator F, which play an important role in our analysis. 

Assumption (B) 

[Bl] There exists Eq > such that for all e G (O,eo]i there exists a locally 



unique solution to (pj|-( |2.11[ )3 with G H^fl) (s > 3). 
[B2] The operator {G'^[a^ ,u'^]Y (the adjoint of G^[(t^,w^]) is an isomorphism 
from H^{n)nH^{n)to {H^{n)nH^in))*. That is for any G {H^{n)n 
H^{n))*, there exists v G H^{n) n H^i^l) such that 

(5.13) (^{G'^[cr^,u^]Y{D^v,v),w'j = {ip,w) Vw e H^{n) D H^in). 

Furthermore, there exists positive constants Kq = Ko{e), Ki = Ki{e), 
such that the following Garding inequality holdtj^ 

(5.14) {F'[a',u']{D%,v),v)>Ki\\v\\l,-Ko\\v\\l2 G Qo, 

where (•,•) denotes the dual pairing of Qo and Qq. Also, there exists 
K2 = K2 (e) > such that 



See Remark 



5.3 



d) for an interpretation. 
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where 



QQ" 



sup 

veQo 



sup sup 

v&Qo weQo 



\F'[<j',u']{D^v,v)\\^_, 



\\v\\H^\w\\m 

Moreover, there exists p > 3 and K]^^ > 0, (m = 0, 1) such that if 
ip G if-'"(rj) and V eVo satisfies ( [5l3l )"'then v e and 

[B3] There exists Banach spaces X, Y with a functional 
\\i;-)\\x>.Y:XxY^Ii+, 
and a constant C > such that for aU cj e X, y £ Y, x G W^, v £ Q 
\\F'[u;,y]{x,v)\\^_, < C\\{uj,y)\\^^^{\\x\\L^ + Mh^)- 



Furthermore , ( • , • ) 



XxY 



is well-defined and finite on W'^ x Q'*. 



[B4] There exists a constant — -f4'3(e) > such that 



XxY 



<K^{e) V7e[0,l]. 



where I^u"^ £ Qg is the finite element interpolant of u^. 

[B5] There exists a constant 6 = (5(e) £ (0,1), such that for any {iih,Vh) £ 
X with |||(n''- cr^ - tih^I^W" - Vh)\\\^ < 6, there holds V(k,„z,j) £ 
W'' X Q'' 



sup 



'F'[cr^,u^] - F'[nh,Vh]){Kh,Zh),Wh 



< Rih){\\a' - fihh^ + \\u'- vhWh^) \\\inh,Zh)\\\^ , 
where R{h) — R{e, h) may depend on e and h and R{h) = o{h'^~^). 
[B6] There exists Kq = Kq{£) > and a > such that for any 
{Xh,Vh) £ Tft := {iKh,Zh) £ Wo x Qo; {Kh,Xh) + b{xh, Zh) = Vxh e T^o'}, 
there holds 

\\F'[a',u']{xh-D\H,0)\\H-i <KGh''\\\{Xh,Vh)\l. 



Remark 5.3. (a) We made an effort in our presentation to state assumptions in 
this section that resemble those in the previous chapter, where conforming finite 



element methods for (2.9|-(2.11)i were studied. It is clear that conditions [Bl]- 
[B6] are similar, but slightly stronger than conditions [A1]-[A5]. For example, the 
inequality (5.14) suggests that the operator — _F"[m^] is uniformly elliptic, which 



rules out degenerate problems. However, assumptions [B1]-[B6] are still not very 
restrictive, and we will show in Chapter [6] that many well-known fully nonlinear 
second order differential operators satisfy these requirements. We also show a simple 
trick at the end of the chapter which makes it possible to incorporate degenerate 
elliptic PDEs (i.e. Ki—0) into the theory. 
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(b) We note that by definition of ct"^, G^, and F' (see Section 1.2) 

G',[a'',u%D^v,v) = G',[D'^u',u']{D'^v,v) = G'^[u%v), 
F'[a' ,u%D'^v,v) ^ F'[D'^u^ ,u^]{D'^v,v) ^ F'[u'']{v). 

It then seems redundant to write G'J^a'^ ,u'^]{D'^v,v) and F'[cr^, it^](£'^u, u) instead 
of G'^{u^\{y) and However, this (longer) short-hand notation naturally 

fits into the mixed method framework, and makes the subsequent analysis easier to 
follow. 

(c) It is obvious that assumption [Bl] is needed, and this assumption is actually 
the same as [Al]; we include it again for consistency and standardization. 

(d) Assumption [B2] is a natural extension of [A2], and F is expected to satisfy 
these conditions provided that —F is uniformly elliptic at u^, and 90 is sufficiently 
regular. We note that (5.14) needs to be understood with care because of the special 
notation we use. The left-hand side should be understood in the distributional 
sense. To derive the inequality, an integration by parts must be used on the second 
order derivative term. Also, Remark |4.8| gives heuristic estimates for the constants 
ifi^o and Kp>^ in terms of e. 

(e) By the standard interpolation theory and (5.12), we have 



h- 



lL2 



(5.15) 

(5.16) /i"^ + -I'^ii'^ll^i <Ch''-^\\u''\\Hi- 

(f) Condition [B5], which is used in the fixed-point argument, states that F' is 
in some sense locally Lipschitz near {a'^,u'^). 

(g) Clearly, if (k, z) e Wq x Qq satisfy 

(k,x) + &(x,^) = VxeWo, 

then D^z = k in a weak sense. However, if {Kh,Zh) G T/j, the analogous equality 
D^Zh = Kh is not necessarily true. Assumption [B6] indicates that the discrepancy 
between Kh and D^Zh under the image of F'[a^,u'^] is small. However, in what 
follows, we show that assumption [B6] holds with a = 1 if F is sufficiently smooth 
at the solution (cr^,u^). 



Proposition 5.4. Suppose 

Then assumption [B6] holds with a — 1 and 

dF{a^,u^) 



i,j = l,2,...,n. 



max 

l<2,j<n 



Or, 



A- max 

1< i,j<n 



Kg = G 

Proof. For any z G Qq, define such that 



_ dF{a-,u-) 



<III(X.,^'.)IIU|A^ 



ai^(c^^u^) 























h \e\ 
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Next, by (5.12), the definition of A"^, and the product rule, we have 



\\X 

< Ch 



V2;||/^2 max 

l<i.j<n 



dri 



\z\\i^6 max 

l<i,j<-n 



dF{a'',u'') 






W '5/ 



Therefore, by Poincare's inequality and a Sobolev inequality 
||A^-n"A^|U2 



< Ch 



max 

l<i.j<-n 



dr. 



max 

l<2,_7<n 



aF(cr^^i^) 






W •i/ 



The result follows from the above inequality. □ 

5.2. Linearization and its mixed finite element approximations 

To derive existence, uniqueness, and the desired error estimates for the mixed 
finite element method (5.8)-(5.9l, we must first study the mixed finite element 
approximations of (4.7)-(4.9), but with an alternative boundary condition: 

(5.17) G;K](w) = ^ mVL, 

(5.18) v = Q on dn, 

(5.19) D'^vvv = Q on do., 

where ip G Qq is some given function. Using arguments similar to the proof of 
we conclude that there exists a unique solution v E H^{n) D -ffo(f^) 



4.2 



Theorem 
to ( |5l7| -( |5.19| ). 

To introduce a mixed formulation for (5.17|-(5.19), we rewrite the fourth order 
PDE (5.17) as the following system of two second order PDEs: 

(5.20) x-D'^v^O, 

(5.21) eAtr(x) + F'[c^^^i^](L>2^;,l;) = ip, 

where tr(x) denotes the trace of x- 

The mixed variational formulation of ( 5.17 )-( 5.19) is then defined as follows: 
find (x, v) e Wq X Qo such that 

(5.22) (x,m) + Km,«) = V/ieW^o, 

(5.23) b{x,w) - £~^d{u'';v,w) = -e~^{p,w) Vw G Qo, 
where for v^w € Q 

d{u^;v,w) := {F'[a^,u^]{D^v,v),w) 

Remark 5.5. We note again that the right-hand side of c?(-; •, •) should be under- 
stood in the distributional sense. 

5.2.1. Mixed finite element approximation of linearized problem. Based 
on the variational formulation ( 5.22 )-( 5.23), we define the mixed finite element 
method for (5.17)-(5.19) as seeking {xh,Wh) G W}} x Qq such that 



(5.24) {Xh,^ih) + h{^ih,Vh)={) 

(5.25) b{xh: Wh) - e^^diu^;vh, Wh) = 



Our objective in this section is to prove existence and uniqueness for problem 



(5.24)- (5.25) and then to derive error estimates in various norms. 
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Theorem 5.6. Suppose assumptions [B1]-[B2] hold. Let v G H^iU) (s > 3) he the 
unique solution to (5.17|-(5.19l and x — D^v. Then there exists Hq = hQ{e) > 
such that for h < hg, there exists a unique solution {xh, Wh) G Wq x Qg to problem 
^1^-^^, where 



tfKo = 0, 



r = min{p, k + 1}. 
Furthermore, there hold the following error estimates: 



(5.26) 
(5.27) 

where 



\\\{X -Xh,v- Vh)\l < Ch'-\Kih + 1)\\v\\h^, 
\\v - vhWl- < K^h'+^-^K^h + l)\\v\\Hi. 



Ki = Cmax{Ki ^ K2e-^ Kr^e^}, K5 - CK, ^ K2KR^e^ , 
£ — inin{s, k + 1}. 



Proof. We first start by showing that the error estimates (5.26)-(5.27l hold 
in the case that there does exist a solution to ( 5.24 )-( 5.25 ). 

Let I'^v denote the standard finite element interpolant of v in Qq. Then using 



(5.11 1, we have for all {fih, Wh) € VFq x Q^, 

(5.28) {xh~l^'''x,l^h)+b{iih,vh~l'''v) = {x~'R''x,l^h)+h{^ih,v~l''v), 

(5.29) h{xh - n''x, Wh) - e-^d{u'- vh - X'^v, Wh) 

— e^^d{u^]I^v — V, Wh). 

Setting Hh = Xh — n''x and wt — vt —I'^v and subtracting ( |5.29[ ) from (5.28) 
yields 

(5.30) (x,. " n'^X, Xh ~ Tl^x) + e-'d{u'- Vh ~ I^, - T\) 

= (x - n"x, Xh - n'^x) + b{xh - n'^x, v - iS) 

Thus, by assumption [B2], 

||(x-n"x,«„-i^)||; 

< ||x - n'^xIL^ \\xh - n'^xIL^ + ||div(x,. - n''x)|L2 ||v(« -t^)||^. 



IL2 



< ||x - n'^xIL^ ||x/. - n'^xIL^ + Chr^Wxh - n'^xL^ ||v(z; -iS)||^, 

where we have used the inverse inequality in the last expression. 
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Using the Schwarz inequality, standard interpolation estimates, and rearranging 
terms, we have 



\{xh-'n\,vh-i''v)\\l 



|2 

1 1,2 



+ K^^Kle-^ \\v - I^v\^jj^ + Koe-^\\vh - T^v\\\^) 
which by an application of the triangle and inverse inequalities yields 



|2 



(5.31) \\\{^-^f^,v-Vh) 

< c{h^'-\Klh' + l)\\v\\l,+K^e-' \\v - vnt..). 

Continuing, we let w & Qof\ HP{fl) {p > 3) be the solution to the following 
auxiliary problem: 

{G',[u']Y{w) = v-Vh inn, 
D^wv ■ V = Q on dn. 

By assumption [B2], there exists such a solution and 

(5.32) \MH.<KRA\v~Vh\\L.. 

Setting K = D'^w G [i?^~^(ri)] it is easy to verify that {n,w) satisfy 

(k. At) + 6(/x, 2;) = V/xeW^o, 
b{K,, z) — e~^d*{u^; w, z) = e~^{v — Vh, z) \/z e Qo, 

where d*{u^; •, •) denotes the adjoint of d{u^; •, •), that is, 

d*{u^;v,w) = d{u^;w,v) \/v,w€Qo- 

We also note that there hold the following Galerkin orthogonality: 

(x - Xh, iJ-h) + KfJ-h, v-Vh) = o e w^, 

Kx - Xh, Wh) - e~'^d{u^; v-Vh,Wh)=0 Vwh e Qq. 
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Thus, choosing z — v — vu get 

V — Vh) + e^^d* [u^; w^v — Vh) 



6(k — n'*K, V — w/i) + e v — Vh, w) 

-b{\l''K,V-Vh) 

b(^K — U^K, V — I'^v) + e^^d{u^; v ~ Vh,w) 

+ (x-x/.,n"«:) 

+ (x - Xh, k) + (x - Xh, n'^K - k) 
- Kx - Xh, w) + {x~ Xh, n'^K - k) 

b{K - n'^K, V ~ I'^v) + e^^d{u^; v-Vh,w-l''' 
~-b{x-Xh,w-l''w) + {x-Xh,W'K- k). 



w) 



Therefore, using (5.32 1, 

-1 II l|2 

< ||div(K-n''K)|| Jv(t;-l''w)|U +K2e-^ 



\V - Vh\ 



\w — 1 wl 



+ ||div(x-X/i)||r2||V(w-X''w)|| 2 + Wx-Xh 



lL2 



L2 



< C{h 



Mh' + K2e-'h^-' llwllff. ||V(« - 



+ /i'^-ldiv(x-X/0 



< CK 



+ /i''"' llx - XhWm + h"'^ \\x - X/i|lL2) IK' - VhW 



L2 , 



and hence 



(5.33) \\v^vu\\l.<CKl^s\h^'+^^-'\\vr^,+ 



2,-2 i.2r-2 



\\^{v-Vh)\\ 



L2 



/^'■'-'llx-x.lll,i+/^''^-1lx-x/.lli.). 



Using estimate (5.331 in (5.31) yields 

\\\{X-Xh,v-Vh)\f^ 

< c{h^'-\Kie + l)\\v\\l. + s'^K.Wv - vnWh) 



< C 



Klh' + l)\\v\\]j, 



2,-2 i.2r-2 



\\V{v-Vh)\\ 



L2 



+ h^'-^\\x~Xh\\\^+h^'-^\\x-Xh\\l^])- 



62 



It then follows that for h < ho, 
\\\ix-Xh,v-Vh)\\\l 

< c{h^'-\Klh^ + l)\\v\\],,_ + KoKl^eh''+^^-'' 
and therefore 

\\\{x~Xh,v -Vh)\\\^ 

< c[h'-\Kih + 1) + KlKR^,e"^h'+'-^ \\v\\^, } 

< Ch^-^{Kih+l)\\v\\Hi, 



where we have used the fact that r > 3. Finally, (5.27) is obtained from (5.261 and 
( [5331 ). 

So far, we have been working under the assumption that there exists a solution 



However, using the Schatz's argument similar to the end of Theorem 4.3 



we can conclude from ( 5.26 )-( 5.27) that ( 5.24 )-( 5.25) has a unique solution for 
h<ho. □ 



Remark 5.7. (a) To obtain optimal order error estimates in the L^-norm (5.27), 
we require p > 4 and A; > 3 in the proof of Theorem |5.6[ 



(b) It is natural to ask why we use (5.21 ) instead of the alternative formulation 
(5.34) eMr{x)+F'[a',u']{x.v)^v. 



As it turns out, it is advantageous to use (5.21) opposed to (5.34), as we now 
explain. 

If we based the mixed finite element method on (5.34), the method would be 
similar to (5.24)-(5.25), but with d{u^;vh,Wh) replaced by d{u^;xh,Vh,Wh), where 

d{u'';xh,Vh,Wh) ■■= {F'[cr^ ,u'']{xh,Vh),Wh)■ 
Notice that by assumption [B2] (cf. ( |5.14 )) there holds 

d{u'';vh,Vh) > Ki\\vh\\Hi - Ko\\vh\\l2 \fvh E Qq. 

However, an inequality such as this one does not hold for d{u'^;xh,Vh,Vh) even if 
iXfnVfi) S T/i, where is defined in [B6]. 

However, if {xh,Vh) & T^, and if we define A"^ € Wo such that 

dF 



A: 



y - (cr^u")!)/, i,j = l,2,...n, 



then 



d{u^;Xh,Vh,Vh) = iXhA^) + {Fp[a'',u'']{xh,Vh),Vh) + {F^[a^ ,u'']{xh,Vh),Vh) 

+ {Fp[<j',u']ixh,Vh),Vh) + {F^[a',u^]{xh,Vh),Vh), 

and after integrating by parts 

d{u^;Xh,Vh,Vh) = iXh,>^^ ~ H'^A") + d{u'';vh,Vh). 

Thus, to obtain any coercivity from the alternative bilinear form •, •, •), we 

need to obtain bounds for \\X^ — H''A'^||i2, adding more complexity to the proof 
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of Theorem |5.6[ We also note the similarities of this derivation and the proof of 
Proposition |5.4| 



5.3. Convergence analysis of mixed finite element methods 

In this section, we give the main results of the chapter by proving there exists a 



unique solution to ( 5.8 )-( 5.9 ) and deriving error estimates of the numerical solution. 
First, we define the bilinear operator M^j : x h- > x such that for given 

{fih,Wh) e W,^ X Qi'g, Mh{fih,Wh) := {Mj^\fih,Wh),Mj^\tih,Wh)) G W^' x is 
the solution to the following problem: 

(5.35) {^J.h - Mjp (nhjWh), Kfi) +b[Kh,Wh- M]^\^ih.Wh)) 

= (m/i7 i^h) + b{Kh, Vh) - G{Kh) Vk/i e Wq, 

(5.36) b{nh - Mjp {fih,Wh), Zh) - e^'^d(u^ ;wh - Mj^^\^j.h,Wh), Zh) 



b{fih,Zh) 



^c{^J.h,Wh,Zh) 



yzh e g(j. 



By Theorem 5.6 M/i is well-defined provided assumptions [B1]-[B2] hold and 
h < Hq. Clearly any fixed point of the mapping M/j (i.e. M/i(/i/i, w^) — (iJ-hiWh)) 



is a solution to problem (5.8|-(5.9) and vice- versa. The goal of this section is to 



show that the mapping Mjj has a unique fixed point in a small neighborhood of 
(n^a^ ,l'^u^). To this end, we define the following sets: 

(5.37) §h{p): = {{f^h,Vh)eW^'xQ'^; - i-^^ iV) |||^ < p}, 

(5.38) Zh : = {{fih, Wh) e x Q^; {nh,Kh) + b{Kh,Wh) 

(5.39) Bh{p) ■.^Mp)riZh. 

For the continuation of the chapter, we set £ — min{s, k + 1}, where k is the 
polynomial degree of the finite element spaces and Q^, and s is defined in [Bl]. 
The next lemma shows that the distance between the center of B/j (p) and its image 
under the mapping M/j is small. 

Lemma 5.8. Suppose assumptions [B1]-[B4] hold. Then for h < Hq, there hold 
the following estimate: 



(5.40) 



whe 



Proof. To ease notation set ujh — H'^a^ 



M. 

jh—s 



(1) 



Sh 



M^^\u''a^,l''u^), r^ = I^u^ - ^i^ and ^ = n"cr^ - a^. By the definition of Mh, 
we have for any (k/j, Zh) e Wq x Qq 

{Ojh, Kh) + b{Kh, Sh) = {n'^a", Kh) + b{Kh,l''u') - G{Kh), 



biiOh, Zh) - e-^d{u'; Sh, z,,) = z,,) - e~^c{n''a\l''-' 



U ,Zh 
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and therefore by (5.5 1-( |5.6[ ), (5.11), and by employing the mean value theorem, 

(5.41) (Uh, Kh) + b{Kh, Sh) = (0": Kh) + b{Kh,r^), 

(5.42) b{uJh,Zh) - e^^d{u'';Sh,Zh) 

= b{0', Zh) - e-^ (c(^V^I''^i^ Zh) - c{a',u',zh)) 

where £,h — Il'V^ — 76*^ and yh = X^u^ — jr^ for some 7 £ [0, 1]. 

Setting Kh — i^h and Zh = Sh, and subtracting ( 5.42[ ) from ( |5.41 1 yield 

Consequently, by [B2]-[B4], and the inverse inequality, 

\\{<^h,Sh)\\l < \\e'\\L2\\uh\\L^ + ||div(c^^)||^,||Vr^||^, 
+ e''\\F'[Ch,yhW,r')\\^_, WshWhi+KoE-' 
< Wh^ + h-' \\uj,,\\^, WVr'W^, 

+ Ce-i||(a,yh)L,^(||01L. + ||r-|lHOP'^llHi+^o£-'||.s,||i. 

+ CK3e-W\9'h2 + \\r'\\Hi)\\sh\\Hi+Koe-'\\sh\\h. 
Using the Cauchy-Schwarz and inverse inequalities, and rearranging terms, give us 



(5.43) 



2 



iiiK,^^oiii:<c 

+ K^'Kie-'{\\e^\\l. 



L2 



< C h 



21-A ||^£||2 

\\(J 



1 7^2,-l7 2£-4|L ,e||2 



Next, we let w E QqD HP{fl) {p > 3) be the solution to the following auxiliary 
problem: 



in Q, 
on dil. 



with 

(5.44) \\w\\hv <KR,\\sh\\L2. 

Setting K = D^w G [iJP-^^^)] "^"^ j^^^g 

(k, ^) +6(^, z) = V^teW^o, 
6(k, z) — e^"'"(i*(u^; w, z) = — e^"'"(s;i, z) Vz e Qq. 
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^^""■"^ - -b{K,Sh)+e '^d*{u^;w,Sh) 



Thus, by ( |5ll] )-([5l2|, 

e~'\\sk\\h 

= -6(n''K,Sft) +e^^d{u'';sh,w) 

= iiJh,n) + {ujh,U''K-n) - (r,n''K) 

- 6(n''K, r") + e^^diu"- Sh,w) 
= -b{ujh,w) + [uh, n'*K - k) 

- , U'^k) - b{U''K, r') + e-^diu"; Sh, w) 
= -b{ujh, w - + (w/j, n'^K - k) 

+ e-i(i^'[a,y,.](r,r-),Z'^z«) 
< ||div(wh.)||^2||V(w-l''w)||^2 + ||cj,JU2||n'*K-K||^, 

'In'^Kll 



K2e ^\\sh\\Hi\\w -X'^w 



+ K:,e'W\e'\\L2 

< c{h'-^\\u:n\\L-^ + K2e-^h--^\\sn\\H^ + if3e"i/i^-2||^e||^^ ^ ||^||^^^ 
Therefore, using (5.44), 

II l|2 ^ riT^2 J2(u'ir-A\\ ||2 , t^2„-2 i 2r-2 ii ||2 , r^2 „-2 7 2£-4 n eii2 



Using this bound in (5.43), we have 

|2 



' 2f-4|i ,e||2 



It then fohows that for h < ho, 



which is the inequahty (5.40). The proof is complete. 



□ 



Lemma 5.9. Let [B1]-[B6] hold and suppose that S H^{n) (s > 3). Then there 
exists an hi = hi{e) > such that for h < mi-a{ho,hi}, the mapping M/j is a 
contracting mapping with a contracting factor of ^ in the ballMh{po), where 

po : = {KjR{h))-' 

hi : = mm {{RrKayi , {K^R{hi)\\umHi)^'} , 

and a > is defined in [B6]. That is, for all {Kh,Wh) G Mh{pa) 

1 



\Mh{^ih - Kh,Vh - Wh)\\\, < 



h - Kh,Vh - Wh)\\\^ 
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Proof. Let [nh.Wh) e B;i(po), and to ease notation we set 



(1) 



r(2)/ 



r(2)/ 



Using the definition of M/j and employing tlie mean value theorem, we conclude 
that for all (xh.Zh) € x gg, 

(5.45) {A4^\xh)+b{xh,A4^^)^0, 

(5.46) 6(M^'\z,,) -£-id(M^;Mf\z^) 

= e"^ (^^(m''; - lOh, Zh) - {c{fj.h, Vh, Zh) - c{Kh, Wh, Zh))^ 

= e^'^{d{u^]Vh - Wh,Zh) - {F'[^h-,Vh\{tJ-h ~ - Wh),Zh)^, 

where = ^.h + 7(^/1 - ^ih) and yn = Vh + 7(wh - w/j) for some 7 e [0, 1]. Here, we 
have abused the notation of and j/^, defining them differently in two different 
proofs in this section. 



Setting Xh = Mj^^ and Zh — M^^' , subtracting ( |5.45[ ) from ( [5.46^ , using as- 
sumptions [B2] and [B5], and the inverse inequality yields 

K'Mf)||[ 



f(2) 



< e 



L2 



- ((F'[a^ u'] {D^wh - Vh) - {tih - 0), Aff ))) + Koe''\\Mj^ 



(2)||2 
(2)11 



+ e ^R{h){\\<T'' -ChlU^ + Wu" -y/illffi) IIK^h - Hh,Wh -Vhjm^ ^^'^h llz/i 
<Ce-'(^KGh'' + R{h){h'-^\\u'\\H^+Po)) 

X \\\{Kh ~ iih,wh - Vh)\\l ||Aff 11^1 + /^oe^iA^f ^11^., 
and therefore 



(5.47) 



A4^A/f) 



< CK^h-^(Klh"- + R\h){e'-^\\u^\]j, + pI) 



X \\\{Kh- p.h,Wh-Vh)\\\^+ Kae \\M)^ 



i|lA,f(2)||2 



Next, we let z <E QqC\ H'p{VI) {p > 3) be the solution to the following auxiliary 
problem: 



(G^K])^(z) = Mf 
D'^ziy • 1/ = 



in n, 
on drt, 



with 



\\z\\H.<Kn,\\Ml^^\\^,. 



Letting x — D^z, we have 

{x,X) + b{X,z)=0 
b{x,y)-e-^d*{u'-z,y) = -e-\M'i^\y) 



VA e Wq, 



and hence by ( 5.45 )-( 5.46), 

-ii|m(2)ii^ 



-b{x.M'^^)+e-^d*{u-:z,M]^^) 



lL2 



(2)^ 



= (M«,n''x)+e-^d(w-;Aff\z) 

= {M\^\x)+e-'d{u';M';^\z) + (M«,n''x-x) 

+ e^^ (^(F'[£,h,yh]{Kh - iJ.h,Wh - Vh),l'^z^ - d[u^;wh - Vh,l'^z)^ 
< ||div(M^'))||^, \\V{z-l''z)\\^,+K2e-^M, 



(1)1 



IL2 



L2 

Ce- 



(2)11 llz-l'^zll 



X - fJ.h,Wh ~ Vh)\\\^ lll'^zl 



/i llifi ^ II /i 



r(i)| 



Kah'^ + R{h){h'-^u'\\Hi + po) 



Thus, 



(2)||2 < (^;^2 2 
h ||l2 - '-'-"-iJo^ 



-X2,-2;^2.-2||^(2)||2^^ 



(1)1 



/i2a ^ ^2(/^)(/^2£-4||^s||2^^ ^ ^2) |||(^^ _ ^^^^^^^ 



lL2 



Using the above bound in inequahty (5.47) yields for h < ho 

K\Mf)IIL 

RHh){h''-'\\urH^ + pif 



X \\\{kh - Ph,WH " Vh)\\\l + Xoe^iMf ^11', 

|2 



2l-A\\ e\\2 I 2 



X \\\{Kh ~ Ph,Wh - Vh) 

It then foUows from the definition of po that for h < min{/io, hi} 
||(M//\Mf)ll. 

< Ko(KGh" + R{h){h'-^u^\H^ + po)) |||(«;. - M?., u'h - vh) 
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Theorem 5.10. Under the same assumptions of Lemma\5.S\ there exists an /i2 = 



h2{£) > such that for h < min{/io, /12} ( |5.8| -(5.9 I has a locally unique solution, 
where /12 is chosen such that 



ho = 



{{KjKgT^ , {2K^KjR[h2)\\u'\\H^ 



Furthermore, there holds the following error estimate: 
(5.48) |||(a^ - al,u^ - ul)\\l < H'-^KsWu^Wh^, 

where 

Ks = CKe = CK^e-^k;-' + K^'Kr,). 

Proof. Let 

pi^2Keh'-^u^\H^. 

Then for h < mm{ho, /i2}, there holds pi < po. 

Thus noting ft,2 < hi, for any {ph,Vh) & B^(pi), we use Lemmas 5.8 and 5.9 
conclude that 



to 



<K^h'-'\\u^H^ 
Pi , Pi 



< 



Pi, 



2 2 

and so M^(/i;j, u/j) € Mh{pi). It is clear that M/i is a continuous mapping. It 
follows from Banach's Fixed Point Theorem '42\ that Mh has a unique fixed point 
((T^,u|) in the ball B/i(pi), which is the unique solution to (5.8|-(5.9|. 

To obtain the error estimate (5.481, we use the triangle inequality to conclude 



< 



< Ch^-^Wu-'WHi + Cpi < CKeh^''^\\u''\\Hi. 

□ 

Note that the error estimates of — //i in Theorem 5.10 are sub-optimal. 
In the next theorem, we employ a duality argument to improve the above error 
estimates and to also obtain error estimates. 

Theorem 5.11. In addition to the hypotheses of Theorem \5.1(^ suppose thatp > 4 
in assumption [B2] . Then there hold the following error estimates: 



HI < KnAK,h'-^+"^^'^i'''^^\u^\H. + KlRih)h 



2f-4||„,£||2 



U 



where 



Kg=CKsUlSi^{K2,KG}. 
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Proof. To ease notation, wc set 

7r= := cr" - al, := - u^. 

We note that by using the mean value theorem, there hold the following error 
equations: 

(5.49) (7r^M/.) + %h,e^) = V/i,, € W^o^ 

(5.50) b{Tr',Vh)-{F'[^h,yh]{Tr',e'),Vh) = Vv^gQ^, 

where = a'^ — 777*^, yh = u'^ — je^ for some 7 G [0, 1]. Again, wc have abused the 
notation of and yh, defining them differently in two separate proofs. 

Next, let Wm e HP-"'{fl)nQo (m = 0, l;p > 4) be the solution to the following 
auxiliary problem: 

{G',[u']y (wm) = (-l)"A™e^ in Q, 
D^Wmi^ ■ u = on dfl, 

with 



(5.51) < ||V"e^| 



L2 



Here, we have used the notation = A, = V, and A°, V° are the identity 
operators on Q. Setting K^n = D'^Wm G [7?^'~'"~^(ri)] we then have 

(k„,^) +6(^,w„) = V^gWo, 
b{Kra,v)-e-^d*{u'';Wra,v) = -£"1 (V™e^ V"?;) Vt; e go- 

Therefore, 

= (tT^, n^Km) + e~^d{u^; Wm) - b{Km - B'^Km, 6^) 

= (tt^, Km) + £~^d{u^; e^,Wm) 
= -b{n^, Wm) + e~^d(if] , Wm) 

- b{K„i - n^Kn, . - l''u') + (7^^ U'^Km " ^m) 

= -b{7r^, Wm - l''wm) + e^^d{u^] , Wm - T-^Wm) 

- b{Km - Ti'^Km, W - I'^u") + (tT^, H'^ Km - Km) 

+ e-^d{u'; e',l''wm) - e-^F'[^h, y^] (7r^ e'),l''wm) 

= -6(7r^, Wm - I'^Wm) + S"'^ d{u^ ] , Wm - I'^Wm) 

- b{Km - n'^K™, - I'^U") + (7^^ U'^Kr^ 

+ s-'({F'[a',u']-F'[^,„yh]){Tr',e'),X''Wm). 
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Bounding the right-hand side in the last expression, we have 



lL2 

< ||div(7r'')||^2||V(w„i + -ftr2e~^||e^||^i 

+ ||div(K„ - n'*K„)||^2||V(u'' + ||7r^||^2||n''K™ - 

+ (f' a'] {D^e' -Tr',0), I^w„, ) 
+ £- 1 ( (F' [a , 2/ J - i^' k^ ^1) ) , 

+ h\\7T^\L2+e-'KGh'^\\\iTr\e')\l 

+ -a^lU. + ||2/„. lll(^^e^)llls)lk™IU-" 

+ h'-'\\u^\\H^+h\\7r-\\L2+e-'KGh'^\\\{7r-,e')\l 
+ s-'R{h)\\\{7:^,e^)\\\l)\\V'-e^\\^, 

<CKn^e-'(^{K2h^-"^ + KGh'^) |||(^^ e^)|||, + |||(^^ e^)|||,') || V^e^l 
Therefore, 

||V™ei^, < [{K,h'"^ + Kah'^) |||«, e^)|||^ + R{h) |||(^^ e^) 

The proof is complete. □ 

5.4. Generalizations: the case of degenerate equations 

In this section, we generalize the analysis of the preceding sections to handle 
cases in which condition [B2] fails to hold, namely when the inequality 

(5.52) {F'[a',u']{x,v),v)>Ki\\vfH^-Ka\\v\\l2 V« e Qo 

does not hold for any positive constant Ki. Thus, in this section we consider 
cases in which the operator F may become degenerate (i.e. has vanishing smallest 
eigenvalue) at the solution u^. An instance of such a case arises when studying 
mixed finite element approximations of the infinity-Laplacian equation (cf. Section 



IL2 



6.3[ ). 

Here, we introduce a more flexible mixed finite element formulation to overcome 



this difficulty. To this end, we rewrite ( 2.9 )-( 2.11)3 into the following system of 
second order equations: 

(5.53) - D'^u" - rlr.y.nu" ^ 0, 

(5.54) ediv(div(CT^)) + ertr(5'^) + F(?^ w^) = 0, 
where 
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Inxn denotes the n x n identity matrix and r is a nonnegative constant that is 
independent of e. Clearly, ( 5.53 H( 5.54 1 is the same as (5.1|-(5.2| with = + 



The variational formulation of ( 5.53 1-( 5.54 ) is then defined as seeking (cr'^, u"^) e 
We X Qg such that 



(5.55) 
(5.56) 

where 

c{fj., V, w) 



(5^M) + 6(/A^i^) = G(m) y^ieWo, 

b{y,v)-e-^c{a',u',v)^0 Vw G Qo, 



{/i G W; jiv ■ V 
(div(/i), Vu) - r(tr(/i),i;), 
{F{fi,v),w) = 2eT{Vv,Vw) 



eriT' 



■{v,w) + rJ„xnW, 



and G(/i) is defined by ( |5.7[ ). We note that (5.55)-( 5.56[ ) is the same as (5.5)-( [5^ 
for the case r = 0. 

Based on the variational formulation ( 5.53 1-( '5.54[ ), we define our mixed finite 
element method of ([2!9|-( |2.11[ )3 as seeking (ct^,0 € x (where := 
W^f^W^ such that 



(5.57) 
(5.58) 



b{al, Vh) - £^^c(ct|, ul, Vh) = Vw/i G Qq. 



The specific goal of this section is to analyze the finite element method (5.571- 



( 5.58 1 and to determine what conditions are sufficient to show existence, uniqueness, 
and error estimates of the solution. Clearly, the finite element method and (5.8|- 
(5.9) have a similar structure, and therefore, one would expect that most of the 
analysis in the previous sections can be inherited in the present case. However, 
one issue of concern is that we have changed the bilinear form 6(-, •) in the new 
formulation, leading to question whether the inf-sup condition (cf. Lemma |5.1[ ) still 
holds. As is now well-known, this is a crucial ingredient in mixed finite element 
analysis, and we have used it copiously in the analysis above (albeit, indirectly). 
We appease these worries in the next lemma, showing that the inf-sup condition 
still holds provided r is small enough. The reason for using the new bilinear form 
&(•,•) will become clear later (see (5.68)). 

Lemma 5.12. There exists positive constants Tq,C depending only on n and 
such that for t < tq there holds the following inequality for any Vh G Qq : 



(5.59) 



sup - — >C\\vh\\Hi. 



Proof. By Poincare's inequality there exists a positive constant Cp that de- 
pends only on fl and n such that for all v G 11^(^1) 



\v\\l2 < C.p\\Vv\\L2. 
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For Vh e Qj} C H^{^), set kh = InxnVh e W^. Then 

b{^ih,Vh) ^ biKh,Vh) _ {dW{Kh),\/vh) - T{tr{Kh),Vh) 



sup 



Vh m 



> 



1 



iin{(l - Cp7iT),Cp} 



Choosing tq = |Cp ^n, we obtain the desired inequahty ( |5.59 ). □ 

Next, we introduce the analogous linearization problem and mixed formulation 
to ( |5.53| -( |5.54| . That is, instead of ( |5.20p -( |5.2l[ ), we write 

(5.60) X ~ — rlnxnV — infl, 

(5.61) ediv{dW{x)) +eTtr{x) + F'[y,u'']{D^v,v) = ip in ft, 

(5.62) 3^1/ . jy = 0, w = ondn, 
where we define 

F'[u!,y]{pi,w) : = -2£tAw - enr'^w + F'[uj - Tlnxny,y]{fJ',w), 
and F'[-, •](•, •) is defined by ^TJ^ . We note that (recaU a" = D^u") 
F'[a'',u'']{^i,w) = ~2eTA'w - enr'^w + F'la" ,u'']{^i,w). 



The variational formulation of ( 5.60 1-( 5.61 1 is then defined as seeking {x,v) G 
Wo X Qo such that 



V}) — e^^d{u^; v, w) = —e^^{^p, w) 

where 

d(w^;w,w) : = (^F'l^" ,u'']{D^v,v),w'^ 

= 2£t{Vv, Vw) - enr^(u, w) + {F'\a'' - 



Vm e Wo, 



,u^■^i^](i:»\,■^;),■^;) 



It then follows that the corresponding finite element method for the linearized 
problem is to find [xhiVh) G W^} x Qq such that 

(5.63) {Xh,t^h) + h{t^h.vu) =0 V^h e l^o^ 

(5.64) 6(x/i, w/i) - e"^d(M^;t;/i, = Vwh £ Qo- 

We now address what conditions are sufficient to show that the finite element 
methods ( 5.57 )-( 5.58) and (5.631 -(5.64 1 are well-posed. As it turns out, we are 
able to obtain results with weaker conditions than imposed in the previous section. 
Specifically, we are able to replace assumption [B2] by the following less-strict 
condition. 

[B2] The operator {G'^[u''])* (the adjoint of G'^[u^] defined in Chapter|4| is an 
isomorphism from H^{n) (1 H^{n) to {H^{fl) n H^i^l))* ■ That is for all 
ip e (-ff^(f^) n H^{n))*, there exists v e H^{n) n H^i^) such that 
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Furthermore, there exists a positive constant Kq = Kq^e) such that the 
following inequality holds: 

(5.65) {F'[a',u']{D^v,v),v) > ~Ko\\v\\l2. 
and there exists K2 > such that 

\\F'[a^u^]\\^^,<K,. 

Moreover, there exists p > 3 and Kj^g > 0, Kj^-^ > such that if (p £ 
(m = 0, 1) and v e Vq satisfies ( [sls] ), then v € RP-'^in) and 

< KliJ\(p\\H-rr.. 

Remark 5.13. We note that the only difference between [B2] and [B2] are the 
inequahties (5.651 and (5.14). Clearly if (5.14) holds, then ( |5.65[ ) holds as well, but 
not vice- versa. 

We now address the well-posedness of the finite element method for the lin- 
earized problem ( 5.63 )-( 5.64). 

Theorem 5.14. Suppose assumptions [Bl] and [B2] hold, r S (0, tq), v € _H'*(r2) (s > 
3) is the unique solution to (5.17|-(5.19) and x = D^v + TlnxnV- Then there 
exists an Hq — /10(e) > such that for h < ho, there exists a unique solution 
{Xh.Wh) G Wf} X Q'^ to problem where 

ho = (min { (KoKlKl^e-'r-') ^ , {K^Kl^e) ^ }) , r = min{p, k -f 1}. 

Here, k is the degree of the polynomial space of and , and p is defined in 
[B2] . Furthermore, there hold the following error estimates: 

(5.66) IIKx - Xh,v- Vh)\l < Ch'-'{K,h + 1)\\u^\h^, 

(5.67) \\v-Vh\\L2<K,h'+--\K4h+l)\\u'\\Hi, 
where 

ki^O (max{if2e"V-5,i^5^^^^i}^ ^ ^ O (k2Kr,t-"^^ , 

i — min{s, k + 1}, 



ll(M,'f)lle : = IImIIl^ +r5||t,||^i, 
\\\M\\l:^hMH^ + \U,v)\l 



Proof. It is clear from the proof of Theorem |5.6| that we only need to verify 
that condition [B2] holds, but with F'[cr^,M^] replaced by F'^^.u'^]. 
By the definition of u^], we have 

{F'W, u%D\, v),v)^ 2eT\\Wv\\l2 - enr^vWl^ + {F'[a', u^]{D\, v),v). 
Thus, if [B2] holds, then 

(5.68) {F'[a',u']iD\,v),v) > k,\\v\\l, ^ kolHl,, 

with 



euT 



Ko, 



Ki 2er. 
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We also notice that 

\\F ,u'\\\no' = sup sup r--^ 

2eT{Vv,Vw) - enT'^{v,w) + iF'la" ,u'']{D'^v,v),w) 
= sup sup r— r 7-- 

veQoweQo \M\m\\'^\\m 

<eT{2 + nT) + \\F'[a',u']{D\,v)\\^^, 

< £t(2 + nr) + K2 =: ^2- 

It then follows that [B2] holds but with F'[a'',u^] replaced by and 
the assertions of the theorem immediately follow. □ 

With the well-posedness results for the linear problem established, we can now 
state and prove the main result of this section (compare to Theorem 5.11). 

Theorem 5.15. Suppose assumptions [B1],[B2],[B3]-[B6] hold, u" G H^{n){s > 
3), R{h) = o{h?~^), T e (0, To), and there exists K3 = Ks{e) such that (5.73) 
holds. Then there exists hi — hi{e) > such that for h < min{/io, /ii}, there hold 
the following error estimates: 

(5.69) - ol^u- - ul)% < KsH'-^u^Ih^, 



(5.70) \\u'-uih2 < KnJK,h'-^+'^'<^^"^\\u^\H^+KiR{h)h 



2£-4||„,e||2 



U 



\H'- h 



(5.71) llu^ - < Kn, (K^h'-^+^'^^'^^^Wu^U^ + KlR{h)h''-^\\u'\\l,), 
where 

K9^CKsma.^{K2,KG}, 
£ = min{s, k + 1}, 

and s is defined in [Bl]. 

Proof. The idea of the proof is to show that [B2]-[B6] hold for F' (and 
F'[ct^m'^]) if [B2],[B3]-[B6] hold for F' (and F'[CT^w^]). The resuh then follows 
using the same techniques as those er nploy ed in the proof of Theorem |5.11[ 

First, from the proof of Theorem 5.14 we know that [B2] holds for F'[a^ , u^]. 
Next, if assumption [B3] holds then 

\\F'[^,y]{x,v)\\^_^ 

2er(Vu, Vz) - enT'^{v,z) + {F'[uj - T/„xny, y](x, w), z) 
= sup 

< £t(2 + nT)\\v\\m + C\\iLj - T/„x«y, y)\\xxYiML^ + Mm)- 
If we define 

(5.72) || (w - r/„x«2/, H^^y ' 
then 

\\P[uj,y]{x,v)\\^^,<er{2 + nT)\\v\\H^+C\\{u:,y)y^y{\\x\\L^ + \\v\\m). 



75 

From the definitions of Q'' and W'\ ||(-,-)lljfxy well-defined on W'^ x 
and if 

(5.73) \\{U^a' ~j<j',I^u' ~ju')y^^<K3 V7e[0,l], 

then it follows that conditions [B3]-[B4] hold for F' with || ('i ') II ifx? ™ place of 

ll(''')llxxF' _ 

Next, for any {iih,Vh) € Wl' x Q^^, (kh, zh) G 1^'' x g'', and G Q(} 



Define /x,- g Wg'' such that 

Mr := M/i ^ InxnVh, 

and notice that if 
then 

|||(n''cr^-/i^,l''u''-W,,)|||^^ <cs. 
Therefore, redefining 6 if necessary, we have 



< R{h){\\a' - ^lr\\L^ + 11"" - vhWh^) \\\{Kh,Zh)\\\^ \\wh\\m 

< CR{h){\\a' - ^ih\\L2 + \\u'~ v^Whi) \\\{nh, zh)\\l \\wh\\m. 

Hence, [B5] holds for 

Finally, we show that condition [B6] holds for F'[a'^,u^]. Suppose that 

(5.74) iXh^iih) + biK,„Vh) ^ Vk,, G W^, 

where {xh,Vh) G Wq x Qq. It then follows that 

iXh - TinxnVh, Kh) + b{Kh,Vh) = 0, 

that is (xh ^ TlnxnVfnVfi) G T;j, whcrc Tf^ is defined in [B6]. Thus, if [B6] holds 
(with Ki = er in definition of |||(-, ■)\\\^) and {xh,Vh) satisfies (5.74) then 

||F'[cr%1t^](X/i - TinxnVh ~ D^Vh, 0)||^_i 
<KGh''\\\{Xh~Tl„xnVh,Vh)\\\^ 

< KGh^[\\\{xh,Vh)\l + V^T{h\\vh\\H^ + |lz;,JU2)) 
<CKGh'^\\\iXh,v^)\l. 
Hence, F" fulfills all [B2]-[B6]. The proof is complete. □ 



CHAPTER 6 



Applications 



In the previous two chapters we have developed two abstract frameworks for 
conforming and mixed finite element approximations of the vanishing moment equa- 



tion (2.9) under some (mild) structure conditions on the nonlinear differential op- 
erator F. The goal of this chapter is to apply the two abstract frameworks to 
three specific nonlinear PDEs, namely, the Monge- Ampere equation, the equation 
of the prescribed Gauss curvature, and the infinity-Laplacian equation. These three 
equations are chosen because they represent three different scenarios categorized by 
their linearizations, which are respectively, coercive, indefinite, and degenerate. It 
is shown that the abstract frameworks of Chapter |4] and [5] are broad enough to 
cover all three scenarios. 



6.1. The Monge- Ampere equation 

The Monge- Ampere equation ( [l.lip is without question the best known fully 
nonlinear second order PDE. It is to fully nonlinear second order PDEs as the 
Poisson equation is to linear second order PDEs. The Monge-Ampere equation 
arises from applications in differential geometry, optimal transportation, geophysics, 
antenna design, and astrophysics. We refer the reader to [191 1421 144) and the 
references therein for more discussions about applications and PDE analysis of the 
Monge-Ampere equation. 

In this section, we consider finite element approximations of the Monge-Ampere 
equation with Dirichlet boundary condition: 

(6.1) det{D^u) = / (> 0) in Q, 

(6.2) u — g on Oil. 

A detailed analysis of conforming finite elements for the Monge-Ampere equa- 
tion was carried out in |39j (also see [61j). where the authors proved optimal error 
estimates in the energy norm. The authors also studied mixed finite element meth- 
ods for the Monge-Ampere equation in [38] (also see |61|) and obtained optimal 
error estimates for the scalar variable. However, we note that the results to be given 
below are sharper than those obtained in |38|, I39j in the sense that weaker regu- 
larities of the solution are required in the error estimates and the dependence 
on of the error bounds is less stringent. 

In the case of the Monge-Ampere equation, we have 

F{D'^u, Vu, u,x)^ f - dct(D^u), 

F'[v]{w) = -cof{D^v) : D^w, 
F'[ii,v]{k,w) — — cof(/i) : K 
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Remark 6.1. The inequality (4.5) implies that F{D^u,Vu,u,x) — f — det(£)^u) 
instead of det(D^M) — /, which is used in most PDE literature [42] . Recall that we 
assume —F is elliptic in the sense of |421 Chapter 17] in this book. 



The vanishing moment approximation ( |2.9[ )-(2.11 ) becomes 

(6.3) -eA'^u" + det{D^u^) = f in Q, 

(6.4) = g on d^, 

(6.5) Au^ on dn, 

and the linearization of 

Geiu") = eA^u" - det{D'^u') + f 

at the solution u"^ is 

G',[u']{v) = eA'^v - : D'^v = eA'^v - div($^V?;), 

where = cof(D^u^), the cofactor matrix of the Hessian D^u^ , and we have used 
Lemma 1.4 to obtain the last equality. 

6.1.1. Conforming finite element methods for the Monge- Ampere 



find u\ e Vg such that 



equation. The finite element method for (6.3 1-( 6.5 1 is defined as follows (cf. ( |4.3[ )): 

dvh ' 



(6.6) - e{Aul, Avh) + (det(i?^0, Vh) = (/, %) - ( e 



dv 



an 



Recall V ~ H^{il), and Vq and Vg'' are the finite element spaces of degree k > A 
defined by (|4^. 

The goal of this section is to apply the abstract framework of Chapter |4] toward 
the finite element method (6.6 1 in two and three dimensions. Namely, we verify 
[A1]-[A5] and determine how the constants, Ci, 5, and £(/i), depend on e. We 
summarize these results in the following theorem. 



Theorem 6.2. Let u'^ G H''{Q.) he the solution to (6.3 1 -(6.5) with s > 3 when 
n ^ 2 and s > 3 when rt = 3. Then for h < h2, there exists a unique solution 
£ Vj}' to (6.6 I. Furthermore, there hold the following error estimates: 



(6.7) 



where 



ul\\H2<Crh'-^\\u'\\H^, 

- ulU. < Cs[e-^^h'\\u'\\H^ + Crmh^'-^u'\\% 



C7 = 0(£^(l-2n))^ 



£ — min{s, k + 1}, 



and hy is chosen such that 



/i2 <c(e-5(i+2")|l7.|l^.L(/i2) 
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Proof. Wc first state the a priori bounds shown in Chapter [s] (also see j36j ): 
(6.9) Wu'Wh. = Oie"^) U = 1,2,3), \\u'\\w.^ = 0{e'-^) (.? = 1,2), 



Il2 



We also note that by interpolation between spaces, we have for p E [2, oo] 

(6.10) pViiiP < pVlll.p^u^iig ^ 0(£^). 

Next, since 

by standard theory for the biharmonic equation, if dfl is sufficiently smooth, then 
e H^{Q) with 

hlk^ <£-^(lldet(I?V)|U. + |l/|U.) 

<£-i(iii?vir2.„ + ii/iu.)- 



Therefore, in view of (6.10^, we have 
(6.11) 



(1 + 2^, 



|u||^f4 = 0[e 2 

Thus, by (6.9), (6.111, and interpolation of Sobolev spaces, we have 

(6.12) " " 

In addition, u'^ is strictly convex. Hence, $^ is positive definite, and therefore, 
there exists C > such that 



ff™ = 0{e"' 2" ") VTOe[2,4]. 

ctly conve 
hat 

($^Vw, Vw) > C\\Vw\\l2 Vw e Vo. 

It then follows that 

(6.13) {G',[u']iw),w)>Ce\\w\\l2 e Fq- 

Next, using a Sobolev inequality 



(6.14) 



= sup sup 



($^Vv, Vw) 



= sup sup 

\\^^\\l4^v\\l4^w\\l^ 

< sup sup 

veVoweVo \\''^\\h^\\w\\h^ 

If dft is sufficiently smooth and v E Vq solves 

{G4u%v),w) ^ {^,w) ^wEVo, 

where ip is some L^(r2)-function, then by standard elliptic PDE theory [321 I42j . 
V G HP{rt) for p > 2. Furthermore, in view of Remark 4.4 and the estimate 

dF 



we have 
(6.15) 



L2. 
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Thus, by ( |6.13| )- ( |6.15[ ), condition [A2] holds with 
(6.16) Co = 0, Ci=0(£), C2 = 0(£~5), 



and therefore (of. Theorems 4.2 and 4.3) 

(6.17) C3 = 0(e-^), C4 = 0(e-'), 

Cs = 0(e-^), /lo = 1. 

To confirm [A3]-[A4], we choose 

r = M/2.2(-i)(i]), Mlr-Mr;>,t.<-.-i)- 

For a smooth function we use Lemma |l.4| and a Sobolev inequality to conclude 



\F'[y]\\vv, ^ {coi{D^y):D^v,w) 

WvWy v(zVaw(^Vo \\y\\Y\\v\\HA\w\\H^ 

(cof (L»2y)Vi;, Vw) 

= sup sup TT— 7^ 



< c 



< c 



\cof{D^y)\\^, 
\ My ) \ 
It then follows from a simple density argument that 

\p'\y\\\vv' 



L2(,i-1) 



\\y\\ 



< c. 



sup 



-<C, 



yGY \\y\\Y 

and therefore condition [A3] holds, and by standard interpolation theory |22l 113] 
condition [A4] holds as well. 



We also note that by (|6.9|)-(|6.10| and Lemma 4.5 
(6.18) 



To verify [A5], we derive the following identity for any Vh G Vg'- 

((cof (ZJ^ye) _ cof (Ll^i;^)) VW, Vz 



\w\\m\\z\\H'^ 



< C\\coi{D^u') - cof{D^Vh)\\^3 



It follows that for n = 2, 

Hence, [A5] holds with L{h) = C. 

For the case n = 3, we conclude by the mean value theorem that for any 
hJ = 1,2,3, 

\\cof{D^u%j - co^{D^Vh)^J\\^^^ = \\det{D^u'\,,)-det{D\h\,j)\\^. 
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denotes the resulting 2x2 matrix after deleting the i row and j 



■th 



-f{D'^Vh\ij - D^m^Ijj)) for some 7 e [0, 1]. 



where D 

colmiin of D'^u^, and A'-' = cof (Z)'^■^i'^|^j 
Noting A*^ e R^^^, we have 

Thus, for any 6 G (0, 1) and Vh G with \\I^'u^ — Vfi\\ff2 < S, we have using the 
triangle inequality, the inverse inequality, and ( 6.10p 

\\F'[u'] - F'[vh]\\yy, < {Wu'Ww^.e + ||i;;,||,^2,6)||w" - VhWn^ 

1^2)11"" - VhWn^ 



< C{\\u^\\w^-e + h 

< C(£"5 +h-^6)\\u^ -Vh\\H^ 

<C{e-l +h-')\\u' -VhWm 
= L{h)\\u'' - Wft||//2. 



Thus, in the three-dimensional case [A5] holds with L{h) = C{e-i + h-^). 

Gathering all of our results, existence of a unique solution to ( |6.6[) a nd the error 
estimates (6.7)-(6.8) follow from Theorem 4.7 and the estimates ( |6.16 )-( [6.18[ ). □ 



Remark 6.3. (a) Estimates (6.7) and (6.8) give the same asymptotic rates in h as 
those obtained in |39| . However, they provide an improvement to these previous 
results in the sense that the constants C7 and /12 have a better order dependence 
in terms of e. 

(b) We require stronger reg ularity in the three-dimensional case to ensure 
L{h) ^o{h'^-^) (cf. Theorem [47|. 



6.1.2. Mixed finite element metiiods for the Monge-Ampere equa- 



tion. The mixed finite element method for (6.3|-(6.5| is defined as follows (cf. 
([5j8])-([5j|): find (c^^w^) € x Q^g such that 

K, Kh) + b{Kh, Ul) = G{Kh) Vk,, G WI^, 



h{al,zn)-e ^c{al,ul, Zh) = ^ 



It 



(6.19) 
(6.20) 
where 



G{Kh) is defined by ([5J|, Q = H^{^), W = [H^iU)] , and W}}, W^, Q(?, and 
are the Lagrange finite element spaces of degree k>2 defined in Section 
We now apply the abstract theory developed in Chapter [5] to the mixed 



b{Kh,u%) = (div(K/i), Vwl), c{al,ul,Zh) = (/ - det(cr^), z/^) , 



5.1 

finit' 



element method (6.19 )-(6.20 1. Similar to the previous subsection, our goal is to 



show that assumptions [B1]-[B6] hold, and to explicitly derive how the constants, 
Ki^ S, and R{h) depend on the parameter e. We summarize our findings in the 
following theorem. 



Theorem 6.4. Suppose € iJ*(f2) is the solution of (6.3 1 -(6.5) with s > 3 when 
n — 2 and s > 4 when n = 3. Furthermore, assume that k > A when n = 3. Then 



for h < h2, there exists a unique solution (ct|,u|) € x Qg to (6.19)-(6.20l 
Furthermore, there hold the following error estimates: 

(6.21) \\\{<j^ ~ alu' ~ u%)\l < Ksh'-^u^\H^ 



.22) 



'IWh^ <Kr 



K,h'- 



KlR{h)h 



2 
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where 



\v\\h^, 

70-5371 , 



\\\{p,v)\\\^ = h\\fi\\Hi + M\l-2 + e- 
if8 = 0(ei (22- 15"))^ Kg^O{e' 

R{h) = 0{\ log/il"^ +{n- 2){e-^h-^ + h-^)^ , 
£ = niin{s, k + 1}, 
and ft,2 is chosen such that 



ho 



r l+in ( 5 



(4-3n 



Proof. First, using the same arguments as those used to show condition [A2] 
in Theorem 6.2 we can also conclude that [B2] holds with 

K 



(6.23) 



0, i^i-O(l), 



K 



Ra 



and therefore (cf. Theorem 5.6 and Lemma [5?9| 

(6.24) Ki^O{e-^), if5 = 0(e~5), KT = 0{e~'^). 

To confirm [B3]-[B4], on noting that w](k, w) is independent of v and w, 
we choose the spaces X and Y as follows: 



X = 



2^(n-l)(n+£(3-n))^f^^ 
ri-1 



F = 0, 

Then using a Sobolev inequality, we have for all w e Jf, y G F, x G W, v € Q 

(cof(a;) : x, w) 



|J^'N,2/](X,^')L-1 



sup 



<C||cof(w)||^„+.<3-„,||xlL^ 

11x11 

<C||(c.,y)||^^y(||x||L^ + lklki). 

Thus condition [B3] holds. 

To confirm [B4], we note that by the inverse inequality, standard stability re- 
sults for the interpolation operator, and (5.12 ) to conclude that if cr^ S [iJ*~2(J7)j 
then for any p G [2, oo\ and £ G [3, min{s, k + 1}] 

(6.25) linVlliP < linV^-iV^lliP + IIi'^'t^IIlp 

< c{h^-"-+'-^\\a'\\H.-. + \\a'\\L. 
Therefore, for any 7 G [0, 1] 

II (^V-7a^IV- 7^^)11^,^ 
|n''cr^-7cr 



lL6(r.-l) 
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For the two-dimensional case, we set £ = 3 and use (6.9)-(6.10) to get 



II (nV- - 7a^IV - ^u-) 11^^^ < C{h^^ \\a'\\H^ + Wa^he) 

< C{hh-^ +e^i) = 0(e"^). 

For the three-dimensional case, we set £ = ^ and use (6.10 )-( 6.12| to get 



7^^)IL,^<c(ikii^t + ikiu-) 



Therefore by Lemma [5 
(6.26) K3 = o(ei 



(26-15n) 



Ke^O 



^£4 



(22-15n) 



To confirm [B5], we have for any (/i/j, v^) e x Q^', (k^, z^) e W^'' x Q'', and 
[F'la^^u"] - F'l^ih, Vh]) {nh, Zh) , w/i^ = (^(cof (cr'') - cof (^h)) : k^, w/j 

< qiog/l|^/ll-t||cof(f7^) -C0f(^;,)|U \\\{nh,ZH)\\l\\wh\ 



where we have used the inverse inequality |131 Lemma 4.9.1]. 

If n = 2, then ||cof (cr^) — cof(^/i)||^2 = 11'''^ ^ M/iIIl^, and so condition [B5] 
holds with R{h) = C\ log/i|i For n = 3, 



I (cof (cr") - C0f(^,,))ij| 



L2 



|det(CT''|ij) - det(/^h|ij)||^2 
|A''J' : {(j'Uj- ^lh\^J)\\L, 

|c \ij ~ M'llij 11^2 



<\W'\\l 



where A'-' = cof(cr'^| 
same notation as in Section 
<5e (0,1) 



7(/i/i| . . — (7^1 . .)) for some 7 e [0, 1], and we have used the 



6.1.1 



Since A'^ € R-^^^ we have for ||n''cr'^ - fj.h\\L2 < 



h „z I 



It then follows that [B5] holds in the case n = 3 with R{h) = C{s-^h-i + 
We note that for the hypotheses in Theorems 5.10 5.11 to hold, we require R{h) = 
o{h?~^) as /i 0+ for fixed e. This requirement is satisfied if £ > 2 in two 
dimension, and this bound is true provided ^ > 4 in three dimensions. 



Next, to verify condition [B6], we first use Holder's inequality and (6.9)-(6.10) 
to conclude that for p e [1, 2] and any i, j, fc = 1, 2, n {n = 2, 3) 
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Therefore, in view of Proposition |5.4| and the estimates 



dF 



dri 
dF 



U IIL° 



O e- 



to conclude that [B6] holds with 
(6.27) a 1, 



Kr=e- 



Thus, the existence of a unique solution (cr|, u|) to ( 6.19 )-( 6.20) and the error 
estimates (6.21)-(6.22) follows from Theorem 5.10| and the estimates (6.23)-(6.27). 

□ 

Remark 6.5. The error estimates in Theorem 16.41 have the same order of conver- 
gence in h as the estimates derived in [38], but the constants' dependence on e in 
(6.21)-(6.22) are sharper than these previous results. 

6.1.3. Numerical experiments and rates of convergence. Extensive nu- 
merical experiments for the finite element methods (6.6) and ( 6.19 )-( 6.20) in the 
two-dimensional setting have already been reported in 39 and 38j, respectively. 
These tests confirmed the error estimates (6.7)-(6.8) and (6.21 1-(6.22), and indi- 
cate that these estimates are sharp. Furthermore, the tests confirm the following 
rates of convergence: 

which are proved in Theorem |3.19| when the viscosity solution u belongs to the 

space w^^°°{n)nH^{n). 

In this section, we expand on these earlier results, performing two and three- 
dimensional numerical experiments and comparing the results with these earlier 
findings. We also show that for certain problems, one must choose an appropriate 
h — e relation in order for the method to converge. The tests below are done on the 
unit square n = (0, 1)" (n = 1, 2, 3). 

Test 6.1.1. For this test, we calculate — for fixed h — 0.02, while varying 
s in order to estimate ||m — We solve the mixed finite element method (6.191- 
(6.20) using the quadratic Lagrange finite element {k = 2) with the following test 
functions: 

(a) u = e^A+A+-l)/\ f={l + xl+xl+ 

[h) u^xl+xl+xl, / = 8. 

After having computed the solution, we list the errors in various norms in Table 
[T] and plot the results in Figures |2||3j The figures indicate that 



\u - 



\\u-uI\\hi = 0{ei) 
Therefore, since h is small, we expect 

\\u-u^\L^^O{e), 
\\u-u^\hi ~ Oiei] 



u 



hWL^ 



\u - u^||l2 



0(e), 



We note that these are the same rates of convergence found in ,39, and 
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Table 1. Test 6.1.1. Error of w.r.t. e {h ^ 0.02) 

E ||u - «^||i,oo(rato) u^||^2(ratc) Ijn - ti^ || „i (rate) || ct - cr'JI ^2 (rate) 

Test 6.1.1a 5.0E-01 1.19E-01( ) 5.71E-02( ) 3.47E-01( ) 3.34E+0G( ) 

2.5E-01 8.91E-02(0.42) 4.63E-02(0.30) 2.88E-01(0.27) 3.08E+00(G.12) 

l.OE-01 5.36E-02(0.55) 3.19E-G2(0.41) 2.09E-01(0.35) 2.72E+0G(G.14) 

5.0E-G2 2.35E-G2(1.19) 1.59E-G2(1.G0) 1.21E-01(0.79) 2.29E+GG(G.25) 

2.5E-G2 1.18E-G2(0.99) 8.95E-G3(0.83) 7.35E-02(0.72) 1.99E+GG(G.2G) 

1.0E-G2 5.57E-G3(0.82) 4.25E-G3(0.81) 3.91E-02(0.69) 1.66E+GG(G.2G) 

Test 6.1.1b 5.0E~G1 1.61E-G1( ) 7.47E-G2( ) 4.27E-01( ) 3.12E+GG( ) 

2.5E-G1 1.36E-G1(0.24) 6.48E-G2(0.21) 3.75E-01(0.19) 2.91E+GG(G.1G) 

l.OE-Gl 7.94E-G2(0.59) 4.17E-G2(0.48) 2.52E-01(0.43) 2.36E+GG(G.23) 

5.0E-G2 4.2GE-G2(0.92) 2.49E-G2(0.74) 1.61E-01(0.64) 1.92E+GG(G.29) 

2.5E-G2 1.99E-G2(1.G8) 1.36E-G2(0.88) 9.7GE-02(0.73) 1.57E+GG(G.29) 

l.OE-02 7.36E-G3(1.G9) 5.76E-G3(0.94) 4.85E-02(0.76) 1.26E+GG(G.24) 

5.0E-G3 3.79E-G3(0.96) 3.1GE-G3(0.89) 2.97E-02(0.71) 1.11E+GG(G.17) 



Test 6.1.2. The purpose of this test is to calculate the rate of convergence of 



— u^ll for fixed £=0.G01 in various norms. We solve problem (6.19 |-(6.2G I using 
the linear Lagrange element (k — 1), but with the boundary condition f 
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Test 6.1.1 




e 



Test 6.1.1 




Figure 2. Test 6.1.1. Error ||u — (top) and ||u — u|||i2 

(bottom) w.r.t. e (ft. = 0.02). 

replaced by cr|j^ • f^jg^ = 0^- We use the following test functions and data: 

{a) = xl + xj + xl, /" = 8, 

^xj + xl + xl 0^ = 2, 

(b) ^xf + xj+xl, f = 720x1x1 - e8640x^, 

g" = xl + xl + xl, (t)" = Uxjiyf + 2vl + iOxlvl- 

After computing the solution, we list the errors in Table [2] and plot the results 
in Figure |4] We note that the mixed finite element theory in the preceding sections 
was only developed for fc > 2. However, our numerical experiments also indicate 
that the method works for the case fc = 1. Indeed, the tests indicate the following 
rates of convergence: 

U - uih- = O {h') , \\u'^ uIWh^ ^0{h). 

Test 6.1.3. The purpose of this test is to calculate the error — || in various 



norms using a fixed h — e relation. We solve the finite element method (6.6) in two 
dimensions with denoting the Argyris finite element space of degree five 



and replace the boundary condition Am|| = e by Au|| = 0'^. We use the 
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Test 6.1.1 




Test 6.1.1 




Figure 3. Test 6.1.1. Error | 
(bottom) w.r.t. e {h ^ 0.02). 

Table 2. Test 6.1.2. Error of 
Lagrange element) 



Ih^ (top) and ||cr - al\\L2 



.r.t. h (£=0.001, linear 



Test 6.1.2a 



1.75E-01 
1.25E-01 
7.50E-02 
6.00E-02 
4.00E-02 
2.00E-02 



11" -«hlli2 

4.65E-02( ) 

2.25E-02(2.16) 
7.95E-03(2.03) 
5.13E-03(1.97) 
1.97E-03(2.35) 
1.13E-03(0.80) 



2.46E 
1.72E- 
1.04E- 
8.07E- 
5.28E- 
4.17E~ 



-01(— ) 
01(1.07) 
-01(0.99) 
-02(1.12) 
-02(1.05) 
-02(0.34) 



7.57E-01( ) 

8.75E-01(-0.43) 
8.39E-01(0.08) 
6.61E-01(1.07) 
5.85E-01(0.30) 
5.28E-01(0.15) 



Test 6.1.2b 1.75E-01 1.04E-01( ) 8.72E-01( ) 

1.25E-01 5.46E-02(1.92) 6.80E-01(0.74) 

7.50E-02 1.97E-02(1.99) 4.26E-01(0.92) 

6.00E-02 1.30E-02(1.85) 3.40E-01(1.01) 

4.00E-02 7.57E-03(1.34) 2.29E~01(0.97) 

2.00E-02 8.43E-03(-0.16) 1.85E-01(0.31) 



3.91E+00( ) 

3.92E+00(-0.01) 
3.75E+00(0.09) 
3.33E+00(0.53) 
3.25E+00(0.06) 
3.04E+00(0.09) 
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Test 6.1.2 




h 



Test 6.1.2 




10-^ 

h 



Test 6.1.2 

10^ p ^ ^ 

Test 6.1.2a 

g Test 6.1.2b 
U . * * 

a 



10' 



Figure 4. Test 6.1.2. Error - u|| 1^2 (top), (mid- 
dle), and \\a^ - a^L^ (bottom) w.r.t. h {e ^ 0.001). 

following test function and data: 



.6\2 



{x{ + x\)[%r^ - x{ - xi) + 8(r^) 
(r-(x?+^i))' 



/ = -^r^-(a;2+.T2), 
= 2 + e. 



2r^ - (xf + xi) 
{r^-{x\ + xl)y 



3 J 
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On the domain = (0,1)^, S C°^{Vl) for any e > 0. However, the hmiting 
function 



u := hm^ = -y 2 - (.T^ + x^) 

is not smooth, and in fact, there only holds u e W^'^ where p € [1,4) (cf. |29) ). 
We solve (6.6) using the following four h — e relations: 



ft. = 2e5, 
h = 0.5£5, 



h — e, 
h = 0.5ei. 



We list the errors of the computed solution in Table|3]and plot the results in Figures 

EHZl 

Since |jit^||//f — >■ oo for any ^ > 2 as e — >■ O"*", we suspect that a stringent h — e 
relation will be needed in order for the method to converge in view of the error 
estimates (6.7|-(6.8). This supposition is verified by the numerical tests, as the 
method does not converge in any norm using the relation h = 0.5ei . Furthermore, 
we observe that the method does not converge in the if^-norm for any h—e relations 
used in the experiments. This behavior is expected since the limiting solution u 
is not in this space. We plot the error of the computed solution in Figure [5] with 
parameters e = h = 0.04. As seen from the picture, the error is concentrated at 
the singularity of u. 



Computed Solution and Error 




Figure 5. Test 6.1.3. Error of computed solution with e = 0.04, 
h = 0.04. 



Test 6.1.4. For our last test, we numerically back up the theoretical results 
given in Chapterjsj that is, we compute the vanishing moment approximation (2.9|- 
(2.11 )i in the radial symmetric case. To this end, we solve (3.6)-(3.9) in the domain 
ri = (0, 1). We use the Hermite cubic finite element to construct our finite element 
space, and we use the following data: 

/ = (l + rV^'/', 5(1) = e^. 

It can be readily checked that the exact solution is u = e*" 

We plot the computed solution and corresponding error in Figure |8]with pa- 
rameters n — 4, e — 10^^, /i — 4.0 x 10^"^. We also plot the computed Laplacian, 
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Table 3. Test 6.1.3. Error of 11 — u? 11 with h — e relation 



2£2 2.00E-01 1.79E-01 3.94E-02 2.02E-02 9.41E-02 6.32E-01 

l.OOE-01 6.32E-02 4.11E-02 2.10E-02 l.OlE-01 7.55E-01 

5.00E-02 2.24E-02 3.45E-02 1.76E-02 8.85E-02 7.84E-01 

4.00E-02 1.60E-02 3.12E-02 1.59E-02 8.17E-02 8.17E-01 
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as well. As shown by the pictures, the vanishing moment 



methodology accurately captures the convex solution in higher dimensions. Also, 
as expected, the Laplacian of u'^ is strictly positive (cf. Theorem 3.9). 

Next, we plot both Ur and u^r in two and four dimensions in Figures [9f |10| 
with e-values, 10"^, 10"'^, 10"^. Recall that the Hessian matrix of u"^ only has two 
distinct eigenvalues u^^ and 
and for both dimensions n 



r'^'v 

= 2 



As seen in Figure |9[ u% is positive for all e-values 
and n — \. 



Corollary 3.7 



This result is in accordance with 
Finally, Figure 10 shows that is strictly positive except for a 
small e-neighborhood of the boundary, which agrees with the theoretical results 
established in Theorem 13. 121 



Test 6.1.3 




Figure 6. Test 6.1.3. Error -it|||L~ (top) and 
(bottom) with various h — e relations. 



Test 6.1.3 




Figure 7. Test 6.1.3. Error - (top) and - u|||//2 

(bottom) with various h — e relations. 




Figure 8. Test 6.1.4. Computed solution of (|3|6])-(|3|8]) (top), 
error (middle), and computed Laplacian (bottom) with n — 4,6 — 
10-\/i = 4x 10-3. 




Figure 9. Test 6.1.4. Computed of (|3!6|-(|3|8]) for n = 2 (top), 
and n = 4 (bottom) with e = 10""'^ (black), e = 10""^ (blue), and 
e = 10-5 (red) {h = 4 x IQ-^). 




Figure 10. Test 6.1.4. Computed w„ of (13J])-(13^ for n = 2 
(top), and n = 4 (bottom) with e = 10~^ (black), e = 10"^ (blue), 
and e = 10"^ (red) (ft. = 4 x lO'^). 
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6.2. The equation of prescribed Gauss curvature 

Let ft C R" be a bounded domain and g G C°(9ri). For a given constant 
K, > 0, the simplest version of the famous Minkowski problem (cf. |47L I42j ) asks 
to find a function u whose graph (as a manifold) has the constant Gauss curvature 
K. and u takes the boundary value g on dQ. The Gauss curvature of the graph of 
u is given by the formula 

det{D'^u) 

(1+ |Vu|2)^' 

and therefore, if such a function exists, it must satisfy 

(6.28) det{D'^u) = /C(l + jVup) in 17, 

(6.29) g on dn. 



The equation (6.28), which is called the equation of prescribed Gauss curvature, is 
a fully nonlinear Monge- Ampere-type equation. 

It is known [4? that there exists a constant /C* > such that for each IC E 
[0, /C*), problem (6.28 1-( 6.29) has a unique convex viscosity solution. Theoretically, 
it is very difficult to give an accurate estimate for the upper bound /C*. This then 
calls for help from accurate numerical methods. Indeed, the methodology and 
analysis of the vanishing moment method works very well for solving this problem 
and for estimating IC* . 

Unlike the Monge- Ampere equation considered in the previous section, we have 
some leeway in defining F{D'^u,'Vu,u,x). For reasons that will be evident later 
(cf. Remark 6.8), we set 

F{D'^u,Vu,u,x) = - 



(6.30) 
and therefore 



dct{D'^u) 

(l-h|Vu|2)^ 



F'[v]{w) = - 
F'[fj,,v]{K,w) = - 



coi{D^v) : D'^w det{D'^v)Vv ■ Vw 

+2 \- ( n + 2) 



(l + |Vf|2)- 

cof(/i) : K 



+ {n + 2) 



(l + |Vt;|2)^ 
det(/i)Vt; • Vw 
(1 + |V?;|2)^' 



Therefore, the vanishing moment approximation (2.9)-(2.11)i is 

det{D^u^) 



(6.31) 



(1 + \Vw 



e|2^ 



(6.32) 
(6.33) 

and the linearization of 

at the solution is 

G',[u'']{v) = eA^v - 



2 

= g 
Au' = e 



detiD^u") 



in n, 

on dil, 
on dil, 



(1 



r|VM^|2)^ (l + |Vu^|2)^ 
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where $^ denotes the cofactor matrix of D^u^ . 

Numerical tests indicate that there exists a unique strictly convex solution to 
( |6.31[ )-( |633l ) with e > (cf. Subsection [6.2.3[ and [371 H]). For the continu ation o f 
this section, we assume that there exists a unique strictly convex solution to (6.31 1~ 
(6.33 1. Furthermore, since the high-order terms in the equation of prescribed Gauss 
curvature are the same as the Monge- Ampere equation, we expect that the a priori 
bounds (6.9 )-( 6.10 1 hold for the solution it^ of the vanishing moment approximation 
(|63Tl)-(|6l3|. 



Before stating the finite element methods for ( 6.31 )-( 6.33) and applying the 
analysis of Chapters |4] and [s] to these methods, we first prove the following identity. 

Lemma 6.6. For all v,w £ Hoi^^) 



.34) 



{F'[u^]iv),w) = 



(1 + |Vm^|2)- 
Proof. Integrating by parts, we have 



{F'[u']{v),w) 



(1 + |Vw=|2)^ 



, Vw 



2 / -P^S/v ■ V(lVu^P) 



(1 + |Vw-|2)^ 



+ (n + 2) 



^det(D2M^)Vu^ • Wv 



(1 + |Vm^|2) 



Noting that (p^D^u" ^ det{D^u'')Iny,n, we conclude 

$^Vw • V(|Vu^p) = $^Vw • 2D'^u^Vu'' = 2 dct {D^u'')(yu'' ■ Vv) 



From this identity, (6.34) immediately follows. 

Since is strictly convex, we arrive at the following corollary. 
Corollary 6.7. There exists a constant C > such that 
(6.35) {F'[u']{w),w) > C\\w\\jJ^ Vw G H^{n). 



□ 



Remark 6.8. It is now obvious why we choose ( 6.30 ) as the definition of F opposed 
to the following choice: 



(6.36) 



F{D^u,Vu,u,x) = -det{D^u)+IC{l + IVmHt^ 



Indeed, if we chose (6.36) instead of (6.30) then 
F'[w^](w) = : D'^w + IC{n + 2)(1 
and a simple calculation shows 



JC{n 



■((1 + iVu^ntAw^ 



where A^cu'^ '■= D^u'^\7u'^ ■ Vu^j^ Since is strictly convex, both Au^ and AooU^ 
are positive terms, and therefore, the linearization of this choice of F is not coercive. 

We note 



^ Throughout this chapter we define A.caV := D^vX/v ■ Vti and Ao 



D-'vVvVv 



that both operators are referr ed to the infinity-Laplacian in the hterature [3l I34| . We shall use 
the latter definition in Section 6.3 
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Nevertheless, the above choice is also valid since it is easy to check that 
satisfies the following Garding inequality 

{F'[u']{w),w) > Cl\\w\\jJ^ - Co\\w\\l, ywGH^in), 

for some positive constants Co ~ Co{e), Ci = Ci{e). Here (•,•) denotes the dual 
pairing between i?o(51) and H^^{rt). In addition, other conditions of Assumption 
(A) also can be verified. We refer the reader to |61| for a detailed exposition. 



6.2.1. Conforming finite element methods for the equation of pre- 
scribed Gauss curvature. The finite element method for (6.31 1-(6.33 1 is to find 
w| e Vl^ such that for any Vh e Vq 



.37) -siAul,Avh) 



(1 



\Vul\^) 



,Vh 



I dvh \ 



The goal of this section is to apply the abstract framework of Chapter |4] toward 
the finite element method (6.371. Specifically, our goal is to show that assumptions 
[A1]-[A5] hold, and as a consequence, we will obtain existence and uniqueness of 
a solution to (6.37), as well as optimal order estimates for the error — We 
also pay close attention on the constants Ci and L(h) and how they depend on the 
parameter e. We summarize our results in the following theorem. 



Theorem 6.9. Let € H'^{il,) be the solution to (6.31) -(6.33) with s > 3 when 
n — 2 and s > 3 when n = 3. Then for h < h2, there exists a unique solution to 
(6.37). Furthermore, there exists positive constants C7, Cg such that 

\\u' ~ ul\\„2 < Crh'-^\ 



(6.38) 
(6.39) 
where 



\\u' ~ ulh2 < Cs{e-^^h'\\u'\\Hi + CrL{h)h 



2£-4| 



C7 = 0(e5(i-2"))^ 

and chosen such that 

h2<c(e-'^^'+^"^\\u'\\H^L{h2) 



Cs - 0(C7£-("+2)) 

£ = min{s, k + 1}, 



Proof. First, (6.35) implies that 
(6.40) {G'^[u']{v),v) > Ce\\v\\l, Vv e Vo, 

and it follows that (G'^[u^])* is an isomorphism from Vq to Vq . 

Next, for any v,w € Vq, using (6.9) and a Sobolev inequality, we have 



.41) 



{F'[u-]{v),w) = 



(1+ |VU^|2)^ 



< ll*'lliil|Vw||L«||Vw|lL6 

<c\\^'\\^i\\v\\h4M\h-, 
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and therefore, 



\\F'[u^]\\vv, = sup sup , , : ^ < C\m\ s. 



In view of Remark 14.41 and the estimates 



dF{u^) 
dpi 



e\2\ 



detpV)^/(l + |Vu^P)"*^ 



= {n + 2) 
< C\\Wu'\\l^ II det(L>2M=)||ioc ^ O 
we conclude that ii v E Vq is the solution to 

(6.42) ((G^KD* {^,w) £ Vo, 
for some (p E L^(r2), then 

(6.43) \\v\\hs < Ce-^ifih^ \\v\\h^ < Cs-'^^+'^Ml^ 



Thus, by ( |6.40[ )-( |6.43| ), [A2] holds with 
(6.44) Co = 0, Ci=0(e), C2 = 0(e-5), 



p = 4, CR = 0{e 



-{n+2)\ 



and therefore (cf. Theorem 4.3 1 

(6.45) C3 = 0{e-'), 

C5 = 0(e-(*+")), 

To confirm [A3]-[A4], we take 



C4 = 0(e-5), 
ho - 1. 



We then have the following bound for any v,z E Vq, y E Y: 

/ C0f(i^^y)V^ \ ^ ||,of(D2y)|| 3||V.||Le||Vz|Ua 

V(i + |vy|2)^ ; 

<C\\cof{D''y)\\^,\\v\\H44H- 
< Cp^yll'-L., IklkHklk- 

L 2 



It then follows that 



sup 



\F'[y]\ 



vv 



\\yh 



- <c. 



and thus, [A3]-[A4] holds. We also note from ([6J|)-( |6.10[ ) that 
(6.46) ||u^||y = o(e5(3-2n)^ ^ 



£3(1-2") 
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To verify condition [A5], we first make tlie following calculation for any w, z € 



(6.47) 



F'[vh]){w),z) 

" I (l + |Vzi^|2)^ " (l + |Vz;^|2)^' 



(<I>'^ - cof(D2w,,))Vw 



, Vz 



, Vz 



(l + |Vi;,|2)^ 

[1+ |VU^|2) 

Bounding the first term in ( |6.47p , we use a Sobolev inequality to conclude 



(l + |Vt;„|2)^ 



(6.48) 



(1 + |V«;,|2)^ 



,Vz < C\\<i>' ~coiiD'vh)\\^i\\w\\H2\\z\\H2. 



To bound the second term in (6.47), we first use the mean value theorem 



(i + |v«1^) 



1 + 2 



(1 + |V^;.P) 
-(n + 2)(l + |Vy^|2 



1 + 4 



where yn = + l^h for some 7 £ [0, 1]. Therefore, for any 5 £ (0, 1) and Vh € V^'* 
with Wl'^u" - VhWm < S 



-,Vz 

(1 + |VU-|2) ^ (l + |Vt;;,|2) ^ 

< C\\Vyh\\L4^iu' - Vh)\\L4'^'\\LA\^w\\L<^\\^4L0 



It then follows from this calculation and (6.48) that in the two-dimensional case 

\\F'[u'] - F'[vh]\\yy^ < Ce-i\\u' - v4h2 = L{h4u' ~ v4h2, 

that is, condition [A5] holds with L{h) — Ce~^. 

In the three-dimensional setting, using arguments similar to those for the 
Monge- Ampere equation, we have 



F'[u^]-F'[v,]\\<C{s-i+h-') 



v4h2 = L{h)\\u'' - w/i||_H-2, 



and therefore [A5] holds with L{h) = C(e"5 + h^^y 



Gathering up these results, and applying Theorem 4.7 with estimates (6.44|- 



(6.46), we conclude that there exists a unique solution to the finite element method 
(6.37) and that the error estimates ( 6.38 )-( 6.39) hold. □ 
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6.2.2. Mixed finite element methods for the equation of prescribed 
Gauss curvature. The mixed finite element method for ( 6.28 )-( 6.29) is defined 
as follows: find (cr|, m|) G x Q'^ such that 



(6.49) 
(6.50) 

where 



((t|, Kh) + b{Kh,ul) = G{Kh 

b{Kh,ul) = (div(Kft), Vu|), 



(1 



dct(a^) 

-|v<P) 



,Zh 



and G{Kh) is defined by (5.7). 

In this section, we apply the results of Chapter [5] to the mixed finite element 
method ( 6.49 )~( 6.50). Namely, we verify that conditions [B1]-[B6] hold, and from 
these results, we obtain existence and uniqueness of a solution to (6.49)^(6.50) as 
well as its error estimates. We summarize our findings in the following theorem. 



Theorem 6.10. Let it^ £ 77* (f2) be the solution to (6.31 )-( 6.33) with s > 3 when 
n — 2 and s > 5 when n — 3. Suppose fc > 3 when n = 2 and fc > 5 when 
n — 3. Then for h < h2, there exists a unique solution (cr^,M|) G x Qg to 



(6.49) -(6.49). Furthermore, there hold the following error estimates: 
\\\{a^-alu'-ul)\l<Ksh'-^u^\H^, 
\\u' - uI\\h. < Kn,{K,h'-^u^H' + KiR{h)h''-^\\u^l,)^ 



(6.51) 
(6.52) 

where 



lll(M,«)llle = ^Mh^ + y\\L^+£''\Mm, 

76-49T1 



if8 = 0(£i(28"19«)), 

logh\{h-^ 



Kg = 0{e- 



R{h) = C 
i = min{s, fc + 1} 



n = 2, 



and ft,2 is chosen such that 

h, = Cmin{e^, [s'^^^'-''-^ Rih,)\\u^\H.) , [e-'^ Rih,)\\u^y.)'^'}. 

Proof. First, using the same arguments as those used to show assumption 
[A2] in Theorem 6.9 we can conclude that [B2] holds with 

Ko = 0, Xi = 0(l), K2 = 0{e-^), 

p = 4, Kji„ = 0(e-("+2)), K^^ = 0(e-2), 



(6.53) 



and therefore, (cf. Theorem |5.6| and Lemma 5.9) 

_ 3 



(6.54) 



3(2"+5)^ 



O e-r 
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We now turn our attention to condition [B3]. To show that this condition holds, 
we set (n = 2, 3) 



X 



^n(n+e(3-n))(^^^) 



||(t^,2/)||xxy = ||w|r2i;.-i)(.+.(3-„), + ||w||2„(,.+.(3-„),|lVy|lL~ € X, y € y. 
Then for any a; g X, y £Y, x e VF, and v £ Q, z £ Qo, we have 



cof(w) : X 



(l + |Vy|2)^ 



+(n + 2) 



det(a;)Vt/ • Vu 

(l+|Vj/|2)^^ 



< C(^||cof(w)||i„+e(3-„)||x||L2||2||ffi + II det(a;)||i„+.(3-„)||Vy||L~||Vw||L2||z||/fi^ 
<C( 



n-l 

i(n-l)(„ + E(3- 

n-1 



.))IIxIIl2 + ||w||2„(,.+e(3-„))||Vy||L=o||v||Hl)||2| 

< c(iic.ir-i,,(„^,(3_„„ + iic.ii2„(„+e(3-„,)iiv2yiu»)(iixiiL2 + mihOii^iih- 

It follows from this calculation that 

\\F'[Lo,y]ix,v)\\^-^ < C\\ioj,y)\\x^Y{\\x\\L-^ + \\v\\h^), 

and therefore condition [B3] holds. 

To confirm [B4] , we use ( [6^ to conclude that if cr^ G [H'-^in)] , then for 
any 7 S [0, 1] and £ € [3, min{s, k + 1}] 



XxY 



< c 



(l|nV|r, 



in-l 

|x,(>.-l)(" + e(3-„)) 
1 1 1 1 



2;,n(,i + e(3-n)) 



i>i(n + e(3->i)) 



< c(h'-"^-^\\a'\\H^-2 + ||a^|L„(„+.(3-„„)' 



For the two-dimensional case, we set i — 3 and use (6.10) to get 



XxY 



For the three-dimensional case, we set € = and use (6.121 and (6. 10 1 to 



conclude 



< c 



XxY 



Combing these two estimates, we have 

(6.55) K^ = (ei(30-i9»)) , K, = O 
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As a first step to confirm [B5], we note that for all {^j.h,Vh) G x Qg and 

(6.56) ((F'[(7^ u^] - F'[pLh, Vh]) {^h, Zh),wh) 
( I Cof(/l/i) Cof(cr'') 



: Kh,Wh 



+ (n + 2) 



To bound the first term in (6.56), we use the mean value theorem to conclude 
/ / C0f(/ift) cof(cr^) 



_(l + |Vl;^|2)=^ (l + |Vu^|2)^ 

_ / (cof(^,J - C0f((T^)) : Kh 



(l + |Vi;„|2)=i 
^ I cof(cr'^) : Kh 

V(1 + |VU^|2)^- 

(C0f(^/J - C0f((T^)) : Kh 
(l + |Vi;;,|2)=*^ 



cof(o-^) : Kh 

(1 + |V«^|2)=^ 



,Wh 



,Wh\ -{n + 2) 



Vyft • V(w/i - ■u^)cof(cr^) : Kh 

i+4 ,Wh I , 

(l + |Vy,.P) 2 



where j/h = vj^ + ju'^ for some 7 G [0,1]. Therefore, by (6.9) and the inverse 
inequality, we have 



cof(/i/i) 



cof(cr^) 



< c(||cofO^;,)-cof(a^)|U2 + ||V2/,,||Loo||u^-i7^||Hi||cof(a^)||Lo 

X ||k/i||l2||w,,||l^ 

< C|log/i|T^/ii-t(||cof(^ft.)-cof(a^)||i2+£-i||V2/,,||L-||u^-t'ft||Hi) 

X \\wh\\m- 

If G Qg with — UhUffi < <5 G (0, 1), then by the inverse inequality 

W^yhh^ < IIVw^IIloo + ||VZV||l=o + |log/i|^/ii-t||V(lV-?;,0||Hi 
< C|log/i|^/ii-t. 

Furthermore, from the mixed finite element analysis for the Monge- Ampere equa- 
tion, we have 



llcof(a^) - C0f(/i/.)|U2 < C(£(2-") + _ ^^11 



L2. 
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Using these two inequalities in (6.57), we arrive at 
cof(/i/i) coi{cr^) 



(6.58) 



< C\ log/ll^/l^-t ((£(2-") + /if (2-"))||a- _ ^,11^, 

< C{e-'\ \ogh\f-"{e-'h-'2 + h-'y'-''{\\a' - fij^h^ + \\ 



u -Vh\\H^ 



Using a similar strategy to bound the second term in (6.56), we add and sub- 
tract terms and use the mean value theorem to conclude 



(6.59) 



det(<T^)Vu^ 



det(/i/i)Vu/i 



^ (1 + |V«,P)^ 



^(1 + |VU^|2) 

_ / (det(g^) - det(/x/,))Vu" ■ Vz/, det(/i/,)(Vu'^ - Vt;/,) -Vz,, 
\ (l + |Vz;^|2)=*^ ^ (n_|Vt;^|2)^ 

V(l + |Vu-|2)^ (l + |Vt;„|2)^ J 
( {coi{S,h) ■■ {cj' - t^h))yu' ■ det(/i?,)(Vw" - Vvh) ■ Vz,, 



Wh 



(1+|V«„|2)^ 



, ^ /det((T-')(Vw^- Vzh)(Va;,,- V(ii^-w,0) \ 

where = cr*^ + ^liih-, Xh = u'^ + j2Vh for some 71,72 G [0, 1]. Bounding the first 
term in (6.59), we use the inverse inequality to conclude 



,'Wh 



(cof(a) : (a^-M/»))V""-V^/. 
(l + |V«,,|2)-r 

< C||cof($ft)||i,2||cr^ - /i/i||L2||V-U^||L==||Vz,i||Loc||w,J|i^ 

<C\\ogh\^h^-"\\coi{^h)\\ 



If ||n (t'^ — /ihlU^ l£ 6 £ {0, 1), then by (6.10) and the inverse inequality 

n-l 



„ / 3-2n 



+ /ii(2-n)^ 



Therefore, 



(6.60) 



/ (cof(a) : ((7^-a^,0)Vm^- Vzfe 
1^ (l + |Vt;,.|2)=*^ 

< C|log/i|^/.i-"(e^ +/ii(2-^))||a^ -Ai,|U2 |||(/«,,z,)|||Jlu;;,|l^.. 
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Bounding the second term in ( |6.59 ), we have 

[ (i + |v..P)^ '"V 

< II det(/i/j)||i2||VM^ - Vw,J|L2||||Vzft||L~||wh.||L- 

< C\\ogh\^-^h^-^det{fih)\\L4u' - VhWmW^hWmlhhWm- 
If lln'^o''^ — M/iIIl^ < € (0, 1), then by (6.10) and the inverse inequahty, 

II dct(M,)|U. < c||a*.||2- < c{\\a'\\i.. + iinV^ - /.;,||2 



0{s"-^ + /ltd-")) ^ 0(e^ + /l-2n+3)^ 



Therefore, 



/ det(AtO {Vu' - Vvh) ■ Vzh ^ 
(6-61) ^ ^f^h 

V (i + |v«,p)^ , 

< C\\ogh\^-"h^-"(e'-^ 



we 1 

det(g-^)(VM^ ■ \7z!,){\/xh ■ \/{u' - Vh)) \ 
(l + |Vx.P)^ '"V 

< ||det(CT")||i~||Vu"||L»||Vz;,.|U2||Vx,J|i~||V(u^-i;,,)|U2||u;;,||i 

< C\\oghf-^ h^-'^Wa'Wl^WWxhh^Wu' - VhWm |||(«;,.,z/.)||U|u^h 
<C|log/i|'^/ii-te-"||Vx„|Uoe||u^-i;„||^i \\\{Kh,Zh)\l\\wh\ 
<C\\ogh\'-"h^-"e-''\\u' -VhWm |||K,z;,)||y|«;,||^i. 



Next, using similar arguments as above, we bound the third term in (6.591 as 
follows: 



(6.62) 



Applying the bounds (6.60)-(6.62) to (6.59), we have 

< c(|log/i|^/ii-"(e^ +;it(2-n)) 

+ |log/i|3-"/i2-"(£^ +/i-2»+3) + |log/i|3-"7i2-«e-«' 

X {\W - M/iIIl^ + l|w^ - Whllffi) IIK^h, Zh)\\\, \\Wh\\H^ 

< C\ l0g/ip-"(£-" + /l-2"+3) (11^- _ + \\u-- y^W^,) IIK^,, ^,)|||^ ||^,J 

Finally, we combine this last inequality with (6.56)- p.59[ ) to get 

< c(£-l|l0g/l|3-"/l2-„^ |iQg^|3-n^g-n_^/j-2n+3)^ 

X (Ik" - m/iIIl2 + l|w" - vhWh^) lll('«/i,2;/i)lll£ Ikhllffi 

< C\ l0g/lp-"(£-l/l2-n _^ ^-n _^ /i-2"+3) 

X (||ct^ - ^,,11^2 + Ijii"^ - vhWh^) lUiih^Zh)]]]^ WwhWn^- 
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Therefore, condition [B5] holds with R{h) = C\\ogh\{e^^ + h^^) in the two- 
dimensional case and R{h) — C[e^^h^^ + e^^ + h^^) in the three-dimensional 
case. 

To establish condition [B6] , we use similar arguments to that of the mixed finite 
element analysis of the Monge- Ampere equation to conclude 



OF 


< 


drij 





dF 



< C 



(1 + |VW^|2)^ 

(l + lVu^P)-^ 



< ||cof(f7^)||ioc, =0(e-i) 



|cof(a'^)|| 



< C\\u'\\w2.^\\cofia')\\^,,e = O . 



Therefore by Proposition 5.4 condition [B6] holds with 

a = Kg ^ Ce--^^^+"\ 

Wrapping things up, we apply Theorem |5.10| and |5.11| to obtain existence and 
uniqueness of a solution ((t|, u^) to the mixed finite element method (6.49|-(6.50). 
The error estimates (6.51 1-( 6.52) also follow from these results and by the defini- 
tions 

Ks = CKr, Kg - Ksmax{K2,KG}. 

□ 

6.2.3. Numerical experiments and rates of convergence. In this sec- 
tion, we provide several two-dimensional numerical experiments to gauge the effi- 
ciency of the finite element methods developed in the previous two subsections. 

Test 6.2.1. In this test, we fix /i = 0.01 in order to study the behavior of m^. 
Notably, we are interested whether ||m — it'^|| — as e 0+. To this end, we solve 
the following problem: find ii| G Vg such thatj^ 



_(1 + |Vm^|2) 

Here, we take V'^ to be the Argyris finite element space 



dv 



of degree k — 5 and 



set r2 = 


(0, 1)2. We use the following 


test function and parameters: 


(a) u — 


e 2 , 


/C 0.1, 


f = 




^1 + ^2 

5 = 62 . 


0.1(1 + (x2-Kx2)e^?+-i)2' 




cos(Y^7r) + cos(-y/a?27r). 


/C = 0.025, 


/ = 


TT^ (xi ^ sin(ya?]'7r) — Xi "'"tt cos(Y^7r)) ^ sin(^a;27r) — """tt cos(^a;27r)) 


16 0.025(H-^(a; 


I ^ sin2(y^7r) + X2^ sin2(y/a?27r)))^ 


9 = 


cos(Y^7r) -f- cos(-yiE27r). 





The computed solution, whose values are given in Table |4j is compared to the 
exact solution in Figure 



11 



As seen from Figure 



11 



the behavior of ||u — 



We note that it is easy to see the finite element methods and their convergence analyses of 



Section 6.2.1 and 6.2.2 also apply to the case / > but / ^ 1. 



107 



Table 4. Test 6.2.1: Error of 
estimated rate of convergence 



w.r.t. e {h = 0.01) and 



II" - «hlL2(rato) ti;;jl„i(ratc) ||u - m^|1^2 (rate) 



Test 6.2.1a 



l.OE-01 
5.0E-02 
2.5E~02 
l.OE-02 
5.0E-03 
2.5E-03 
l.OE-03 
5.0E-04 



6.12E 

4.27E- 
2.88E- 
1.64E- 
1.03E- 
6.35E- 
3.18E- 
1.82E- 



-02(— ) 
-02(0.52) 
-02(0.57) 
-02(0.62) 
-02(0.66) 
-03(0.70) 
03(0.75) 
03(0.80) 



3.34E 
2.59E- 
1.97E- 
1.34E- 
9.75E- 
6.92E- 
4.24E- 
2.85E- 



-01(— ) 
-01(0.37) 
-01(0.39) 
-01(0.42) 
-02(0.46) 
-02(0.49) 
-02(0.53) 
-02(0.58) 



3.04E+00( ) 

2.80E+00(0.12) 
2.54E+00(0.14) 
2.20E+00(0.16) 
1.94E+00(0.18) 
1.70E+00(0.19) 
1.41E+00(0.21) 
1.21E+00(0.22) 



Test 6.2.1b 



l.OE-01 
5.0E-02 
2.5E-02 
1.0E~02 
5.0E-03 
2.5E-03 
l.OE-03 
5.0E-04 



2.84E 
1.87E- 
1.20E- 
6.34E- 
3.78E- 
2.19E- 
1.02E- 
5.56E- 



-02(— ) 
-02(0.60) 
-02(0.64) 
-03(0.70) 
-03(0.75) 
-03(0.79) 
-03(0.83) 
04(0.87) 



1.95E 
1.47E- 
1.08E- 
6.92E- 
4.80E- 
3.24E- 
1.87E- 
1.20E- 



-01(— ) 
-01(0.41) 
-01(0.44) 
-02(0.49) 
-02(0.53) 
-02(0.56) 
-02(0.60) 
02(0.64) 



2.51E+00( ) 

2.27E+00(0.15) 
2.02E+00(0.17) 
1.70E+00(0.19) 
1.47E-f00(0.21) 
1.27E+00(0.22) 
1.03E+00(0.23) 
8.74E-01(0.24) 



behaves similarly to that of the Monge- Ampere equation, that is, we observe the 
following rates of convergence as £ — > 0+: 

\\u-ul\\L2~0{e), \\u-uI\\hi ~0{ei), \\u ~ uI\\h2 ~ 0{ei) . 

Since we have fixed h very small, we expect that ||u — behaves similarly. 



Test 6.2.2. In this test, we calculate using the Hermann- Miyoshi mixed 
finite element method developed in the previous subsection to calculate the rate 
of convergence of — with respect to h for fixed e. We also compare the 
numerical tests with Theorem |6.10 Since u'^ is generally not known, we solve the 
following problem (compare to (6.49|-(6.50)): find ((t|,m^) S W^e x Q^e such that 



(6.63) 
(6.64) 



[al, Kh) + (div(K,i), VO = G{Kh) 



(div(afJ,Vz,,) 



{i + \^ui\^y 



where 



Q^. : = {vh e Q''; 



dn 



Test 6.2.1 



Test 6.2.1a 




10^ 10^ 10"' io' 

e 



Figure 11. Test 6.2.1. Change of \\u - w.r.t. e {h ^ 0.01) 
We use the following test functions and data: 

(a).^ = e^ fl + x[+^fe^?+^ 

' 0.1(1 + (a;2 + a;2)e-?+-i)' 

- £(4(1 + xl + xl) + i2 + xl + xlf)e-^ , 

.g" = , (j)" = (^(1 + xl)uf + 2xiX2i^it/2 + (1 + xj 

JC = 0.1. 

{h)u' ^l{xl + xl)\ /C = 0.1, / = J(x?+a;2)4, 
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Table 5. Test 6.2.2: Error of - u^|| w.r.t. h {e = 0.01) and 
estimated rate of convergence 

h li"' - «hlL2(rato) ||«' - ti^ 11^1 (rate) Her' - cr^|1^2 (rate) 

Test 6.2.2a 2.00E-01 2.04E-04( ) 5.98E-03( ) 4.40E-02( ) 

l.OOE-01 2.60E-05(2.97) 1.52E-03(1.98) 1.68E-02(1.39) 

5.00E-02 3.28E-06(2.98) 3.72E-04(2.03) 6.07E-03(1.46) 

2.50E-02 4.16E-07(2.98) 9.25E-05(2.01) 2.19E-03(1.47) 

1.25E-02 5.24E-08(2.99) 2.31E-05(2.00) 7.87E-04(1.48) 

Test 6.2.2b 2.00E~01 2.05E-03( ) 4.72E-02( ) 3.64E-01( ) 

l.OOE-01 2.77E-04(2.89) 1.19E-02(1.99) 1.46E-01(1.32) 

5.00E-02 3.66E-05(2.92) 2.89E-03(2.04) 5.44E-02(1.42) 

2.50E-02 4.72E-06(2.95) 7.09E-04(2.03) 1.97E-02(1.47) 

1.25E-02 6.02E-07(2.97) 1.76E-04(2.01) 7.04E-03(1.48) 



7(6xlxUxf + xl) + I5xfxj{xf + x^) + 20xlx^ + x}^ + xl^ 

J , s 2 



0.1(^1 + x\(x\ + a;2)6 + xl{xl + x\f 
- 2%%e{x\+x\f, 

0^ = {lx\ + X^^(x\ + X^^'^v'l + \2{x\ + xl^)^ XyX2VxV2 

+ {7xl+xl){xj+xl)\i. 



We record the results in Table [5] and plot the results in Figure [T2j The data 
clearly indicates the following rates of convergence: 



\\u'-ul\\H^=Oih'), \\u^-ulU.=Oih^). 



These are exactly theoretical rates of convergence proved at the beginning of this 
section, indicating that our theoretical estimates for u^ — u| are sharp. On the other 
hand, we note that the numerical rate is better than the theoretical estimate for 
a'^ — af^ which is expected because the theoretical rate of convergence for cr^ — af^ 
is clearly not optimal from the approximation point of view. This phenomenon 
also occurs when approximating the linear biharmonic equation by the Hcrmann- 
Miyoshi finite element method (cf. f35|). 



Test 6.2.3. For this test, we use our numerical method to approximate /C* and 
compare our results with those found in [4] , where the method of continuity (which 
was used to prove the existence of classical solutions to the equation of prescribed 



Gauss curvature) was implemented at the discrete level. We compute ( 6.31 )-( 6.33) 
with the following Dirichlet boundary conditions and domains as used in [4]: 
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Test 6.2.2 




10-^ 

h 

Test 6.2.2 




10-' 

h 

Test 6.2.2 




Figure 12. Test 6.2.2. Change of Wu" - w.r.t. h {e = 0.01) 



(a) 5 = 


V 




■''21 






= 


(-0.57,0.57)2. 




(b)5 = 


1 


-4- 








?t = 


(-0.57,0.57)2. 




(c)g = 


1 


- {xi - 


' 0.075)2 _ 


(2^2 


-0.015)2, 


n = 


(-0.57,0.57)2. 




(d) .9 = 


V 


'l-xl 


- t2 
■''21 






n = 


(-0.72,0.72) X (- 


-0.36,0.36) 


(e)ff = 


1 


~x\- 


•'''2J 






n = 


(-0.72,0.72) X (- 


-0.36,0.36) 


(f) .9 = 


1 


- (a;i - 


' 0.075)2 _ 


(2^2 


-0.015)2, 


n = 


(-0.72,0.72) X (- 


-0.36,0.36) 



Ill 



Table 6. Test 6.2.3. Computed K* with e = -0.001 and h = 0.031 

Computed K.* IC* in |4| 

Test 6.2.3a 2.07 2.10 

Test 6.2.3b 2.20 2.24 

Test 6.2.3c 1.95 1.85 

Test 6.2.3d 2.68 2.61 

Test 6.2. 3e 2.71 2.73 

Test 6.2. 3f 2.20 2.27 



We remark that for the above choice of data, the solution of the prescribed 
Gauss curvature equation is concave, and so we set e < in order to approximate 
the solution (see |37^ I61j for further explanation) . Table [6] compares our results 
and those of [J. Table |6 shows that our numerical method gives comparable values 
to those computed in 1^. Finally, we plot the computed solution of Test 6.3a for 
/C- values 0,1, and 2 in Figure [Tsl We also compute and plot the corresponding 

I — X2 and e = 0.001) for comparison. 



1, and 2 in Figure 13 
convex solution (with g = 



'1 
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Com BUled Solution: K=0,0 Man: 0.B50 Com puled Sol jtiori K=0.0 Max: -0.592 




Figure 13. Test 6.2.3a. Computed concave solution (left) and 
convex solution (right) with JC = 0.0 (top), JC — 1 (middle), and 
/C = 2 (bottom), h = 0.025 and e = —0.001 to compute the 
concave solution, where as h — 0.025 and e — 0.001 to compute 
the convex solution. 



113 



6.3. The infinity-Laplacian equation 

In this section, we consider finite element approximations of tfie infinity-Laplacian 
equation: 



.65) 
.66) 



u = g 



in ri, 
on do,. 



where 



A, 



D^uVu ■ Vu 



1 



d^u du du 



and g £ C(9ri). We note that unlike the PDEs considered in the previous two 
sections, the infinity-Laplacian equation is not fully nonlinear, but rather quasilin- 
ear. Still, its non-divergence form, degeneracy, and strong nonlinearlity in the first 
order derivatives makes the PDE difficult to study and approximate ([H I33L 163) ). 
In particular, the linearization of the operator Aqo gives a degenerate linear dif- 
ferential operator which serves as a perfect example for testing the mixed finite 



element theory developed in Section 5.4 



Remark 6.11. As pointed out in Section 6.1 both AqoW ■— D^vVv ■ Vw and 



|Vl>|2 



are called the infinity-Laplacian in the literature [3l I34j because 



Aoo^ 

they give the same infinity-Laplacian equation. Here we adopt the latter definition 



for a reason which will be clear later (see Remark 6.14). 



The infinity-Laplacian equation (6.65 1 arises from the so-called "absolute mini- 

v\ 



mal problem" which is stated as follows: Given a continuous Junction g 
find a function it : f2 R such that for each V <Z VL and each v G C{V) u^^y — ■^^gy 
implies esssupy|Vu| < esssupy|Vw|. The equation finds applications in image pro- 
cessing and many other fields, we refer the reader to two recent survey papers 
[51 I28j for detailed discussions on the latest developments on PDE analysis and 
applications of the infinity-Laplacian equation. 

Like the equation of prescribed Gauss curvature, we have some flexibility in 
defining F{D^u, Vw, u, x). One possibility is to define F{D^u, Vu, u, x) :— — AqoM, 
but this leads to difficulties in the linearization (see Remark 6.14). Here, we define 



.67) 



F{D^u,Vu,u,x) 



AqqU 

|Vm|2 +7 



|VuP+7 ' 



where 7 > is a positive parameter that will be specified later. The reason for 
introducing 7 is to avoid dividing by zero in the expression. 
It is easy to check that 



r [v\(w) = ,„ ,^ h 2 



|Vt;|2+7 ■ (|V'y|2-H7) 



F'[^i,v]{k,w) = ,_ ,„ ^ h2^ ^ 



2 _1_ -,,^2 ' 

w 



if 
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The vanishing moment approximation becomes 



(6.68) 

(6.69) 
(6.70) 

The hnearization of 
at the solution is 



e|2 



= 



+ 7 

= g 



2„,£ 



in ri, 

on dfl, 
on dfl. 



|V-u^|2 +7 



|VW=|2 +7)2 



Numerical tests indicate that there exists a unique solution to ( 6.68 )-( 6.70) (cf. 
Subsection 6.3.3 and |37j ). and therefore, for the continuation of this section, we 
assume that there exists a unique solution to ( 6.68 )--( 6.70). 

Before formulating and analyzing finite element methods for ( 6.68 )-( 6.70), we 
first state the following two identities. 

Lemma 6.12. Suppose that n = 2. Then there holds the following identity: 
iVwpdAwp - \D'^w\^) ^ {AwVw - D'^wVw) ■ V{\Vw\^). 



The proof of of Lemma 6.12 is a straight-forward (and tedious) calculation, 
so we omit it. Next, with the help of Lemma |6.12[ we are able to establish the 
following identity. 

Lemma 6.13. Suppose that n — 2. Then for any v G Hq{Q,), there holds 



.71) 



{F'[u%v),v) 



,e|2 



L2 



(|Vw=|2+7)2 



Proof. Integrating by parts we get 



|Vu'^|2+7 ' 

Vv ■ Vu" 
iVu-^p +7 



+ 



Vu^P +7)2 



|Vm'^|2+7 ' 



|Vu^|2 +7' 



iVu'^p +7 



(|Vu^|2+7) 



Ar^u'^v 



|Vu'^|2 +7 



|V-U^|2 +7 



Au^v 
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Thus, 

{F'[u^]iv),v) 



— , VW • VU^ - — — ^ ,W + — r:, ^AU^V 



£|2 



Noting that 



L2 
2 



L2 



|Vu^|2+7 ' 
|Vu^|2 +7 



div 



\ |V-u^|2+7 

|Vm'^|2+7 (|Vm^|2+7)2 
_ (|Vu'^|2+7)(|Au^|2 - |d2u^|2) _ (Aw=-Vw^' - i:'2wSVii'^) • V(|Vw^'|2) 



(|Vm^|2+7)2 



we have by Lemma |6.12[ 

{F'[u^]iv),v) ^ 




2„e|2 



7 /|Au^|2 - \D^U' 

2 (|Vu^|2+7)2 

^ U|Vm^|2+7)2' 



□ 



Remark 6.14. (a) UnUkc the two PDEs analyzed in the previous sections, the 
operator F'lu"^] is not uniformly elliptic, that is, there does not exist constants 
Kq, Ki > such that 

{F'[u']iv),v) > KiMhi ^ Ko\\v\\l2 yveH^in). 



Thus, when constructing and analyzing mixed finite element methods for (6.681 



(6.70), we must instead use the abstract analysis of Section 5.4 which is developed 



exactly with such a case in mind. 

(b) If we set F{D^u, Vu, u, x) = — AqoU := D^uVu ■ Vu, then the linearization 
of F would be 

F'[v\{w) ^ -D^wVv ■ Vv - 2D^vVv ■ Vw, 
and it is an easy exercise to see that 

,2 1 



{F'[u^]{v),v) = II Vw • Vu" 



(|Au^|2-|Z?2y^|2,^;2) 
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Thus, the reason we use the definition (6.67) is so that we are able to control the 



zeroth order term in the linearization as shown in the following corollary. Neverthe- 
less, numerical experiments of |37L 161) indicate that the vanishing moment method 
with F(D^u, Vu, It, x) = —AqoU also work well for the infinity-Laplacian equation. 



Corollary 6.15. Suppose n — 2. Then there exists a constant 70 
such that for 7 G (0, 70] , there holds 

(6.72) {G',[u']{v),v)>Cs\\v\\l, V« € I/q- 



70(e) > 0, 



Proof. If UVu-^Hl^ ^ 0, then by ( |6.71[ ), we have 

2 

{G',[u%v),v)>Ce\\v\\l, 



m'^P + 7 



>Ce\\v\\l.~^ 



L2 



(|Vu-|2+7)2'^ 



Vw^P + 7)2 



M 



> C 



,e||2 



71 



Choosing 70 — ^ 11 '^'ii 2"^°° ' '^^ desired result. 

On the other hand, if ||Vm^||l=° — 0, then = const and F'[u'^] = 0, then we 
can choose 70 to be any positive number to obtain 



{G',[u^]{v),v)=e\\Av\ 



2 

L2- 



□ 



6.3.1. Conforming finite element methods for the infinity-Laplacian 
equation. The finite element method for ( 6.68 1-( 6.70) is defined as finding u| G 
Vg^ such that 



(6.73) 



e{Aul,Avh) 



A. 



-,Vh 



dvh \ 
dvj, 



where we assume that 7 S (0, 70] for the rest of this subsection so that the inequality 
(6.72) holds. Furthermore, we assume that ||Vm^||l°° > 1- This assumption is not 
necessary in our analysis, but is does simplify our presentation (cf. ( |6.76| ). 

The goal of this section is to apply the abstract framework of Chapter |4] to the 
finite element method (6.73). Specifically, we now show that conditions [A1]-[A5] 



hold, which will then gives us the existence, uniqueness, and error estimates of the 



solution to (6.73). Of particular interest is the constants' exphcit dependence on e 



in the error estimates. We summarize our findings in the following theorem. 



Theorem 6.16. Suppose n — 2, and let g H'^{VL) he the solution to (6.68)- 
(6.70) with s > 3. Then there exists an = /13(e) > such that for h < /13, (6.73) 



has a unique solution. Furthermore, there holds the following error estimates: 

«-2| 



(6.74) 
(6.75) 

where 

C2 -- 



- u^iWh^ < Cih ~\\u- 



\\u' ~uiU2 <Cs(C2h'\\u'\\He+CjLih)h^'-^\\u'\\Hc), 



'W '9-1 



C7 



Ce-^j-^\u'\w2,a.\u'\ 



w '9-1 



Cg = CC'jCr. 
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where q is a number in the interval (l,oo), Cr is defined by (6.771, L{h) is defined 
by (6.871, £ = min{s,A: + 1}, and k denotes the polynomial degree of the finite 
element space. 



Proof. First, Corollary 6.15 implies that (G'^[u'^])* is an isomorphism from 
VotoV„*. 

Next, we note that 



(6.76) 



(|Vu^|2 + 7)^ 



I 'p—2rn 



< 1, 2m>p>l. 



Thus, for any v,w £ Vq, we have by using Sobolev inequalities for any q € (1, oo) 



{F'[u%v),w) 



< 



|Vu^|2 +7 



\Vv\\l2\\Vw\\l'^ 



,w 



iVu-^p +7 



\Vv\\Lq\\w 



Li- 



|Vm^|2 + 7 
7v\ 

< C\\D^u^ 



\Vv\\Tq\\w 



< C{ \\Vv\\l2\\Vw\\l2 + \\D'u'\\^_^ ||Vw||l.||u;||l= 



Next, by the standard PDE theory, if we assume that u'^ and dft are sufRciently 
smooth, and if v £ Vq solves 

where ip is some L^{n) function, then v e HP{n) for p > 3. Furthermore, in view 
of Remark 4.4 and the inequalities (which come from (6.76l) 









du" du" 
dxi dxj 




< 










|Vm^|2 +7 




|Vu^|2 + 7 



^p^ 



< 2 



iVu-^p +7 



<C{ WD^u'l 



+ 2 



(|Vu^|2+7)2 



iVu-^p +7 



l^viii. 



\Vu' 



e|3 



(|Vu^|2+7)2 



<C||W||i^, 
we have that in the case p ~ 4 

\\v\\H^<Ce-^D'um^ML.. 
It then follows that condition [A2] holds with 
(6.77) Co = Ce, Ci = Ce, C2 = C\\D'^u'\\ _^ 



It then follows from Theorem 14.31 that 
(6.78) C4 = CC2S~\ C5 = CCiCRE-^, 



ho^C{C2CRy 
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To confirm [A3]-[A4], we set 



W 1 



(6.79) 
where q G (1, oo). 

Using a Sobolev inequality and the inequality (6.76), we have for any y G 
Y, v,w€ V„ 



< c 



/ 



IV2/P + 7 



\Vv\\l2\\Vw\\l^ 



\D'y\\^_^^\\^v\\LA\wU^ 



< C-i --\\D^y\\^^Jv\\nA\w\\n^. 

Here, we have used the fact that for p < 2m and a; > 0, 

some constant that only depends on -p and m. 
It then follows from this calculation that 

\F'{y\\\^y. 



p — 2 m 

< C7 2 for 



sup 



WvW 



< c. 



Thus, [A3]-[A4] holds. 



To verify condition [A5] , we first note for any Vh G Vg and w & Vq 



(6.80) 



{F'[u^]^F'[v,])i 

12 



^ V|Vw/ip +7 ^ |Vu^|2 +7 



• Vw - 



|VWft|2 +7 



|Vu' 



e\2 



-7 



^ (|Vu-|2+7)2 (|V«„|2+^)2 



|VW,,|2 +7 



e|2 



2Aooit^Vii'' • Vw 



(|Vu^|2+7)(|Vv,,|2+7) 



(|V?/^|2+7)(|Vz;ft|2+7) 
1 



1 



(|V?/-|2+7)2 (|V«;,|2+7)2 



(|Vi;„|2 + 7)2 

Bounding the second, fourth, and sixth term on the right-hand side of ( |6.80 1, 
we use Sobolev inequalities to conclude that for WT'^u"^ — w/i||ff2 < (5 G (0, |) and 
for any q G (1, 00) 

|V^,|2^|V^^|2 
^(|VW^|2+7)(|V«,,|2+^)^ 

/ (Vt;/,-VM'^)-(Vi;fe+VM") ' 
(|Vu-|2+7)(|V«„|2+7) 



(6.81) 
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Here, we have used that fact that if — w/i.||if2 < (5 £ (0, |) and || Vu'||loo > 1, 

then ||Vwft||ioo > C for some positive constant C that is independent of h, £, and 
7- 

Similarly, 
(6.82) 

|V7;,|2-|Vu-|2 



(|Vu-|2+7)(|Vi;;,|2+7) 



and 
(6.83) 



1 1 

(|VU^|2 +7)2 ~ (|Vv|2 +7)2 



(6.84) 



\ (|Vm^|2+7)2(|Vi;|2+7)2 J 

< CpV||^_^ \\u'\\%2\\u' - Vh\\H4M\H-- 

Bounding the first term in ( |6.80 ), we use similar techniques to conclude 
D^VhVvh - D^u^Vu") ■ Vw 



|VW;,|2 +7 



< C-/-' \\D^Vh - D^u'')Vu' ■ Vw\\li 



To bound the third term in (6.80), we use the identity 

D'^wVvh ■ Vvh - D'^wVu" ■ Vit" = D'^w{Vvh + Vu") • {Vvh - Vu"), 
to obtain 



(6.85) 



< C||M^||/i-2||u'' - Vh\\m\\w\\H2. 



|Vt;;,|2 + 7 

Next, we write 
= (Aoou" - Aoo'y/i)Vv/i • Vw + AooU^(Vu^ - Vw/j) • 

= - D'^Vh)Vvh ■ Vvh + D^u'Wu' ■ Wu' - D'^u'^Vvh ■ Vv,^Vvh ■ Vw 

+ Ar^u^iWu" - Vvh) ■ Vw 
= [{D^u" - D^Vh)Vvh ■ Vvh + D^u''{Vu'' + Vvh) ■ {Vu" - Vvh)^Vvh ■ Vw 

+ K^u'{yu' -Vvh) - Vw, 

so that 
(6.86) 



AooW^Vu^ - AooVfiVvh) ■ Vw 

(|V«„|2+^)2 
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Applying the bounds (6.81)-(6.86) to the identity ( 6.80[ ), we obtain 

{{F'[u-]-F'[v^]){w),z) 



\F'[u']~F'[vh] 



\vv* 



sup sup 

weVo zeVo 



< sup sup 

weVo zeVo 



\\w\\h^\\z\\h^ 
\F'[u^]^F'[vh]){w)\ 



Hence [A5] holds with 
.87) L{h) ^CWD^u'l 



\u^\\UW' 



vh\\m 



for any q G (1, oo). 

Gathering all of our results, existence and uniqueness of a solution to the finite 
element method (6.731 and the error estimates (6.74)~(6.75) follow from Theorem 
Ijland the estimates (|6.77|)-(|6.79|). □ 



6.3.2. Mixed finite element methods for the infinity-Laplacian equa- 
tion. As noted in the previous subsection, is possibly degenerate, and there- 
fore we need to resort to the abstract formulation and analysis of Section |5.4| for 
mixed finite element approximations of the infinity-Laplacian equation. 

The mixed finite element method for ( 6.68 1-( 6.70) is then defined as follows: 
find (ct|,u|) e Wl' X such that 



(6.89) 



where t e (0,to) (tq is defined in Lemma 5.121 

b{Kh,ul) = (div(K/i), VO, 
c{(yl,ul,Zh) = {F{(jl,ul),Zh), 



2 , CT^Vu|-Vu| 



Vu||2 + 7 |Vu||2-H7' 



We also recall that 



on 



The goal of this section is to apply the abstract analysis of Section [5^ to the mixed 
method (6.88 )-( 6.89 1. We summarize our findings in the following theorem. 



Theorem 6.17. Let £ H^^n) be the solution to ( |6.68[ )-( [6?70| a nd let a" = 
D^u^ + rlny.nu'^ With T € {0,Tq), where Tq is defined in Lemma 5.12. Then there 
exists hi — ft-4(e) > such that for h < h^ there exists a unique solution to (6. J 
(6.891. Furthermore, there holds the following error estimates: 

(6.90) 111(5^ - alu' - ul)\\\, < Ksh'-^u^lH^, 



(6.91) 



IWh^ < Ki,jK9h'-'\\u'\\Hi+KlRih)h 



2f-4|L,e||2 



\h' h 



121 



where 



h\\ti\\m + IIa^IIl^ +t^\\v\\hi, 
£ — min{s, k + 1}. 



K. = CK.K 



K'i is defined by (6.93), Kn-^ is defined by (6.92), and Kq is defined by (6.99). 



Proof. First, by (6.76 1 for any v,z € Qo and for any q e (2, oo) 



{F'[a',u']iD'v,v),z) 



(|Vm^|2+7) 



-,w 



< 



\Vv\\l2\\Vw\\l2 



L 9-2 



L '!-2 

From this calculation, we conclude 

ll^'[-^-1llQQ•< ^11-11,^ 

for some q € (2, oo). 

Therefore, using the same arguments as those used in the proof of Theorem 



6.16 



we can conclude that condition [B2] holds with 
(6.92) Ko = Ce, 



-2 



Kh, = Ce 
Next, to confirm [B3]-[B4], we set 



X 



K2=C\\CT-L^., 

p = 4. 



\\{uj,y)\\xy:Y = 1 ^W^W Vlu e X,y eY. 

L 9-2 

where q is any number in the interval (2, oo) We then have for any uj G X, y G 

Y, x<^W, veQ, zeQo, 

{F'[u:,y]ix,v),z) 



< 



/ xVy • Vy — 2wV?/ • V 

(|Vy|2 + 7) 
|Vy| 



z] +2 



llxlU2||z 



L2 



(cjVy ■ y)Vy ■ Vt) 
(|V2/|2+7)2 ' 



(|Vy|2+7)^ 



(iVyP+J) 

I|Vi'||l2||^IL ^II^IIl. 



\\Vv\\l2\\z\\l4u;\\ 2^ 

L 9-2 



L9-2 



<C7-5|lc^|l 2, (||x||L2 + |^'||H0ll^llff^- 

L 9-2 
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It then follows that 

\\F'[uj,y](x,v)\\H-^ < C||(^,y)l|xxy(||xlU^ + I^II^O' 
and thus, assumptions [B3]-[B4] hold with 



(6.93) 



7" 3 \\Tl^a' _ T& 



\XxY 



Next, for ifih,Vh) € W,^ x Q^^ with \\\{n''a' - fih,l''u' - Vh)\\\^ <S e (0, i) 
and {Kh,zh) eW^ xQ'' 



(6.94) {F'[a', u'] - F'[fih,Vh]) {i^h, z^) 



|Vw,,|2+7 |Vw^|2+7 \\yvh\^+-/ |Vw^|2+7 



(|Vu^|2 + 7)2 



To bound the first term in (6.94 1, we add and subtract terms to deduce 
|Vw,,|2 +7 iVw-^P +7 

l^-h^Vh ■ Vf/j — KhS/u^ ■ Vu"^ _ , _ ^ / 1 1 



and therefore by the inverse inequality. 



|Vi;,,|2+7 |Vu'^|2+7^ 
I (|V«^|2+7)(|Vm-|2+7) 



(6.95) 



|Vu,,|2+7 |Vw^|2+7 
< c(||Vw,, + Vu"||l2||Vu,, - Vii^||L2 

+ ||Vu^||i=. ||Vu^ - Vvh\\LA\'^u' + Vi;,,|li2 ) WkuWl^ 



Using a similar technique to bound the second term in (6.94), we first write 



V|Vt;h|2+7 |Vu^|2+7 



^Zh 



Vzh 



|Vt;,J2+^ 
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It then follows that 
(6.96) 



/ Hh'^Vh cr^Vu^ 



< 



— o'^||z,2||Vu/(||i:,oo||V2/(||i,2 



+ ||(T^||L~||Vt;,, - Vu'\\L2\\VZh\\L^ 

+ ||ct^||loo||Vu^||i,oo||Vw^ + Vt;ft||z,2||Vw^ - Vt;^||L2 ||V2:^||loc) 

< C(^h'^\\flh - cr^||L2||u^||Hi + ||cr^||Loo||u'^ - Vh\\H^ 

+ h-^\\(7'\\L^\\Vu'\\L^\\u'\\Hi\\u' - Vh\\Hi)\\zh\\m 

Next, we write 

^ {a'^Vu' ■ Vu'')Vu^ ■ Vzh - {nh^Vh ■ VvhWvh ■ ^ zn 
(|V^;,,|2+7)2 

+ (cT^Vu^ • Vm^)Vu^ • Vzh 



(|V«-|2+7)2 (|Vt;ft|2+7)2 
Noting 

(ct^Vm^ • Vm^)Vu^ • Vzh - {iihSvhSvh)Vvh ■ Vzh 

+ {i^h^Vh ■ S7vh)iVu^ - Vvh) ■ Vzh, 
we conclude 

(cT^Vu^ • Vu^)Vu" • Vzh - ifihVvh ■ Vh)Vvh ■ Vzh 



(6.97) 



(|Vi;;,|2+7)2 

<c(\\a' - fih\\L4'^u'\\l^\\Wzh\\L^ 

+ yh\\L4^u' - VVhhAl^u' + VVh\\L^\\Vu^\\L<^\\VZh\\L'> 

+ Wl^hhAl^^hWlooWVu' - Vi;ft||i2||V^/,||ioc) 
< c(||Vu^||ioc +/i-^||a^||i,2||M^||Hi||Vu^||Loc 

+ h-^\\a'\\L2\\u'\\jji'^\\u' - VhWrnW^hWrn- 

Wc also have 

1 1 



(ct^Vu^ • Vw^)Vu^ • Vzh 



(|Vu-|2+7)2 (|Vt;,,|2+7)2 



(cr^Vu^ • Vu' 



JVM • yZh 12 , .AO/|v7..ffl2 , -A2 



(|Vt;/,|2+7)2(|V«-|2+7)2 
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and therefore, 



(6.98) 



1 



1 



X|Vu^|2+7)2 (|Vz;^|2+-y)2 

<C{\\a'\\L^\\Vu^\U\\Vzh\\L^\\\Vvh\^ 



< Ch-^\\a'\ 



e||3 IL.eilS 



m\\u - VhWH^WZhWm- 



Gombming ( 6.95 )-( 6.98), we have 



3 



\F'[a' ,u'] - F'[^lh,Vh]\\^, < c(^h-^a'\\L^\\Vu'\\l^.\\u'\\ 
It then foUows from the inverse inequahty, that 

{F' [(T^ ^u'']- F'[flh,Vh]){Kh,Zh),Wh 



sup 



< 



C\\ogh\-^(^h-''\\a'\\L^\\Vu'\\l^\\u'fH2 



X (||cr" - fJ-hh^ + Wu" - Vh\\H^){\\l^h\\L^ + \\zh\\H^), 

and therefore conditfon [B5] holds with 



R{h) = C|fog/i|5( 



Finally, we confirm assumption [B6]. First, we note that 



du' du" 
dxi dxj 



iVu^P +7' 



and 



d ldF{a^,u^ 
dxk V drij 



dxidxk dxj dxi dxjdxk n dxi dxj ^ 



e\2 



and therefore by (6.76) 



max 

l<i,j<2 



dF{a^,u'') 



max 

l<i,j<2 



dri 



< c. 



(|Vm^|2 + 7)2 



(II 



2„,e| 



and therefore by Proposition 5.4 condition [B6] holds with 
(6.99) 

Finally, we apply Theorem |5.15| to obtain existence and uniqueness of a solution 
(a^,u|) to the mixed finite element method ( 6.68 )-( 6.70) as well as the estimates 
(|6J0|-(|6JT|). 

□ 



6.3.3. Numerical experiments and rates of convergence. 
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Test 6.3.1. In this test, we numerically solve the infinity-Laplacian equation 
using the Argyris element of degree fc = 5 for fixed h — 0.015 while varying e. The 
purpose of these experiments is to estimate the rate of convergence of ||w — in 
various norms, where u is the viscosity solution of (6.65)-(6.66|. To this end, we 
solve the following finite element method (compare to (6.731): find u| G Vg such 
that 
(6.100) 



e{^ul,^Vh) 



dvh \ 
dvj, 



We set n = (-0.5,0.5)2, ^ 
(a) u 



and use the following two test functions: 



4/3 



4/3 



(b) u = xl +X2, 



/ = 0, 



8(a;2 + y2) 



We note that the second test function is smooth, but the first does not belong 
to C'^{Vt) since its second derivatives have singularities at xi = and X2 = 0. 
After computing the solution for different e-values, we list the errors in Table [7] 
with their estimated rate of convergence and plot the results in Figure [15] The 
numerical experiments indicate the following rates of convergence as e — > 0+ : 



o 



u 



o 



(-). 



Since we have fixed h small, we expect that 



\\u-uI\\h2 « O (es) . 
has similar rates of convergence. 




Max: 6.8586-3 




6.3.2 



Figure 14. Test 6.3. la. Computed solution (left) and its error (right) 
with e = 0.001 and h = 0.015. 

Test 6.3.2. For our last test, we verify the theoretical results derived in Section 
To this end, we solve the following problem: find ((T|,ufJ G Wic such that 



(6.101) 
(6.102) 
where 



K^L Vh) - e~^c[al, ul, Vh) = if, Vh) ^Vh e , 
W^. :- {^^h e fihiy ■ iy\g^ - 0" + rg}. 
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Table 7. Test 6.3.1. Error of - it||| w.r.t e {h ^ 0.015) 



e 




II" - «hlli2(ratc) 


II" - "hllffi (rate) 


II" - "hll/f2 (rate) 


Test 6.3.1a l.OE- 


-03 


2 15E-03C ) 


3 90E-02(' ) 




5.0E~ 


-04 


1 52E-03fO 50) 


3.29E-02(0.25) 




2.5E- 


-04 


1.06E-03(0.52) 


2.75E-02(0.26) 




l.OE- 


-04 


6.54E-04(0.53) 


2.15E-02(0.27) 




5.0E- 


-05 


4.51E-04(0.54) 


1.77E-02(0.28) 




2.5E- 


-05 


3.09E-04(0.54) 


1.45E-02(0.29) 




l.OE- 


-05 


1.88E-04(0.55) 


1.10E-02(0.30) 




Test 6.3.1b l.OE- 


-03 


1.01E-02( ) 


7.20E-02( ) 


1.36E-h00( ) 


5.0E- 


-04 


6.02E-03(0.75) 


5.10E-02(0.50) 


1.21E-H00(0.16) 


2.5E- 


-04 


3.61E~03(0.74) 


3.70E-02(0.46) 


1.08E+00(0.16) 


l.OE- 


-04 


1.86E--03(0.72) 


2.50E-02(0.43) 


9.36E-01(0.15) 


5.0E- 


-05 


1.15E-03(0.70) 


1.89E-02(0.41) 


8.44E-01(0.15) 


2.5E- 


-05 


7.14E~04(0.68) 


1.44E-02(0.39) 


7.63E-01(0.15) 


l.OE- 


-05 


3.87E-04(0.67) 


1.01E-02(0.39) 


6.70E-01(0.14) 



We use the following test function: 

= cos(a;i) — cos(a;2), 0'^ = v\ cos(a;2) — cos(a;i), 

cos(xi) sin^(xi) — cos(a;2) sin^(a;2) 



= e(cos(xi) — cos(x2)) -f 



sin^(a;i) + sin^(a;2) -|- 7 
We compute ( |6.101 )-(6.102 ) for fixed e — 0.01, while varying h with il = 
(—0.5, 0.5)^ and 7 — =lE-4. We list the error of the computed solution in Table 
|8] for both t = and r = 1. As expected, for the case r = 1, we observe the 
following rates of convergence: 

||u^-<|U2 =0(/i3), \\u' -ulU^^Oih"), W -al\\v^^O(h). 

We also observe that the same rates of convergence appear to hold for the case 
T = 0, although our theoretical results of Section 6.3.2 do not cover this case. 

Table 8. Test 6.3.2. Error of - u\\\ w.r.t h {e ^ 0.01) 



11"^ - "hlli2(i-ate) - <|lj:ji(rato) |lcr^ - o-,^|1^2 (rate) 

T = 2.0E-01 8.74E-06( ) 3.92E-04( ) 9.93E-03( ) 

l.OE-01 1.14E-06(2.94) 1.03E-04(1.93) 4.11E-03(1.27) 

5.0E-02 1.38E~07(3.05) 2.53E-05(2.02) 1.54E-03(1.41) 

2.5E-02 1.67E-08(3.05) 6.22E-06(2.03) 5.08E-04(1.60) 

l.OE-02 1.33E-09(2.76) 9.98E-07(2.00) 1.29E-04(1.49) 

T = 1 2.0E-01 8.74E-06( ) 3.92E-04( ) 9.93E-03( ) 

l.OE-01 1.14E-06(2.94) 1.03E-04(1.93) 4.11E-03(1.27) 

5.0E-02 1.38E-07(3.05) 2.53E-05(2.02) 1.54E-03(1.41) 

2.5E-02 1.66E-08(3.05) 6.22E-06(2.03) 5.08E-04(1.60) 

l.OE-02 1.22E-09(2.85) 9.98E-07(2.00) 1.29E-04(1.49) 



Test 6.3.1 




Figure 15. Test 6.3.1. Error \\u~u'jJ\l2 (top), — 
die), and — uf^Wn^ (bottom) w.r.t. e (h — 0.015). 



CHAPTER 7 



Concluding Comments 



In this final chapter, we give some conchiding comments about the vanishing 
moment method and its finite element and mixed finite element approximations for 
fully nonlinear second order PDEs. In particular, we point out some main issues 
accompanying with the methodology. 

We recall that the vanishing moment method and the notion of moment so- 
lutions are exactly in the same spirit as the vanishing viscosity method and the 
original notion of viscosity solutions proposed by M. Crandall and P. L. Lions in 
[24j for the Hamilton- Jacobi equations, which is based on the idea of approximating 
a fully nonlinear PDE by a family of quasilinear higher order PDEs. The vanish- 
ing moment method then allows one to reliably compute the viscosity solutions of 
fully nonlinear second order PDEs, in particular, using Galerkin-type methods and 
existing numerical methods and computer software (with slight modifications), a 
task which had been impracticable before. As a by-product, the vanishing moment 
method reveals some insights for the understanding of viscosity solutions, and the 
notion of moment solutions might also provide a logical and natural generaliza- 
tion/extension for the notion of viscosity solution, especially, in the cases where 
there is no theory or the existing viscosity solution theory fails (e.g. the Monge- 
Ampere equations of hyperbolic type [20j and systems of fully nonlinear second 
order PDEs.) 

7.1. Boundary layers 

As pointed out in Chapter [2] in order to approximate a second order PDE by 
a quasilinear fourth order PDE, we must impose an extra boundary condition such 



as those given in (2.11). Because the extra boundary condition is artificial, it is 
expected that a "boundary layer" ought be introduced in an e- neighbor hood of dQ. 
For example, in the case that Au'^ = e is used as the extra boundary condition on 
dfl, and since we do not know a priori the true value Au on Sfi (note that Au may 
not even exist if the viscosity solution u is not differentiable) , Au^ and Am take 
different values on dft in general, and the discrepancy between Au"^ and Au could 
be large although this can only occur in a very small region (i.e., an e-neighborhood 

of dn). 

Since the convergence of to it as e \, 0+ is only expected and proved in low 
order norms (cf . Chapters [2] and |3| , the error Am^ — Au in an e- neighborhood of 
dfl does not cause any problem for the convergence. Our numerical experiments 
do confirm this conclusion. Moreover, as expected, our numerical experiments also 
confirm that \\u^ — it|j//2 does not converge in general (cf. Test 6.1.3). On the 
other hand, a closer look at the error of computed solution in Figure [5] shows that 
the error is concentrated in an £-neighborhood of dfl and at the singularity of the 
solution u. 
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To improve the accuracy and efficiency of the vanishing moment method, we 
propose the following simple iterative surgical strategy, which consists of three steps. 



Step 1: Solve numerically ( 2.9 )-( 2.11 )i as before for a fixed (small) e > 



Step 2: Find Am^ on the inner boundary of the e-neighborhood of Sfi, and 
extend the function to dft by any (convenient) method. We denote the extended 
function by c^. 

Step 3: Solve numerically (2.9)-( 2Jl] )i again with A-u^|ao — £ being replaced 



by Altaian = 

Remark 7.1. (a) can be obtained by an interpolation technique, or by doing a 
ray tracing along the normal on dVL, or simply by letting be the maximum value 
(a constant) of Aw| on the inner boundary of the e-neighborhood of d^. 

(b) Clearly, Step 2 and Step 3 can be repeated, although one iteration is often 
sufficient in practice (see numerical experiment below). 

(c) The above iterative surgical strategy is a "predictor-corrector" type strategy, 
where the prediction and correction are done on AmU^o. 

(d) To make the algorithm more efficient, the solution computed in Step 1 can 
be used as an initial guess for the nonlinear solver in Step 3. 

As a numerical example for the iterative surgical strategy, we solve the Monge- 
Ampere equation using the conforming finite element method developed and ana- 



lyzed in Section 6.1.1 that is, we numerically solve (2.9 1 with F{D'^u'^, Vu^, u^,x) = 
f{x) — det(D^'u^) using the finite element method (6.6). Here, we use fifth de- 
gree Argyris elements to construct the finite element space, and set f2 = (0, 1)^, 
f = {\ + x\ + X2)e^^~^^^, so that the exact solution is m = e(^i+^2)/2^ 

In Step 2, we extend Aw| in the neighborhood of dO,, to by linear interpo- 
lation to construct c^. After performing Steps 1-3, we repeat Steps 2 and 3 four 
more times to determine whether repeated iterations make a significant impact on 
the error. We use the parameters e — 0.01 and h = 0.01 for all computations. 

After computing the solutions in Step 1 and 3, we record the errors in Tables 
[l||2j We also plot the cross-section of the computed Laplacian Au^ at X2 — 0.8 
in Figure [l] after each iteration. Tables [T]-[2] clearly indicate the iterative surgical 
strategy decreases the error at each step. In fact, the error in every norm is de- 
creased by nearly a factor of ten by performing Steps 1-3 just once. However, the 
error decreases only modestly after repeated iterations and has no impact on the 

and errors after two iterations. Figure [l] also indicates that the boundary 
layer is greatly reduced after the first iteration, and improves modestly after each 
subsequent iteration. 

iteration # |ju-Mf||i2 Wu-uIWh^ 



1.48E-02 l.OOE-01 1.79E-H00 

1 1.88E-03 2.11E-02 4.63E-01 

2 1.51E-03 1.23E-02 2.53E-01 

3 2.15E-03 1.18E-02 1.77E~01 

4 2.51E-03 1.24E-02 1.42E-01 

Table 1. Errors of u — m| using the iterative surgical strategy 
(e = 0.01, = 0.01). 
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iteration # — ||u — wfjl 11^* — ^illln'^."" 



2.02E-02 4.25E-01 2.93E+01 

1 3.94E-03 9.50E-02 6.14E+00 

2 3.52E-03 5.06E-02 3.83E+00 

3 4.66E-03 3.93E-02 2.79E+00 

4 5.21E-03 3.31E-02 2.50E+00 



Table 2. Pointwisc errors of u — u| using the iterative surgical 
strategy (e = 0.01, = 0.01). 



i - 




Figure 1. Cross-section plot of Am| at X2 — 0.8. Black solid line 
is exact solution, and dotted blue lines are the computed solutions 
for iterations 0, 1, 2, 3, and 4. 



7.2. Nonlinear solvers 



After problem (2.9)-(|2.11 )i is discretized, we obtain the (strong) nonlinear 



algebraic system ( |4.3[ ) or ( |5.8| )-(5.9 1 or ( 5.57 H( 5.58 1 to solve. To this end, one has 
to use one or another iterative methods to do the job. In all numerical experiments 
given in Chapter [6j we use the ILU preconditioned Newton iterative method as our 
nonlinear solver. Since Newton's method often requires an accurate starting value 
to ensure convergence, hence generating a good starting value for Newton's method 
is also an important issue here. So far we have used two strategies for the purpose 
in our numerical experiments in |37L 1381 139] and in Chapter |6] The first strategy 
is to use a fixed point iteration to generate a starting value for Newton's method. 
However, this strategy may not always work although its success rate is pretty high. 
The second strategy, which is more involved, is the following "multi-resolution" or 
"homotopy" strategy: first compute a numerical solution using a relatively large 
e, then use the computed solution as a starting value for the Newton method at 
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a finer resolution e. The process may need to be iterated in e for more tlian one 
step. Our experiences tell that 1 — 3 steps should be enough to generate a good 
starting value for Newton's method at the finest resolution e at which one wants 
to compute a solution. 

It is expected that for 3-d simulations and for time-dependent fully nonlinear 
PDEs (see Section 7.3 below), more efhcient fast solvers are required. It is well- 
known that the key to this is to use better preconditioners for the linear problem 
inside each Newton iteration because solving (large) linear systems inside each New- 
ton iteration costs most of the total CPU time for executing the Newton's method. 
One plausible approach, which will be pursued in a future work, is to use more 
sophisticated multigrid or Schwarz (or domain decomposition) preconditioners (cf. 
|74j ) to replace the ILU preconditioner. With help of the better preconditioners, 
Krylov subspace methods |71J can be employed as the linear solver inside each 
Newton iteration. Put all pieces together, we arrive at a global nonlinear itera- 
tive solver which can be called the Newton-Schwarz/Multigrid-Krylov method (cf. 

])■ 



7.3. Open problems 

As the vanishing moment method was introduced very recently, there are many 
open questions concerning with the method. The foremost one is to generalize 
the convergence results of Chapter [3] to the general problem (2.9)-(2.11) under 
some reasonable structure conditions on the nonlinear differential operator F. The 
convergence rate is probably hard to get unless the viscosity solution of the limiting 



problem (2.7)-(2.8) is sufficiently regular (cf. Theorem 3.191 



Another interesting but completely open problem is to develop a vanishing mo- 
ment method for fully nonlinear second order parabolic PDEs. Unlike the situation 
for quasilinear PDEs, going from fully nonlinear second order elliptic PDEs to fully 
nonlinear second order parabolic PDEs is far from straightforward. One reason 
for this is that there are several different legitimate parabolic generalizations for 

[Ml ES [76]). 



equation (2.7) (cf. 
parabolic PDEs are 

(7.1) 
(7.2) 



Two best known fully nonlinear second order 



F{D^u, Vu, u, X, t) — ut = 0, 
-utdet(L'^u) = f{Vu,u,x,t) > 0. 



Extensive viscosity solution theories have been developed for both equations (cf. 
[4511551 F75l 176) and the references therein). However, to the best of our knowledge, 
no numerical work has been reported for these equations in the literature. 

Formulation of the vanishing moment method for ( |7.1[ ) is straightforward (see 
|37j ). By adopting the method of lines approach, generalizations of the finite el- 
ement and mixed finite element methods of Chapter [4| and [5| should be standard. 
However, the convergence analysis of any implicit scheme is expected to be hard, in 
particular, establishing error estimates which depend on polynomially instead of 
exponentially will be very challenging. Furthermore, we note that numerically solv- 



ing equation (7.2 1 using the vanishing moment method is expected to be difficult. 



In fact, it is not clear how to formulate the method for (7.2). 



Finally, another interesting open question is to explore the feasibility of ex- 
tending the notion of moment solutions and the vanishing moment method to 
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degenerate, non-elliptic, and systems of fully nonlinear second order PDEs (cf. 
[IIl[IIl|40l[20l[58]). 
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